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Compressibility Effects in Two-Phase Flow 


R. F. TANGREN, C. H. DopGe, AND H. S. SEIFERT 
Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California 


(Received September 13, 1948) 


A gas-water mixture, when expanded through a de Laval nozzle, acts as a compressible fluid. 
A general expression for the flow process in an idealized mixture is developed. The general thermo- 
hydrodynamic relations for flow through a de Laval nozzle are developed. The magnitude of the 
velocity of an energy pulse or signal through the mixture is determined and is compared with the 
velocity of flow. Experimental results are compared with those predicted using the equations de- 


veloped. 





I. INTRODUCTION 


HE behavior of a two-phase system, such as a 

mixture of small gas bubbles in water, has 
attracted interest for a considerable number of 
years. As early as 1930 Ackeret' derived an ex- 
pression for the velocity of sound in a gas-water 
mixture, and the dynamics of such systems has 
been studied by Heinrich? (1942), Charyk* (1943), 
and Zwicky* (1944). References to publications of 
these men, as well as more recent workers in the 
field, are given in the bibliography, Section IX of 
this paper. Such a gas-water system is of interest 
because of the possibility of applying it to the 
propulsion of underwater devices in which a gas is 
injected into water, the working fluid. 

The procedure used in the analysis to follow is 
to apply the basic laws of continuity, momentum, 
energy, and ideal gas equations of state to a mixture 
in such a way as to derive: 


(a) an equation of state for the mixture, 
(b) an equation of motion for the mixture. 


1J. Ackeret, “Experimental and theoretical investigations 
on Cavitation in water’’ (Translation from Technische Mechanik 
— acces (Forschung), Vol. 1, No. 1, Berlin, 

0). 

2G. Heinrich, ‘“‘Uber Stromungen von Schaumen,” Zeits. f. 
angew. Math. Mech. 22, 117-118 (1942). 

3 J. Charyk, Jet Propulsion Laboratory Progress Report 2-2 
(California Institute of Technology, November, 1943). 

*F. Zwicky, Aerojet Engineering Corporation Report R-10 
(January, 1944). 


These equations are then examined for the special 
cases in which the mixture approaches: 


(a) the all gas state, 

(b) the all liquid state, 

(c) an intermediate state corresponding to a 
practical experimental mass ratio. 


The flow of the ‘“‘practical’’ mixture through a 
de Laval nozzle is calculated and plotted and some 
comparisons are made with experimental measure- 
ments. 

The present analysis is greatly simplified by the 
assumptions tabulated below. The results, there- 
fore, are only indicative of what might follow from 
a more detailed consideration involving viscosity, 
turbulence, heat transfer, surface tension, vapor 
pressure, resonance, and the like. Surprisingly good 
correlation is found in spite of the simplicity of 
these assumptions. It is assumed that: 

1. The liquid is an incompressible fluid, and 
effects due to its viscosity, surface tension, and 
vapor pressure are unimportant. 

2. The gas is an ideal gas, with negligible vis- 
cosity, constant specific heats, and is insoluble in 
the liquid. 

3. The mixture is “homogeneous” in the sense 
that the bubbles of gas are so small and uniformly 
distributed that an arbitrarily small sample con- 
tains the same mass ratio of gas to liquid as the 
whole mixture. 

4. The gas and liquid are always at the same 
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8+ 2.0 N= p/p* 


yu* p/p 
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Fic. 1. Dimensionless graph of supersonic nozzle flow and 
thermodynamic parameters, referred to values at throat as 
functions of the relative area and for the particular value of 
59=0.50. 


temperature, and the flow is insulated (i.e., adia- 
batic). 

5. The flow is one-dimensional, laminar, and 
expands or compresses slowly enough that inertial 
transients may be ignored. 

6. A slowly varying pressure change or signal is 
transmitted through the mixture with a definite 
critical or “‘sonic’’ velocity, as distinguished from 
the anomalous dispersion and attenuation of high 
frequency pressure waves. 

The nomenclature is defined here in alphabetical 
order for reference. Symbols will be redefined in 
the text at the places where they first appear. 


Roman Letters 


a—sonic or critical velocity, ft./sec. 
A—area, ft.*. 
Cys—specific heat of liquid, ft.?/sec.*°R. 
C,—specific heat of gas at constant pressure, ft.2/sec.?°R or 
ft. Ib./slug°R. 
C.—specific heat of gas at constant volume, ft.?/sec.?°R. 
M;--mass of liquid, slugs per ft.* of mixture. 
M,—Mass of gas, slugs per ft.* of mixture. 
p—pressure, lb. /ft.*. 
Q—heat added per slug, ft.?/sec.*. 
R—gas constant, ft.?/sec.?°R. 
T—temperature, °R. 
u-—velocity, ft./sec. 
’—work done per slug, ft.?/sec.?. 
( )o—value under initial conditions. 
()*—conditions at the throat. 


Greek Letters 


y—tatio of specific heat of gas at constant pressure to the 
specific heat of gas at constant volume. 

é—ratio of volume of gas to volume of liquid per unit 
volume of mixture. 


638 


u—ratio of mass of gas to mass of liquid per unit volume of 
mixture. 
p—mass density of mixture, slugs/ft.’. 
ps—fluid mass density, slugs/ft.*. 
Po—gas mass density, slugs/ft.*. 


II. EQUATION OF STATE FOR A TWO-PHASE FLUID 


A unit volume of mixture at any instant is 
occupied by a mass M, of gas and M, of liquid, 
which are in the mass ratio 


u=M,/ My. 


Since the volume of gas plus the volume of liquid 
must equal the total volume, this fact may be 
used to define an average mixture density as follows: 


M,/ po t+ M;/p;=(Ma+My;)/p. (1) 
As the mixture expands or contracts, » remains 
fixed, but the volume ratio 
5=pupys/ pg 


varies inversely with the gas density. By com- 
bining the definitions of u and 6 with Eq. (1) above, 
one arrives at two useful relations among the 
densities 


p/ py =(1+n)/(1+8) (2) 
and 
i+u 1 
1/0 (1/0)| —-—| (3) 
p Ps 


The equation for conservation of energy, when 
written in terms of energy per unit volume, with the 
assumption that gas and liquid temperatures are 


$8 < 2.00 —— 
$>200—-— |] / fy 


8, + 0. sf 8,70.25 


8, = 0.05 8, =0.50 


= 1.00 


8, * 2.00 


msM" 


8, = 0.05 


8, 8, = 0.10 
8, = 0.25 





° O O02 04 06 O08 10 O08 O06 O4 O02 
CONVERGENT = DIVERGENT 
A/A 


Fd 


- 


Fic. 2. Relation between Mach number ratio (M/M*) and 
area ratio (A*/A) for various values of do. 
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always equal, takes the form 


(M,+M;)dQ= (M,C.+M;C;)dT 
+(M,+M,)dW (4) 
or 


dQ=[(uC.+C;)/(i+n) dT +dW. (5) 


The work done in expanding a unit mass of the 
fluid against the external pressure p is dW = pd(1/p) 
and if this is substituted in Eq. (5), together with 
the fact that in an insulated flow dQ=0, we have 


(p/p*)dp=([(uC.+C,)/(1+u) ]dT. (6) 


The equation of state for the gas component 
alone is 


b/po=RT. (7) 


One makes use of the equation of state (7) and 
the density relation (3) to eliminate the unwanted 
variables p, and T and to arrive at the differential 
equation of state involving only pressure p and 
mixture density p as variables. Inserting (3) in (7), 
differentiating, and combining the result with (6) 
yields the adiabatic state equation 


dp/p= [(uC.+ Cy+uR)/(uCy+ Cs) 


1 
| do/o( _-—- ~)}. (8) 
1+u py 


which may be integrated (noting that C, = C,+R) to 
PL(1/p) —1/ps(1 +m) JeCrt er (“ert Cr = constant. (9) 
If the new and significant parameter is defined 
T= (uCp+Cy;)/(uC.+Cy) 


and Eq. (9) simplified with [ and the density 


$<2.00 


8>2.00——— =2.00 


7 

v8" 1.00 
/ 20.50 
ou 


0.25 
20.10 


20.05 
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Fic. 3. Relation between velocity ratio (u/u*) and area ratio 
(A*/A) for various values of do. 
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Fic. 4. Relation between pressure ratio (p/p*) and area ratio 
(A*/A) for various values of do. 


equation (3), then 
PLu/(1+)p, ]* =constant. (10) 


To express the adiabatic state equation in terms 
of temperature, insert (7) in Eq. (10), from which 


pPL{u/(i+u)}{RT/p}]*=constant (11) 


(12) 


It may be seen from Eq. (12) that I has a signifi- 
cance for adiabatic mixed flow similar to that of y in 
adiabatic gas flow. 

From the definition of T as p> 


r-C,/C.=y7. 


Thus the limiting condition in which ny, (pure 
gas) gives the familiar adiabatic equation for a gas 


pr-)!1/T =constant. 


or 
p-)/T /T =constant. 


On the other hand, as p—->0 (pure liquid) [1.0 and 
the process is isothermal with T =constant. 


Ill. EQUATION OF MOTION FOR A 
TWO-PHASE FLUID 


By applying to Newton’s equation of motion the 
equation of state just derived (Eq. (9)), there 
results an expression stating bulk fluid velocity u in 
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terms of;pressure and known fixed quantities. The 
assumption of no viscosity and one-dimensional 
flow permits the use of the following simple form 
of Euler’s equation : 


udu = —dp/ p. (13) 


Eliminating 1/p by means of Eq. (9) and integrating 
between the limits u, p and uo, po gives the general 
result 


u?=uy?+(20/(0—1) [u/(1 +4) J[1/ poo | 
X [po— p(po/p)"" ]+2(po— p)/ps(Atu). (14) 


Equation (14) may be expressed in terms of 
volume ratio 59 (where the subscript 0 refers to the 
initial point of the flow) by using the relation, 
derivable from Eq. (10), that 


(p/ po)" © = 60/5. (15) 
Hence, we have finally 


u?=Uy"* +2 po pp(1+y){ (Pr (I'—1) ]do 
XL1—(p/po)?-?'P]+(1—p/po)}. (16) 


It is of interest to examine the general relation 
(16) for the cases of pure gas and pure liquid, re- 
spectively. By definition, 59 =yp;/poo. If we use this 
and allow u—-~, then (16) becomes 


47 = Uy? +2(po/ poo) [¥/(¥ = 1) ] 
X[1—(p/po)-”/7], (17) 


8570.05 
20,10 

70.25 
Lg, =0.50 
21.00 

8, 22.00 


8 < 2.00 —— 
$ >2.00 —— 





8+ 0.25 
8,2 0.10 
8, 0.05 
© 9 02 04 06 08 10 08 06 04 a2 0 
}»——CONVERGENT ——DIVERGENT 
7B 
- 


Fic. 5, Relation between volume ratio (6/5*) and area ratio 
(A*/A) for various values of do. 
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which is the familiar equation for the adiabatic 


- expansion of a pure gas. 


If we let w—-0, Eq. (16) reduces to the well- 
known Bernoulli equation for an incompressible 


fluid 
u? = Uy? + 2(po/p;s)[1— p/ po]. (18) 


IV. FLOW OF A TYPICAL MIXTURE 


In order to continue the analysis of mixture flow 
in a quantitative manner, a particular practical 
case will be considered. This makes it possible to 
take advantage of certain simplifying approxima- 
tions. The example now to be discussed is the 
water-helium system, accelerating from rest. 

In order that the gas bubbles be able to do work 
during expansion on the liquid matrix in which 
they are imbedded, the ratio of gas to liquid volume 
must be limited to some maximum value, or the 
liquid will become dispersed in such a way that it 
no longer encloses the gas. The volume ratio 6 has 
been arbitrarily limited to 2.0. It should be noted 
that this is the maximum value of 6 and that 6p is 
usually smaller than this. 

Using water as the liquid and helium as the gas, 
with 6=2.0, typical values of » and I for several 
pressures are calculated below: 


p PF r 
1 atmos. 4.4x10-* 1.00022 
2 atmos. 8.9 10-4 1.00044 
3 atmos. 17.8107‘ 1.00088 


At the end of Section II it was shown that when 
nis very small, ! approaches unity and the mixture 
behaves as if the gas process were isothermal. It is 
seen from the above table that this is the case for a 
practical gas-water mixture. Physically this means 
that the mass and specific heat of the liquid are so 
large relative to the gas that the gas may remain in 
thermal equilibrium with the liquid during expan- 
sion without appreciably lowering the liquid tem- 
perature. This is certainly a close approximation so 
long as 6<2.0. 

In view of the fact that [T—-1.0 and p—-0, Eq. 
(16) may be somewhat simplified. The expression 
(p/po)-/? can be expanded in an exponential 
series, of which only the first two terms need be 
used. 


(p/po)?-Y/F=1+[(("—1)/T] Inp/po. 
Neglecting y relative to unity in (16) gives 


u* = U9? — 2(po/ ps) 5o(Inp/ po) + 2po/ ps1 bases b/ po] 


to a close approximation. 
Assuming expansion from rest, u%=0, we may 
write a dimensionless form of the velocity equation, 


psu?/2po= — bo(Inp/po) +1 — p/ po. (19) 
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Because of the approximations involved, Eq. (19) 
indicates that as the free stream pressure p ap- 
proaches zero, the velocity « approaches infinity. 
This would be true if the liquid acted as an infinite 
heat source and maintained the gas expansion 
strictly isothermal. In an actual physical case, 
§ would exceed the limiting value of 2.0 long before 
the heat of the liquid was perceptibly depleted. 

It is interesting to note in Eq. (19) that the only 
property of the gas which enters into the calculation 
of velocity is its volume ratio 69. In other words, 
the molecular weight, specific heat, etc., do not 
affect u in this approximate case, and a given 
volume of one species of ideal gas is as effective as 
any other. 


V. VELOCITY OF A PRESSURE WAVE 
THROUGH A MIXTURE 


A basic assumption, made at the beginning of 
this analysis, was that a slowly varying pressure 
pulse is transmitted through a mixture with a 
perfectly definite “‘sonic’’ velocity. If this be so, 
then it is possible to apply the results of de Laval 
nozzle theory and show that the bulk velocity of 
the fluid equals this ‘‘sonic’”’ velocity at the throat 
of a nozzle through which a mixture is flowing, that 
is, the Mach number is unity at the throat. 

It is realized that ordinary sound waves or rapidly 
fluctuating pressures are transmitted through bub- 
bles with a complicated array of dispersion, attenu- 
ation, scattering, and other phenomena. It is 
assumed here that the pressure variations existing 
in nozzle flow, unlike high frequency sound waves, 
are transmitted as if the medium were essentially 
homogeneous. 

It can be shown! that the basic equation for the 
velocity a of a wave or signal through an isotropic 
homogeneous elastic medium is 


a’=dp/dp. 


The term dp/dp may be evaluated from the pre- 
viously derived equation of state, Eq. (8). By 
operating on Eq. (8) with Eqs. (2) and (3) and the 
definitions of 6 and 69, one may obtain the ex- 
pression for the speed of a signal through a mix- 
ture as 


a?=(Up/p,)(u/(i+u)J1+1/6}. (20) 
Equation (20) may be verified for extreme values 
of wu. As po 

a’—yp/p,= RT, 
which is the equation for the adiabatic speed of 
sound in the gas. As u—-0, a2, or the speed of 


sound in an incompressible fluid. 
One may write an approximate form of Eq. (20) 


®A. B. Wood, A Textbook of Sound (G. Bell and Sons, 
London, 1930), p. 51. 
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Fic. 6. Relation between density ratio (p/p*) and area ratio 
A*/A) for various values of 0. 


for the practical case where u—-0 and 1.0 as 
follows: 


a? = (po5o/ps)(1+1/6]*. (21) 
With the help of Eq. (15), in which the approxi- 


mation [1.0 has also been used, Eq. (21) may 
be put into the same form as Eq. (19): 


pj0?/2po={60/2)(1+p/podo }’. (22) 
VI. FLOW THROUGH A DE LAVAL NOZZLE 


If the mixture flow now be considered to be 
confined to a duct of varying cross section, the 
equation of continuity requires that 


dA/A+dp/p+du/u=0. 


If the duct has a minimum crpss section or throat, 
as in the de Laval nozzle, at that point 


dA/A=0. 


Combining the continuity equation (with the throat 
condition) with Euler’s Eq. (13), one has the 
basic result that the velocity of flow at the throat 
must be 


u?=dp/dp. 


This, however, is just the magnitude of the ‘“‘sonic’”’ 
velocity as stated in a preceding section. It is thus 
concluded that the flow velocity at the throat equals 
the local signal velocity. This ‘critical’ velocity is 
denoted by a superscript *, and at the throat 


“u=a=4* =a", 


In order to evaluate as a function of 59 the critical 
pressure ratio p*/po, at which sonic flow velocity 
is achieved, the approximate expression (22) for a? 
is equated to the expression (19) for u*. Since 
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Fic. 7. Relation between velocity of a pressure wave (a) 
and 6 with constant pressure (p). 


u*? = q*?, 


(50/2)(1+p*/pod0)? = (1—p*/po) — 50 Inp*/po. (23) 


Equation (23) may be solved graphically to give 
the pressure ratio corresponding to ‘‘sonic’’ velocity 
as a function of the initial volume ratio 69. This 
equation is limited in accuracy by the approxima- 
tion u—0 and r—1.0. 


VII. GRAPHICAL PRESENTATION OF NOZZLE 
FLOW PARAMETERS 


The equations derived in the preceding sections 
are sufficient to supply data on the variation of flow 
velocity, sonic velocity, area, density, gas volume, 
and pressure throughout a de Laval nozzle. Most of 
the parameters have been expressed in terms of 
pressure ratio p/po as the independent variable. 
It will be more instructive and convenient, how- 
ever, to plot them as a function of the nozzle area 
ratio A/A* as the independent variable. This is 
done by calculating A/A* as a function of p/p» and 
plotting all quantities against area rather than 
pressure. 

In dimensionless form the quantities of interest 
are (note that * refers to throat conditions) : 


Local Mach number ratio... M/M*=u/a. 
Velocity relative to throat. u/u*. 


Pressure ratio............ p/p*. 
Volume ratio............. 5/6*. 
Density ratio............ p/p*. 
Eh dibewskicimcy A/A*. 


These quantities are plotted together in Fig. 1 
for a particular value of 59 and separately for a 
range of values of 59 in Figs. 2-6. The sequence of 
calculations necessary to arrive at these curves are 
now described in some detail. The numerical values 
have been tabulated in Table I. 


A. Local Mach Number u/a 


The ratio of u/a may be determined by com- 
puting u/(2p0/p;)* and a/(2p0/p;)! from Eqs. (19) 
and (22). 
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This ratio is the local Mach number M, and 
since the throat Mach number M*=1, M=M/M* 


also. It is computed for a series of values of p/po 
and do. See Fig. 2. 


B. Velocity Relative to Throat Velocity u/u* 


Since u=u*=a at the throat, expressions (19) 
and (22) for u and a may be equated to determine 
graphically or numerically the pressure ratio p*/po 
which exists at the throat; p*/po is, of course, 
a function of 59. We use this value of p*/po in (19) 
to find u*/(2p0/p;)', which is then divided into the 
previously determined u/(2p0/p;)'. u/u* is shown 


in Fig. 3. 


C. Pressure, Volume, and Density Ratios 


Pressure relative to the throat, p/p*, may be 
simply found once the critical pressure ratio p*/po 
is known, since p/p* = (p/po)(po/p*). The volume 
ratio relative to that at the throat, 6/4*, is approxi- 
mately the reciprocal of p/p*, as may be seen from 
Eq. (15). The density relative to that at the throat 
may be found by using Eq. (2) with uw neglected 
relative to unity to obtain 


p/p* = (1+6*)/(1+54), 
p/ p* =(1+650/(b*/po) 1/1 + 50/(b/po) J- 


D. Area Ratio A/A* 


The area relative to that of the throat may be 
calculated from the density and velocity ratios 
using the equation of continuity 


Apu=A*p*u*, 
written in the form 
A/A*=(p*/p)(u*/u). 


The values of A/A* as a function of p/po are corre- 
lated with the preceding variables M, u, p, 6, and p, 
which are then plotted against area ratio in Figs. 2-6. 
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Fic. 8. Relation between u* and pressure (po) 
for various values of do. 
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Fic. 9. Convergent nozzle with flow of water alone. 


E. Variation of “Sonic” Speed a with Volume 
Ratio 5 and Local Pressure p 


Using Eq. (15), the relation (21) for a* may be 
written in the approximate form, valid for 6 positive 
but less than some small integer: 


a* = (p5/p,)[1+1/6 }?. 


It is interesting to note that if the local pressure p 
is held constant, a* has a minimum at 6=1, which 
can be shown by differentiating the above equation 
with respect to 6. This minimum is surprisingly low. 
In Fig. 7, a is plotted against 6 for p=1 atmos. and 
p=2 atmos., and has a minimum value of 66 ft./sec. 
at 6=1 for p=1 atmos. 

It is suggested® that the reason for the much 
lower sonic speed of the mixture when compared 
with that in either the gas or the liquid is that the 
gas acts as a weak spring when coupled with a 
large mass, the liquid. This point can be illustrated 
by considering the speed of a signal through a pipe 
containing a liquid-gas mixture. This mixture must 
be considered homogeneous for this simplified 
analysis which cannot consider the effects of dis- 
continuities in the medium. The speed of sound 
through such a tube may be written as a?=E/p. 
Where E is the coefficient of elasticity, and, for 
the mixture in the tube, E becomes 


E=-—dp/(dV/V). 
Since for an isothermal gas dV,/V,= —dp/p, then 


d V/ V= —= [6/(1 +5) |(dp/p) 
and 


E=((1+6)/6]p. 
If the mass of the gas is neglected, p becomes 
p=p;/(1+4) 


* Anderson, Rush, Jr., and McClellan, M.S. Thesis in Aero- 
dynamics, California Institute of Technology, June 1, 1947. 
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Fic. 10. Convergent nozzle with flow of air-water mixture 
showing expansion after leaving nozzle. 


and the sonic speed squared is 
a? =(p/p;)(1+6)?/6. 

Using Eq. (15), considering the process as iso- 
thermal (['=1), the sonic speed squared may be 
written as 

a? = (podo/ps)(1+1/8)?. 
This equation is identical to Eq. (21) derived earlier. 
Thus it may be seen that this model, combining the 
elasticity of the gas and the density of the liquid, 
may be used to obtain a physical picture of the 
compressibility effects in two-phase flow. 

The specific value of the sonic velocity at the 
throat, a*, is of interest because the actual flow 
velocity u* equals a* at the throat. In Fig. 8, 
u*=a* has been plotted for a range of values of 
the quantities pp and 49, which are the parameters 
which are most easily measured and controlled. 
It should be noted that for any given po, u* becomes 
a minimum at 6 9=0.10, which corresponds to 
6* = 1.0. 





Fic. 11. Mixture flow with convergent-divergent nozzle. 
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VIII. COMPARISON OF THEORY AND EXPERIMENT 


A limited amount of experimental work was 
carried out at the Jet Propulsion Laboratory, 
C. 1. T., on the flow of gas-water mixtures through 
ducts. The object of this work was to study the 
thrust produced by such a mixture as it expanded 
through a nozzle. 

The gas was injected in a “‘mixing”’ section where 
the duct cross section was large enough that the 
liquid velocity was negligible. Figure 9 shows the 
test equipment with water alone flowing through 
it and Fig. 10 shows the same equipment with gas 
being injected. Since this duct had only a con- 
tracting section, considerable expansion was evi- 
dent after the mixture left the duct at the point 
of narrowest cross section. In Fig. 11 a short 
expanding section had been added to the duct. 

So long as incompressible liquid only is flowing 


' through the duct into the atmosphere, the boundary 


conditions require that the static pressure at the 
exit of the duct be atmospheric. If gas bubbles are 
added, the mixture becomes compressible and above 
the critical flow speed the medium becomes in- 
capable of transmitting pressure changes upstream 
against the flow. At this point all pressures inside 
the duct are proportional to ‘‘reservoir”’ pressure po 
and.are no longer influenced by the atmospheric, 
or discharge, pressure. As a consequence, pressures 
higher than atmospheric may be obtained at the 
exit section of the nozzle, and expansion to atmos- 
pheric pressure may occur after leaving the nozzle, 
as is implied in Fig. 10. This figure is therefore 
taken as qualitative evidence that compressible 
flow at the critical speed is occurring. 

A series of runs was made in which increasing 
amounts of gas were added to a constant volume 
flow of water, thus increasing both 69 and po. 
Although it was desired to measure po/p*, the 
chamber to throat pressure ratio, it was actually 
necessary to measure po/p1, where p; was the 
pressure at a point slightly upstream from the 
minimum cross section such that A*/A, was 0.993. 
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TABLE II. Comparison of measured and calculated 
pressure ratios. 











Percent 

Calculated disagree- 
Measure values values ment 

Po pi 
abs. abs. 

bo (p.s.i.)  (p.s.i.) (po/p1) (po/p1) (po/p1) 
0.193 29.7 14.7 2.02 2.33 13.3 
0.270 32.7 15.7 2.08 2.15 3.3 
0.492 34.7 19.7 1.80 2.00 10.0 
0.84 40.7 23.7 1.71 1.93 11.4 
1.36 48.7 29.7 1.66 1.73 4.0 
1.67 6.0 


2.01 54.7 34.7 £57 





This shift in location was necessary because it was 
not possible to locate a pressure tap precisely at 
the end of the converging nozzle. 

Since the static pressure is very sensitive to 
area ratio in the vicinity of the throat, the accuracy 
of the experimental measurement is limited. In 
Table II the measured and calculated values of 
po/p. are compared for a series of measurements 
made with the pressure tap slightly upstream from 
the throat (A*/A,=0.993). 

Perhaps it should be pointed out that the experi- 
mental work, which was directed toward other 
objectives than the verification of the theory, was 
terminated before this present analysis was carried 
out; thus no opportunity to investigate discrep- 
ancies has been available. 
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Fluid Flow Patternst 
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A new photographic technique is described, for observing velocity fields in water, induced by the 
entry of solid missiles. The results of a number of observations are given, and compared with theo- 
retical formulas. The corrections involved in the method are analyzed theoretically, and the theory 
of the corrections confirmed experimentally. It appears that the method is well adapted to velocity 
fields 5-100 feet per second, created by objects several inches or more in diameter, but may lead to 
troublesome corrections at other speeds or with smaller objects. 





I. EXPERIMENTAL TECHNIQUE 
1. Slichter Multiflash Bubble Technique 


HE velocity field of a moving liquid may be 

observed through time-exposed photographs 
of illuminated particles! suspended in the liquid, 
or on its surface. This method is often called 
“chronophotography.”’ 

We wish to report on a new technique for ob- 
serving transient velocity fields in water, with a 
vertical plane of symmetry, such as occur in water 
entry and cavity motion problems. In this tech- 
nique, which is due to Professor L. B. Slichter,’ 
multiflash pictures are used in place of time ex- 
posures; this permits a greater concentration of 
light energy. Moreover, the particles photographed 
are located in a vertical plane by the following 
device. 

A thick-walled tube of soft rubber, perhaps 0.5 in. 
in diameter, is punctured at regular intervals by a 
wire needle about 0.01 in. in diameter. When the 
tube is inflated under water at an air pressure of 
10-20 lb./in.?, air bubbles of diameter d=0.2-2.0 
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Fic. 1. 


"+ We report here work done under Task B of Contract 
NOrd-8555 with the Bureau of Ordnance. We wish to thank 
Professors L. J. Hooper and F. S. Findlayson of the Alden 
Hydraulic Laboratory of the Worcester Polytechnic Institute, 
where the experiments were performed, for their advice and 
cooperation. 

1 See Bibliography at end of this paper: [5], [6], and Et 
For a general discussion of photographic techniques, see [2 
pp. 280-94, and [4], Vol. 2, pp. 265-74, and Plates 1-27. 

2 In collaboration with the Morris Dam Group at the Cali- 
fornia Institute of Technology, working under OSRD con- 
tract. 
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mm are emitted by each hole. These rise under 
gravity, in a predetermined, approximately vertical 
plane. (See Fig. 1.) If a symmetric missile moves 
through the water in this plane, the bubbles will 
also move in it, giving a plane section of three- 
dimensional motion. 

If the position of a bubble is located by two light 
flashes at a time Af apart, and Ax is the vector 
displacement of the bubble as measured directly 
(by simple proportion) on a photographic plate, 
then the average vector velocity will be 


v=Ax/At. (1) 


A sample photograph is shown on the cover. 

The complete experimental set-up is shown in 
horizontal plan in Fig. 2. When the light source is 
about 30° to the rear* of the bubble plane, the light 
is reflected from a small spot on the bubble, whose 
diameter is often as little as one-third of the bubble 
diameter d. The center of the spot can then be 
located with an error less than d/12, so that Ax 
can be measured to about 0.06 mm in favorable 
cases.‘ 

However, we should remember that what we are 
measuring with 0.06-mm accuracy is the displace- 
ment of bubbles in a finite container, whereas we 
wish to infer the corresponding displacement of an 
infinite extent of water. The “‘corrections’’ needed 
to make this inference are discussed in Sections 9-13. 


Glass Walled Tank 





Flash Tubes (Edgerton) 
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000 
Air Bubbles 

















* This was also essential for obtaining enough illumination 
and is incompatible with the use of a plane of light to select 
a plane of bubbles. 

* An accuracy of 0.025 mm is claimed in [4], using a pre- 
pared emulsion, but the droplets were not confined to a plane. 
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Fic. 3. Vertical cavity due to 8 X2-in. prolate spheroid. Model speed 29 ft./sec. at entry. Time interval about 0.003 sec. 


2. Experimental Difficulties 


One of the main experimental difficulties comes 
in obtaining a uniform array of “good’’ columns of 
regularly spaced bubbles. Bubbles more than 0.9 
mm in diameter spiral up irregularly, and thus do 
not stay in a predetermined plane. Bubbles much 
less than 0.4 mm in diameter rise so slowly as to 
drift out of the plane with small convection cur- 
rents. 

Most synthetic rubbers we punctured gave ir- 
regular individual holes. Commonly, a continuous 
stream of air would spurt out, and form irregu- 
larly into bubbles of various sizes. Each individual 
hole, carefully punctured in natural soft rubber, 
gave a “‘good” column at a suitable air pressure, 
consisting of bubbles of uniform diameter, spaced 
perhaps 2 cm apart. At higher air pressures, bub- 
bles would be too large; at lower pressures, they 
would be too small or not form at all. 

The most difficult trick was to make all the bub- 
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ble columns ‘“‘good’’ at the same air pressure. Best 
results in this regard were obtained by puncturing 
the tube when wet, inflated with air about 1/5 
atmosphere overpressure and stretched about 20 
percent. The tube was later used in the same state 
under about 1-atmosphere overpressure, with the 
holes about 120° from the top of the tube. By 
twisting the tube slightly, more uniformity between 
bubble columns could be achieved—larger bubbles 
being formed by any hole at angles @>120° with 
the top of the tube. 

Two General Electric FT-2 1 flash tubes with 
2000 volts and 9 microfarads were used as light 
sources, about 30° to the rear of the bubble plane. 
A separate tube was used for each flash (of the 
double flash), and the light was concentrated on 
the bubble plane by individual cylindrical lenses. 
The flash duration was of the order of 10~5 second. 
Since we photographed a field about 24’ x18”, 
we were limited by the total available light inten- 
sity. We obtained best results with Super Panchro 
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Fic. 4. Deep seal of vertical cavity due to 4 X2-in. oblate spheroid. Model speed 29 ft./sec. at entry. 
Time interval 0.0025 sec. 


Press Sports Type film; Super XX film was also 
satisfactory. A high contrast developer should be 
used. Best results were obtained with Von-L, and 
D-72 1:1 was satisfactory. 


3. Use of Oil Drops 


In order to reduce the corrections for virtual 
mass and ‘‘bubble rise’’ discussed in Part III, we 
tried replacing air bubbles by droplets of liquids 
whose specific gravity o was nearer unity. We found 
(cf. Section 2) that when o>0.9, such droplets 
rose so slowly as to drift out of the vertical plane. 
We had best results with hexane, for which o =0.69. 

The technique used for air (i.e., punctured rub- 
ber tubing under pressure) was not satisfactory for 
oil droplets, because the oil clogs the holes. When 
holes large enough to remain open were used, the 
resulting flow of hexane was too great to be useful. 
Moreover, when a metal pipe with 0.002-in.— 
0.02-in. holes was used, the fluid would cling to 
the sides of the pipe until an enormous bubble was 
formed. Although some study of fluids, detergents, 
and metals was made, we were unable to get good 
results with this technique. 
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The technique used most successfully was the 
excretion of hexane from small hypodermic needles. 
For this purpose #27 hypodermic needles (i.d. 
0.0075 in., o.d. 0.016 in.) were used. The needles 
were mounted vertically on observation tubes 
which were in turn mounted or a brass tube by 
means of short links of rubber tubing for shock 
absorption. An excess pressure of about 50 dynes/ 
cm? was enough to make bubbles of about 2.4-mm 
diameter rise 3 cm apart. Less pressure would form 
the same size bubbles but spaced farther apart. 
Greater pressure would produce a jet from the 
needle. This would form into bubbles less than 1 
mm in diameter and would be difficult to retain in a 
plane. 

The tip of the needle was left intact. A bubble 
would form around the opening in the side of the 
tip, rise vertically to the tip, and then pull free. 
Thus the effect of the surface tension between the 
hexane and the metal was greatly reduced. 

The needles clogged readily if the pressure was 
taken off and water allowed to enter the system. 
The system was cleaned by an internal blast of 
air and the small wire cleaners that are furnished 
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with the needles. However, if air got into the hexane 
line while in operation, the needles would tend to 
clog. 


II. DISCUSSION OF OBSERVATIONS 
4. Vertical Entry of Spheroids 


We applied the experimental technique described 
in Sections 1 and 2 to the vertical entry into water 
of various spheres and spheroids at about 29 ft./sec. 
entry speed, the unstable vertical orientation being 
maintained by threading a taut piano wire through 
the axis of the model. 

A typical velocity field is displayed in Fig. 3. It 
will be noted that the motion of the cavity wall is 
mainly radial, or perpendicular to the cavity axis. 
This was predicted by Mallock® and Bauer. It 
would also occur in an ideal fluid if the cavity wall 
were a “free boundary,” and the spheroid moving 
at constant velocity—the velocity being y/2 for- 
ward of the radial direction, where the free bound- 
ary makes an angle y with the direction of motion. 
In Plate 1 the boundary layer seems to be percep- 
tible. Just outside of this, the radial component of 
the fhuid velocity seems about one-fourth the ve- 
locity v of the spheroid, and the rearward com- 
ponent perhaps v/10, near the point of separation 
of flow. 

As the spheroid descends further, the water is 
accelerated inwards by gravity, so that the cavity 
collapses on the axis in a “‘deep seal,’’ described by 
Mallock. Fig. 4 illustrates the velocity field just 
after deep seal; a thin jet is rising inside the cavity 
along its axis. 


5. Vertical Entry of Wedges 


We made studies of a 90° wedge similar to those 
for spheroids of Section 4. We did this because 
there was known® a complete theory of the steady- 
state cavity (i.e., wakes) flow past a symmetric 
90° wedge in a nonviscous fluid. For a wedge with 
verticals at (1,0), (0,0), and (0,1), the complex 
position coordinate z is related to the complex 
conjugate velocity vector ¢, by the formula 


4% 3 
1.7442 = ———_ + —_ 
(i+g*)? 1+¢4 
1 t+Ww c+W 1 
+- (Wein + Win +—, (2) 
4 ‘-W c—-W*s 2v2 


where W=(1=72)/vV2=(—1)*. In Fig. 5 we have 





5A. Mallock, Proc. Roy. Soc. A95, 138-43, (1918); W 
Bauer, Ann. d. Physik 82, No. 7, (1927). 

63), pp. 104-5, formula (27), from which Eq. (2) can be 
deduced. Flat plates are mathematically simpler, but de- 


celerate too rapidly under free fall to approximate steady- 
state conditions at all. 
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Fic. 5. Comparison between experimental and theoretical 
flow around a wedge. Experimental data taken from Plate 35; 
4-in.-90° wedge, 2-ft. long. Theoretical flow normalized to 
produce experimental water motion at bottom of wedge. 


compared observed velocities with steady-state 
theory. The observed velocities are 20 percent-— 
40 percent less than the predicted velocities, and 
are directed 0°-20° more towards the surface. 

We conjecture that both these effects are due to 
the presence of the horizontal free surface, but the 
second may be due to the correction of Section 14. 


6. Two-Dimensional Jet 


The two-dimensional cavity behind a 4-in. diam. 
wedge entering at 29 ft./sec., collapses in a ‘‘deep 
seal,’ much like the axially symmetric cavity 
behind a spheroid (cf. Section 4). We display in 
Figs. 6 and 7 (see cover) both the pinching off of the 
cavity before deep seal, and the approximate ve- 
locity field of the up-jet, which is formed after deep 
seal above the point of deep seal. Thus, in Fig. 6, 
the water is moving inwards towards the cavity 
(which appears light gray), and also downward in 





Fic. 6. Deep seal of 90° 2-ft.x4-in. wedge. Model speed 25 
ft./sec. at entry. Time, interval 0.0025 sec. 
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Fic. 7. Deep seal of 90° 2-ft.x4-in. wedge. Model speed 
25 ft./sec. at entry. Time interval 0.0025 sec. Shows flow 
accompanying up-jet formed at surface. 


the upper portion of the picture. Other photographs 
reveal the velocity field associated with the down- 
jet formed below the point of deep seal. 

We were able to make a mathematical model of 
an infinite two-dimensional jet with a free surface 
as follows. Suppose the velocity § at infinity is 1, 
we have a physical z-plane, a complex velocity po- 
rential W=U+iV, and a conjugate velocity 
¢=dW/dz. One can map a complex half ¢-plane 
onto the shaded region of the W-plane (see Fig. 8) 
by the generalized Schwarz-Christoffel transforma- 
tion 


dW =tdz=(tdt/t+1). (3) 


One can map the complex half of the r-plane onto 
the shaded area of the ¢-plane (hodograph plane) 
by setting 27=¢+1/t¢. Then 4(0)=0, 7r(0)=@—; 
t((B)=—1, 7(b)=1; t(A)=e@, 7r(A)=-—1. These 
points correspond under the (unique) linear frac- 
tional Moebius transformation, t= —2/(7r+1). Sub- 
stituting, we get ¢+1 = (¢—1)?/(¢+1)’, dt=4(¢—1)/ 
(§+1)8, tdt/(t+1) = —16¢/(¢ —1)(¢ +1), and so 





dz= —16/(§—1)(¢+1)*. (4) 
This can be integrated in closed form, as 
4 4 t+1 
2=— - +2 In——+C. (5) 
(¢+1)? ¢+1 ¢—1 


In Fig. 9 we have compared the flow pattern (20) 
with Fig. 7; the agreement is surprisingly good, in 
view of the facts that gravity is not negligible, 
that the flow is not stationary, and that we have a 
horizontal free surface. 
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7. Oblique Entry of Spheres 


We also studied the velocity fields associated with 
1.5-3.0-in. diam. spheres entering water obliquely, 
at speeds of 80-150 ft./sec. 

The instantaneous streamlines in the vertical 
plane of motion are approximately concentric 
circles, centered in the point of impact. 

In Fig. 10 we have compared the vector bubble 


displacements, as observed in one of our photo- | 


graphs, with the field of a uniform line of sources, 
and reflected line of sinks, as sketched in Fig. 11. 
(The graphical construction of this field is greatly 
facilitated by use of Ex. 3 on p. 10 of Kellogg's 
Potential Theory.) 

It will be noted that the velocity field is far from 
being ‘‘axially symmetric’ about the axis of mo- 
tion—greater above the cavity than below it. 







Experimental Cavity 


Theoretical Cavity 


/ ———->Bubble Motion from 

Plate 39 Corrected for Bub 

ble Rise and Wall Effect. 
Theoretical Two- 

Dimensional Flow with In- 
finite Jet Normalized to 
V4 Have Same Magnitude as 
Bubble Motion at Point * 1°. 








4 


Fic. 9. Comparison between experimental and theoretical 
z-dimensional jet flow. 
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Theoretical Water Motion due to Line of Sources and 
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Point 1) Source Line is Taken from Surface to Mean 
Position of Center of Sphere. 





Fic. 10. 


This asymmetry of the velocity field is consistent 
with the observed fact’? that the surface effect 
on hydrodynamical forces is small, if a sphere is 
four diameters below the surface. A sink m diameters 
above the surface would cause velocity fields at 
points »—} and m+} diameters below the surface, 
differing only in the ratio 1+ $n. If we consider 
that most of the velocity is due to sources below the 
surface and compute pressure variation from 
Bernoulli’s equation (3) with dU/dt=0, we are led 
to expect an unbalanced force of at most 1 percent 
at four sphere diameters submergence. 


8. Surface Seal and Ricochet 


It is known*® that at high Froude numbers 
F =y"/gd>5000 or so, the ‘‘deep seal’’ of Section 4 
is replaced by a “surface seal’’ which may be de- 
scribed somewhat as follows. The splash is blown 
inward by air rushing into the cavity. This forms a 
constriction or “‘neck.”’ The air rushing through 
this neck, causing a local underpressure, is }p,v*, 
where p; is the air density, and v is at least the 
missile velocity. This airflow probably separates 
from the cavity wall at the edge of the splash, 
leaving a turbulent air-filled or foamy region at 
mixed pressure between the inrushing air and the 
cavity wall. Some water is thrown against the 
cavity wall, which it makes turbulent and opaque. 
The underpressure causes the neck to contract 
more and more, choking off air from the cavity, 
until a complete ‘‘surface seal”’ is formed. 

*This was observed by C. Ramsauer, Uber den Ricochet- 
schuss, Thesis (Kiel, 1920). The rough computation below 


seems to us clearer than Ramsauer’s “pressure zone”’ ex- 
planation. 


SSee [11], [12]; also Mallock (reference 5). Roughly, 
surface seal supersedes deep seal when 4$9,v* exceeds the hydro- 
Static overpressure 3pgd at the depth (3 diameters) of deep 
seal—or when F>6p/p:=4400, roughly. 
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Line of Sources 


We were fortunate in photographing the ap- 
proximate velocity field in the neck of the cavity. 
In Fig. 12, the rear side of the neck is seen to be 
moving inwards horizontally; this we found to be 
true in all cases. The front (upper) side of the cavity, 
on the other hand, usually moved upwards at first, 
but later reversed its direction of motion, at least 
with sufficiently oblique entry. The ‘‘dipole”’ field 
of secondary flow, clearly revealed in Fig. 13, is just 
beginning to form under the cavity. 

More recent photographs show that, in the case 
of shallow entry angles (35° or less with the hori- 
zontal), after surface seal is complete, a dipole-like 
velocity field travels vertically downward, away 
from the cavity. This is illustrated in Fig. 13, and 
is presumably due to the down-jet accompanying 
surface seal. This down-jet is especially visible in 
photographs taken at steep entry angles, and can 
be distinguished in these cases from the indenta- 
tions of the cavity wall caused by inward-blown 
splash. 

Another interesting series of pictures concerns 
the ricochet of a 3’’ aluminum sphere at 7° impact 
angle (see Fig. 14). It is remarkable that the water 
velocity tends to be tangential to the cavity at the 
end of the long plough, and not radial, as in most 
cavity motion. It is true that this is necessary if the 
circulation around the closed curve of Fig. 15 is 
zero (as required by Kelvin’s theorem for non- 
viscous fluids), if the flow along the free surface is 








Fic. 12. Surface seal behind 1.5-in. diam. sphere, which 
entered at 150 ft./sec., with trajectory 30° below horizontal. 
Time interval 0.004 sec. 
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Fic. 13. ‘‘Dipole’’ moving down from point of surface seal, behind sphere which entered at 250 ft. /sec., 


with trajectory 


15° below horizontal direction. Time interval 0.004 sec. 


mostly normal to the surface, and if the interior 
streamlines escape to the free surface, as they do. In 
Figs. 14 and 15, the missile has moved out of the 
photographic field to the right. 


III. CORRECTIONS 
9. General Remarks 


The observations described above are subject 
to three corrections. 

The first is the ‘‘wall effect’’ due to the finite 
size of the container. We estimated it on the basis 
of the Dirichlet flow around a 4-in. diameter sphere 
in a 32-in. diameter spherical tank, as given in [3], 
Section 93, pp. 124-5. We did this because we 
worked at least 16 in. from the wall of a 52 52-in., 
34-in. deep tank. The estimated effect seemed small. 

The other two involve the relative motion of the 
bubble in the water, due to the difference between 
bubble density and water density. This difference 
causes the bubbles to rise under their own buoyancy. 
The rate of rise » of bubbles in each column can 
be measured; it is normally between 3 and 8 
cm/sec., and consistent to within +10 percent in 
each “good” column. We simply subtracted a 
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vertical “bubble rise correction’’ of v,At from each 
observed displacement. This was usually small 
(cf. the bubbles near the edge of each photograph) ; 
however, because of it, the multiflash bubble tech- 
nique for observing velocity fields is only useful for 


fluid velocities above 30 cm/sec. 





Fic. 14. Ricochet of 3-in. aluminum sphere. Entry angle 8°, 
velocity 115 ft./sec. Time interval 0.0018 sec. Sphere has left 


the picture on the right. 
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The difference between bubble density and water 
density also makes the bubbles respond more 
violently to pressure gradients (i.e., accelerations) 
than the surrounding water. The resulting ‘‘virtual 
mass” correction is the main difficulty with the 
multiflash bubble technique; we shall now discuss 
it directly. It should be observed that the problem 
has arisen before, in connection with the theory 
of suspensions.? However, the treatment sketched 
below, and first given in [13], offers some novelty. 





10. “Virtual Mass’”’ Correction 


We estimated the “virtual mass’ correction, 
which might be expected in the absence of gravity 
(buoyancy) and walls, by assuming equilibrium 
between the following three forces on each bubble. 

Force 1. A pressure thrust ma, where m was the 
mass of water displaced by the bubble, and a 
was the acceleration of the surrounding water. 

Force 2. An inertia “‘force’’ of moa,+3m(a,;—a), 
where o was the specific gravity of the bubble, and 
a, was the bubble acceleration. Here me is ordinary 
inertia, and 3m is the “induced mass’’ inertia. 
({3], p. 124). 

Force 3. A resistance D to the relative bubble 
velocity vi—v(v; being bubble velocity, and v the 
velocity of the nearby water). We estimated this 
from empirical data, for steady relative motion. A 
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reasonably good fit (cf. Fig. 16) for an undistorted 
air bubble is given by 


Cp =0.4+24/R, (6) 


for R<10°. Here R=|v—v,|d/v is the Reynolds 
number. For distorted bubbles, we used (7) 
(below) and 

Cp=0.4A+24/R. (6’) 


A is the ratio of actual cross-section area to undis- 
torted cross-section area. 

There is no theoretical justification for assuming 
that Forces 1-3 are independent, but we could 
think of no better assumption. If the assumption is 
made, the condition of equilibrium becomes 


m(o+})ai—1.5m=Da=}p(v1—v)?ACp. (7) 


where A = zr’ is the cross-section area of the bubble, 
and Cp is the drag coefficient. 


11. Case of Impulsive Velocity: Bubbles 


In case the water is given an impulsive velocity v, 
which is thereafter held constant, it is easy to 
deduce from Section 10 that the initial relative 
velocity v:—v of an air bubble should be V =2v, as 
in Section 3. This relative velocity should then die 
out according to the differential equation dV /dt 
=aV*+B8V, where the constants a, 8 are determined 
by formulas (7)—(6’) with a=0 in (7). This is 
easily integrated, giving a predicted ratio s;/s of 
bubble displacement s; to water displacement s, 
for each time interval At and sphere density o. 
Two curves obtained in this way are shown in 
Figs. 17-18. 
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bubbles nolonger spherical. 





100 6.8 
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*A. B. Bassett, A Treatise of Hydrodynamics (Cambridge University Press, London, 1888), Vol. 2, Chapter V; C. W. Oseen, 
Hydrodynamik (Leipzig, 1927), p. 132; T. Chan-Mou, Mededelung 51 Lab. Aer. Hysr. Delft; G. I. Taylor and M. Glauert, 
Reps, and Mem. 2025 Aer. Res. Comm; E. Brun and M. M. Vasseur, Belleveue No. 3 (1947); Irving Langmuir and K. B. 
Blodgett, Office Tech. Services, U. S. Dept. Commerce, No. 27, 565. 
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We checked this prediction experimentally as 
follows. Suppose a half spheroid (ellipsoid of revolu- 
tion), is floating with its equatorial plane on the 
surface of the water. Let the hemispheroid be given 
a downward impulsive velocity.’ 

Since the pressure on the water surface is con- 
stant, the pressure gradient will be normal to it. 
Hence, so will the acceleration. Integrating, the 
impulsive velocity imparted to the hemisphere will 
also be vertical; this was also confirmed experi- 
mentally by our photographs. We conclude that the 
velocity potential U satisfies 


U = Uo=const. on the (x, z)-plane initially. (8) 


But by the well-known principle of reflection, this 
determines mathematically a unique solution of (1), 


given zero velocity at infinity, and the normal 
velocity 9U/dn on the boundary. This solution 
corresponds (cf. reference 9) to supplementing the 
half spheroid by its mirror image in the water 
surface—i.e., by taking the corresponding flow 
induced by a spheroid, relative to fluid at rest at 
infinity. In the case of ideal fluids, this flow is 
well known ([{3], p. 141), and reduces to a dipole 
field in the special case of a hemisphere. 

In this special case, the theoretical angle 4, 
between the vertical and the direction of fluid 
motion at a point P, subtending a polar angle ¢ 
(see Fig. 19 where 0 is the mean sphere center), 
satisfies 0:,=¢+tan-'3 tang. In Fig. 20 we have 
plotted as a function of ¢, the difference 0.—4 
between @, and the experimental flow direction @., 
corrected for bubble rise and wall effect. This is 
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© This was accomplished by threading the model on a 0.040” vertical steel wire under tension and dropping a 10-Ib. weight 


down the wire. The acceleration of the model took less than 0.0002 sec. and its velocity was about 25 ft./sec. We had trouble 
with wire rupture and horizontal vibration until we stabilized the weight and placed a 1/16-in. shock absorber below the weight. 
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one of the first" quantitative comparisons of pre- 
dicted velocity fields with photographed velocity 
fields. It will be seen that there is a scatter of +5°, 
but that a definite disagreement is present for ob- 
served points near the surface (¢ nearly 90°). This 
may be due to the finite displacement. If we re- 
member that the actual water surface is above the 
mean sphere center, and that theoretically, at all 
times (see Section 12), 6=180° at the surface, we 
are led to plot also 6.—6,* against ¢, where @,* is 
the angle between the vertical and the velocity 
vector as P for flow around a sphere with center at 
the surface. 

In view of Sections 3-4, the ratio s./s of the 
photographed bubble displacement s, (corrected 
for bubble rise) to the theoretical water displace- 
ment s, ‘should vary between one and three. In 
Fig. 17, we have compared the theoretical curves 
described above, with observed values of s,/s, 
plotted against s.. Note the data tends to lie above 
the theoretical curve for s.<0.15 and below for 
s->0.15. This may be explained by the fact that 
the distortion of the air bubbles is not constant but 
greater for areas which are nearer the spheres, and 
thus have greater motion. 

On the other hand, the evident general agree- 
ment between the predicted curve and the ob- 
served points confirmed our view regarding the 
nature and magnitude of the ‘‘virtual mass”’ cor- 
rection. 


12. Impulsive Velocity: Check with Oil Drops 


Analogous experiments were made using hexane 
bubbles. In Fig. 20b, 6.—6; is plotted against ¢ for 
four spheroids. There is a dispersion of about +4° 
in most regions, presumably due to errors in ob- 
servation. As in the case of air bubbles (Section 7) 
one notes a discrepancy between the observed and 
theoretical data near the surface, although below 
the surface (6 <90°) this is not so noticeable. 

As in Fig. 17, we plot in Fig. 18 the resulting 
“theoretical” curves of s:/s and the experimentally 
observed ratios s./s for the case of the four spheroids 
with hexane bubbles. Here we took A in (6’) to be 
2.5 (i.e., the distorted cross section was 2.5 times 
the undistorted area). Again the scatter is of the 
order of +10 percent, but the general agreement 
between the predicted curve and the observed point 
is apparent. As in Fig. 17, one notes that for high 
values of s,, the data lie below the theoretical 
curve. This may be explained as before in terms of 
the distortion of the bubbles. 

We note that in the case of hexane, the ratio of 
s;/s never exceeds 1.26. If the specific gravity of 

1 None is given in [5], [6], or [7]. We have not found one 
of Prandtl, Plate 1(a). A comparison of Pitot tube observa- 


tions with theory has been made by Bryant and Williams, 
Phil. Trans. 225A, 199-245 (1925). 
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the “‘improved’”’ fluid were 0.85, this ratio would 


never exceed” 1.1. 


The case of a disk entering water should be the 
limiting case of an oblate spheroid. Figure 21 is a 
plot of 6.—6, against the polar angle ¢ from the 
mean center of the disk, the motion of the disk 
being less than 0.1d. The discrepancy near ¢=90° 
is more noticeable here than in any other case 
(cf. Figs. 20a—20b). This is probably because the 
water pushed out radially prevents the motion on 
the surface near the disk from being vertical. 

A comparison of the ratio of the magnitudes is 
made in Fig. 22. It is noted that in this one case, 
all the data lie above the theoretical curve and in 
fact, whereas the theoretical curve has as its average 
value about 1.24, the experimental data have an 
average value of 1.45. We are not able to explain 
this. Though the experimental errors were larger 
with oil drops than with air bubbles, we believe 
they were insufficient to explain the preceding 
discrepancy. 


13. Case of Impact 


If an object with principal radii of curvature r, r’ 
at the point of impact is dropped into still water 
with speed v, the rate at which water volume is 
displaced per unit time (i.e., the cross-section area 
at the surface level) is given at first by 2(rr’)*v*t, 
plus terms of higher order in time ¢ after impact. 
It seems reasonable to suppose (cf. Section 4) 
that at large distances, the motion of the water 
corresponds to a dipole with proportional strength. 
Therefore at a fixed point of the water, the in- 
stantaneous velocity should be fixed in direction, 
but have magnitude proportional to t—i.e., the 
vector acceleration a should be approximately constant. 
In this case, one can integrate (7) to get 


t 1 eee ] ) 
" 2kovet+ki—A kit+A 
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12 The factor 3 for air has the unsuspected advantage that it 
magnifies motions far from the missile which might otherwise 
be difficult to observe. 
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where A =[k:?+4(1—o)akz }!, ki = 18vd?, ko =3/10d, 
a is the magnitude of a and d is the diameter of the 
bubble. 

By means of (9), we may plot the bubble velocity 
against time for each value of a. These curves are 
shown in Fig. 23. 

We compared the results of the preceding calcu- 
lations with multiflash pictures of the velocity 
field of a sphere entering water vertically, pre- 
viously taken by Professor Slichter. The compari- 
son of theory with experiment is shown in Fig. 23, 
points associated with the same bubble being con- 
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nected with the same value of a. Again we have good 
qualitative agreement. 

It is to be noted that if there were no relative 
bubble motion, the curves (9) of Fig. 23 would be 
replaced by straight lines. The large discrepancy 
brings out the fact, not previously realized, that 
the multiflash bubble technique is unsuitable for 
studying velocity magnitudes just after impact. 


14. Circular motion and general conclusions 
In case the water moves in a circle of radius R, 


_at constant angular velocity W, the correction 
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predicted by (6’)-(7) can again be found. The 
bubble will move in another circle, also with 
angular velocity W, so that a,;—a will be 90° out 
of phase with v,—v, and will be W times as great. 
Hence, Forces 1, 3 of Section 10 will form the sides 
of a right triangle, whose hypotenuse will be Force 2. 
To a first approximation, this leads to the predic- 
tion that im the case of circular motion, the phase of 
the bubble velocity will be (8r/3CpR)} radians ahead 
of the water velocity vector. Figure 24 gives con- 
tours of this phase advance (in degrees) against 
r/R, for various values of Cp. 

In the case of vertical entry of 1.5-in. diam. 
spheres at 150 ft./sec., this predicted phase advance 
can attain a maximum value as high as 40° or 
more." These maximum values are attained for 
bubbles at some distance from the cavity, and 
therefore having orbits whose radii are small. 
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The predicted phase advances can, in some 
cases, be checked against observation. On many 
double-flash pictures, two cavity outlines are super- 
imposed. Where the outlines cross, the water dis- 
placement is tangential. The measured angle be- 
tween this tangential direction, and the observed 
direction of bubble displacement, agreed very well 
in all cases with the predicted phase advance. This 
agreement confirms further our correction formulas. 

The predicted “virtual mass’ corrections are 
considerably smaller for 4-in. diam. missiles than 
for those of 1.5-in. diam.—perhaps of the order of 
10-20 percent of the velocity near the cavity and 
in front of the missile, and 15-35 percent further to 
the side. In general, pictures taken by the multi- 
flash bubble technique tend to exaggerate motion 
in the direction of acceleration. If Cp=0.4, 99 per- 
cent of the relative velocity is lost in eight (bubble) 
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'8 The measurements and computations underlying these statements were made by Mr. Robert Isaacs. 


VOLUME 20, JULY, 1949 


657 








350 the relative displacement is limited absolutely by a 





“Bubble Displacements few millimeters, except insofar as the instantaneous 

from Multiflash Picture of & | r / : f h ; d 

Sphere Entering Water. F acceleration a cm/sec." 0 the water tends to 
wo Theoretical Curves produce a relative velocity of the order of (v1—») 

Based Upon Constant @ f /9\4 y 

of Water, co*5 + 24 / ~(a/2)* cm/sec. 


APPENDIX 






250° 
y During the course of our work, we made some observations 
Pe on the behavior of bubbles, which may be worth recording. 
* Our measurements of the rate of rise of air bubbles in water, 
200° in the region 40<R<200, gave values of the Cp which were 
2 considerably lower than might be inferred from the standard 
- references (Handbuch. Exp. Physik, Vol. 4, Part II). We dis- 
> covered that a similar disagreement was found in recent inde- 
1SO* = pendent measurements by Bryn and Zdonik,“ and have 
< plotted the data in Fig. 16. (For R>50, the bubbles do not 


rise in a straight line, but spiral upward; for R>1000, they 
are no longer spherical.) 

However, these data must be examined rather critically. 
There is a large mutual influence of two spheres. At very low 
R, Stimson and Jeffery (Proc. Roy. Soc. All1, 110 (1926)), 
have deduced a correction factor 3d/4D, where d is the bubble 
diameter, and D the distance apart. This is corroborated 
visually by watching the apparent “‘attraction’’ of bubbles in 
the same column, and has been estimated empirically as 10 


100° 


(Assumed Diameter of Bubble - .7mm) percent by Mr. Jackson of H.M.U.D.E., Portland. 
o i 2 4 4 5 (6 cf 6 9g We also investigated the ability of bubbles to maintain 
Time - m sec. their spherical shape in a pressure field. In order to have the 
Fic. 23. change in hydrostatic pressure through one bubble diameter d 


be less than the difference between the pressure inside and 


. : : outside the bubble, caused by surface tension, we must have 
diameters of bubble motion, if the water undergoes g < 30 cm; with hexane drops, we get the same result. 


no acceleration in this time. Hence, in practice, Under an acceleration of Ng, the corresponding figure would 
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Fic. 24. 
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4“ T. Bryn, Forschung 4, 27 (1933); Zdonik, Absorption of Gases in Liquids, M.S. Thesis, Chem. Eng. M.I.T. (1942). See 
also O. Miyagi, Tech. Reps. Tohoku Imp. Univ. (Sendai) 8, 587-600 (1929), and references given there. We measured the 
rate of rise by a double flash exposure with a disk rotating at known speed. We estimated the bubble diameter by volumetric 
measurements, by silhouette pictures, and by measuring the distance between two reflecting spots on the bubble, assumed 
to be spherical. 
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be Nd?<30 cm—hence, if d=0.7 mm, bubbles should be able 
to support accelerations up to the order of 6000g, but not 
more. This may explain why we found it difficult to get good 
pictures of bubbles in the region immediately in front of 
spheres moving at 150 ft./sec. Of course, when bubbles are 
distorted, there is a small correction to be made to account 
for the shift in position of the spot on the bubble from which 
light is reflected. 

Silhouette pictures like Fig. 25 show the bubble distortion 
directly; in such cases, we measured the motion of an esti- 
mated centroid of the bubble. Fig. 25 is of oildrops under 
impulsive motion of a spheroid (Section 11); the bubbles near- 
est the spheroid were dished out at the center, and no longer 
even convex. Under similar circumstances, air bubbles simply 
exploded, appearing like bunches of grapes in silhouette photo- 
graphs. These observations place a further limitation on the 
accelerations to which the multiflash bubble technique can be 
applied. On the other hand, empirical evidence indicated that 
this limitation was less severe than the rough calculations 
outlined above had led us to suppose. 
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Fic. 25. Silhouette of 4X8-in. prolate hemispheroid with 
hexane bubbles. Presubmerged model struck by 5-lb. weight 
traveling 29 ft./sec. Time interval 0.0017 sec. Note distortion 
of hexane bubbles due to impact. 
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Dust—An Active Meteorological Factor in the Atmosphere of Northern Africa 


M. G. EL-FANby 
Asststant Professor, Farouk I University, Egypt 
(Received January 31, 1949) 


In subtropical Africa dust is shown to be an active meteoro- 
logical element, particularly below the 7-inversions common 
at 2000-3000 feet under Etesian conditions. Actual instability 
of the lower atmosphere is then not necessary for dust to rise, 
and mere mechanical stirring by day by winds of about 20 
miles per hour or more is sufficient to produce severe sand 
storms over desert areas under such conditions. 

A dust particle suspended in the atmosphere has been 
assumed to possess an “‘air-pocket” or “surrounding eddy” of 
its own. The radius of the “‘air-pocket” of a particle of radius 
“‘a"’ cm (less than 5 micron) is shown to be about 10a cm. 
It is doubtful whether one can assume the persistence of one 
and the same “‘eddy,”’ but with such an assumption made the 
average radius of the “eddy” is found to be 6a by observation, 
but on rejecting the assumption there appears to be a close 
agreement between theory and observation. The particles 
also share their excess of heat with their environment in a 
very small time, of the order of 10-?-10-* second. 

The amount of heat that can be directly supplied to the 
lower atmosphere by baked dust blown in a dust storm may, 
by itself, be sufficient to produce convection up to a height of 


3000 feet or more. The same range of height of convection is 
also attainable in a cloud of dust by the absorption of solar 
radiation scattered or diffusely reflected within the cloud. 
In this respect it is suggested, though not fully verified, that 
the effect of suspended dust should be taken into account 
when discussing problems of upper-layer instability in sub- 
tropical Africa. 

It has also been observed that in the special cases of 
“squalls” the intensity of the storm, notably the poorness of 
visibility, is proportional to the time during which the squall 
reaches its maximum. The feature has been given a mathe- 
matical explanation. ; 

The conclusion is that owing to increased absorption of 
solar radiation within dusty atmospheres, or to direct heating 
of the air by baked dust that may be blown by day in the 
warm season of the region, dust is found to play a part that 
cannot be ignored in modifying the processes of water vapor 
and heat transfer, by altering turbulent and radiative pro- 
cesses, and in modifying the properties of air masses invading 
the deserts. 





INTRODUCTION 


ARTICLES of dust suspended in the atmos- 

phere can have diameters larger than 0.1 
micron but not less. Samples of particles from 
layers of ‘‘dust-haze”’ or “drifting dust’’ at various 
heights, ranging from near the earth’s surface up to 
about 3 km, were collected in Petri dishes and were 
examined under the microscope. The maximum 
diameters of the particles so collected could be 
taken as 4 micron. 

Table | gives the limiting velocity of fall V of 
dust particles relative to still air, using Cunning- 
ham’s formula? and taking the density of dust to be 
2.64. Although in a turbulent atmosphere conditions 
may be different from those assumed relative to 
still air, yet it is clear from the table that the 
particles liable to remain in suspension for a long 
time over a wide area must have radii less than 
5 micron. 

During the warm season of the region (March to 


TABLE I. Limiting velocity of fall of dust particles 
relative to still air. 











Limiting velocity Limiting velocity 





Radius ‘‘a”’ of fall V of fall 
“ cm/sec. meters/day 
0.1 0.0007 0.6 
0.5 0.0095 a3 
1.0 0.0346 30 
2.5 0.2035 176 
5.0 0.797 690 











'G. C. Simpson, Quart. J. Roy. Met. Soc. 67, 115 (1941). 
?E. Cunningham, Proc. Roy. Soc. London 83, 357 (1910). 
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September) the establishment of a T-inversion or 


reduced lapse-rate at about 2000 to 3000 feet above | 


the earth’s surface, together with the increase of 
surface wind above a critical value of about 20 
miles per hour will cause the development of local 
sand storms which may extend vertically to 1000 
feet or more. 

An example of this type occurred over the 
Western Desert on 30th April, 1941 and visibility 
was reduced to 100 yards. Figure 1 gives the dry- 
bulb readings from Sidi Barrani on that day. 
Surface heating by day together with heating of 
the surface air by baked blown dust soon established 
a super-adiabatic lapse-rate to about 2000 feet, as 
shown in Fig. 1, and the intensity of the storm was 
thus increased. 

This type of sand storm is also much activated 
when local conditions are in favor of strong surface 
winds. An outstanding area is the Upper Nile 
Delta and Cairo District, where the ‘‘orographic 
bottle-neck effect’’ between the Mokattam Hills 
and the hills east and northeast of Faium (Fig. 2) 
favors great increase of wind speed near Cairo with 
S and SSW winds. 

It will be seen, however, that owing to absorption 


‘of solar radiation within clouds of dust, or to direct 


heating of the lower atmosphere by baked dust 
that may rise by day, the eddy diffusivity K is 
highly increased. Thus, representing K by the 
usual notation® 

K =}wd, (1) 


3D. Brunt, Physical and Dynamical Meteorology (Cambridge 
University Press, London, 1939), pp. 224-237. 
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where w is the term for the vertical velocity and d 
is the vertical extent since: 


(1) In severe dust storms, such as ‘‘haboobs’’ or severe 
“khamsin’’‘ dust may rise to about 3000 meters in an average 
of about 3 hours, K amounts to-about 10° cm?/sec. 

(2) In the common cases of air layers of ‘‘drifting dust’”’ or 
“dust-haze”’ upper air ascents show that such atmospheres 
usually have sharp tops and that the top of a dusty layer 
usually undergoes an elevation by day. This elevation amounts 
to about 500 m or more in 6 hours during the warm season of 
the region, and is due to increased absorption of solar radiation 
within the dusty layer. The increase of turbulence and conse- 
quent redistribution of dust in the vertical by day also 
produce an apparent elevation of the base of the dusty layer, 
of the type shown in Table II, which gives the hourly surface 
observations of horizontal visibility at Almaza (Cairo District) 
on 3rd July, 1941. Dust was apparently elevated from the 
surface at about 0900 G.M.T., but came down again by 
1630 G.M.T. 

The upper air report from the same district (Heliopolis) on 
the same day stated: ‘‘Cloudless sky. Haze-top between 0300 
and 0400 G.M.T. at 9500 feet. Horizontal visibility good at 
10,000 feet.’”” Between 0855 and 1025 G.M.T. another report 
stated: ‘Poor visibility at 10,000 feet. Haze-top at 11,000 
feet.’ In other words, the haze-top was elevated by 450 m in 
6 hours, so that K amounts to 5 X 104 cm?/sec. in this particular 
case. 

(3) With sand storms below 7-inversions or reduced lapse- 
rates dust rises to about 2000 feet in about 3 hours, so that 
K amounts to about 105. 


ROLE OF DUST PARTICLES IN TURBULENT 
TRANSFER OF HEAT AND MOISTURE 


A solid sphere moving through the atmosphere 
is assumed to be associated with a thin film of air. 
Pressure field is only disturbed within this film, 
and so long as the sphere is rising it must have an 
acceleration greater than the earth’s acceleration g. 
If ‘a’ is the radius of the solid sphere, the pressure 
“pb” at any distance “r’’ from the center of the 
sphere is given by Stokes’ approximation® 


p=constant+ 3uU(a cosé)/r’, (2) 
where —U=velocity of uniform stream of air of 
viscosity p, 

= —2/9(0—p)(ga*/u), 


in which o and p are respectively the density of 
sphere and air. 
Let z=rsin@ and x=rcos@, so that p=constant 
$uU(az/r*), and 


(0p/dz) = (3/2)nU(a/r*) — (9/2) U (az?/r>). 


The boundary layer is defined such that at its outer 
surface 0p/dz must be equal to —gp, so that 


(r/a)* = (sin?6—}3)(o¢—p)/p, (3) 
which represents two double solids (Fig. 3), the 





sae G. El Fandy, Quart. J. Roy. Met. Soc. 66, 323-336 


*L. Hoph, Handbuch der Physik (Verlag Julius Springer, 
Berlin, 1927), Bd. 7, 2, pp. 107-108. 
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— = EARLY MORNING ASCENT 
—<-— = MIDDAY ASCENT 


(SIDI BARRAN! 30 APRIL 194i) 


Fic. 1. Sidi-Barrani dry-bulb readings on tephigram 
(Western Desert 30/4/1941). 


total volume of which is Q and is equal to the 
volume of the boundary layer. 
Therefore, 


o=2f redx+2f rx'dz, 
0 0 


where s* = 4a3(o¢—p)/p. 

If ‘‘b’’ is the radius of a sphere whose volume is 
equal to Q, then b equals about 13a, and for all 
practical purposes 


b=10a, (4) 


is a fair approximation, taking o/p to be about 
2000 for dust in air. 

It will now be assumed- that every boundary 
layer is a physical representation of an ‘‘eddy.”’ 
In other words, it will be assumed that one and 
the same physical process is responsible for turbu- 
lent transfer of water vapor and dust particles. 
Usually layers of air of relatively high moisture and 
dust content exist below temperature reduced 
lapse-rates or inversions, such as those formed by 
subsidence within the northerly ‘‘Etesian’”’ winds 
of summer. 

Let n be the average number of dust particles 
per cc of air (=concentration) at a certain level 
which has become rich both in water vapor and 
dust content. The mass of water vapor transferred 














Fic. 2. “Orographic bottle-neck effect’’ between Mokattam 
Hills and the hills E and NE of Faium. 
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TaBLe II. The dusty atmosphere of Almaza, 3rd July, 1941. 











Time G.M.T. Visibility in code figure Present weather in code 
0345 6 = (4-10 km) 05 = (dust haze) 
0445 6= 05 
0545 6 05 
0645 6 05 
0745 6 05 
0845 6 05 
0945 8 = (20-50 km) 00 = (clear) 

1200 8 00 
1500 8 00 
1800 6 05 








mechanically to the layer by the vertical diffusion 
of dust will be 4/31b’p.nx, gram per unit volume, 
where p; is the air density and x, the average value 
of the humidity mixing ratio (in grams per gram of 
air) near the earth’s surface during the rise of dust. 
The net increase of water vapor within the level 
will be 4/3rnb*(x1p1—x2p1) gram per unit volume, 
where x2 is the original average value of x within 
that level (i.e., before the addition of dust). Let 
the increase per gram of air be denoted by Ax, 
so that 

Ax =4/3rnb*(x1—x2)p1/p2. (S) 


To compare this theory with observation, con- 
sider the dusty atmosphere of 3rd July, 1941, which 
prevailed over Egypt after sand storms of the above 
described reduced lapse-rate type over the Western 
Desert on July 2, 1941. Heliopolis’ early morning 
tephigrams of the 2nd and 3rd together with the 
mid-day dry-bulb readings of the 2nd are reproduced 
in Fig. 4. From these tephigrams the curves of 
Fig. 5 have been deduced to find the vertical 
distribution of x (per kilogram of air in this case) 
and the wet-bulb potential temperature 6’. Con- 
tinuous curves such as ABCD represent the early 
morning distribution of the 2nd, and dotted lines, 
such as A’B’C’D’, represent the corresponding 
distribution of the 3rd. These curves show clearly 
that although the same air mass prevails between 
2000 and 10,000 feet on both days, yet there is an 
average increase in x by 0.8 gram and an average 
increase in 6’ by about 3.5°C within 24 hours. 

The turbulent transfer of water vapor by air and 
by particles of dust takes place in proportion to 
their masses, so that on allowing for proper units 
and applying Eq. (5) to the above example we have: 


Ax =0.8(4/32a'n)o/p, 





z Z 





WH 
@D 
KY 

















Fic. 3. The two double solids of Eq. (3). 


662 


where o/p=3000 in this case, x,;=13.5, x2=4.5, 
pi/p2=1.08. Therefore, b?=240a* or b=6a approxi- 
mately. 

On the other hand, if instead of giving x; the 
value near the earth’s surface only it is given an 
average of 7, which in our example corresponds to 
the whole path in the turbulent atmosphere, we 
get for a fair approximation 


b=10a, (6) 


as before, which shows a close agreement between 
theory and observation. 

Again, when a dust particle at any temperature 
is blown into the air a rapid process of sharing 
heat between the particle and its environment takes 
place. The thermal capacity H, of a dust particle 
of radius ‘‘a’”’ and specific heat c together with its 
accompanying pocket of average radius “‘10a”’ is 


H,=4/32a*(oc+Cp), 


where Cp=sp. ht of air at constant pressure 
= 0.24 cal./g. 


1000p = 1 approximately, 
and 
c=0.188 cal. per gram. 


sé 9 


As an example let “a 
value 10-* cm, so that 


H,=7-10-" unit. 


be given the convenient 


Let the particle be originally in an environment 
of temperature Ty, and be suddenly projected into 
another environment of temperature 7,, where T; 
is greater than 7,. (This is almost always the case 
when dust is blown by day.) In a small interval of 
time dt, the amount of heat conducted away from 
the particle to its new environment will be H, 
where, for a first approximation 


H2.=4na?k(q/a)dt unit, 


in which is the air-conductivity (=k’pCp; k’ being 


40° So” 60° 70° 8° 90° 








-$cRAMS 
HELIOPOLIS 


o—o = 2/7, 
~~ = 3/4 \enr 
wee = 2/7/41 0850 amt 


Fic. 4. Early morning Heliopolis tephigrams of 2nd and 3rd 
July, 1941, and mid-day dry-bulb readings of 2nd. 
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the molecular conductivity), which can be taken as 
6X10-* c.g.s. unit, and g is the difference between 
the temperature of the particle and that of the 
outer atmosphere during the time interval dz. 

The average temperature difference between the 
particle and its new environment during the time ¢ 
required for thermal equilibrium is (Ta—T.)/2, so 
that in our example above the total heat that can 
be transferred by conduction is 


t 
f H.dt=2(T4—T,.)t- 10-8 
“9 


unit approximately. 

The amount of heat carried out from the bound- 
ary layer by the process of replacement during the 
time dt is 


H;=4nb?FpCpq'dt, 


where gq’ is the number of °C of the boundary layer 
above the new environment during the interval dt, 
and F is the rate at which air within the boundary 
layer may be replaced by normal air from the 
turbulent region outside (in cc per unit area). 
After time ¢ the particle and its pocket become 
in thermal equilibrium with their new environment, 
and thé total heat conducted by replacement will be 


t 
f H;dt=1/2(Ta—T,.) Ft-10~ 


0 


approximately, where q’ is given the same average 
value as q. 

Equating heat lost and heat gained during such 
processes we get 


t 
Hi(Tz—-T?') = f (Hy-+H;)dt (7) 
0 


where 7,’ is the final temperature attained by the 
dust particles and environment. 
In the example under consideration we have 


(Ta—T.’)-10-? =(Ta—T.)t- 10-8 
+3(Ta—Ta) Ft-10~, 


and on collecting terms in ¢ we get, for a first 
approximation, 


t=(Ta—T.’)[(Ta—Ta)(10+ -SF)}"'-10-* (8) 


which shows that for measurable temperatures in 
the atmosphere and for all values of F, ¢ is only a 
very small fraction of a second. 

Aitken® also made the observation that dust motes 
are not affected by solar heat focused by a large 
lens, while larger objects burn up at once under the 


* J. Aitken, Trans. Roy. Soc. (Edinburgh) 42, 489 (1902). 
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Fic. 5. Vertical distribution of x and 6’ from Heliopolis 
early morning ascents on 2nd and 3rd July, 1941. ABCD= 
curves for 2nd. A’B’C’D’=curves for 3rd. 


same conditions. The explanation was that the 
mote lost heat very rapidly to the environment. 


ROLE OF DUST PARTICLES IN DEVELOPMENT 
OF CONVECTION 


The relation between the concentration m at any 
level z and the limiting velocity of fall V to the 
degree of turbulence in the air is given by’ 


on @ On on 
—-_(x-) +V—. (8) 
Ot dz 02 02 

Upper air observations show that in dusty 
atmospheres practically dry-adiabatic lapse-rates 
are established, except for reduced lapse-rates or 
isothermal distributions within the first 1000 to 
2000 feet, so that K can be treated as a constant 
independent of height within dusty atmospheres 
and ‘ 


dn/dt = K(d°n/dz*) + V(dn/dz). (9) 


7C. W. Normand, Indian Met. Office (Calcutta) 22, VI 
(1921). 
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TaBLe III. Lower Egypt monthly actinometric 
averages in °C.* 








March April May June July Aug. Sept’ 
(a) 62.3 67.2 70.0 72.9 72.4 71.9 69.4 
(b) 41.3 46.2 49.6 52.1 52.0 $1.3 49.3 
(c) 31.0 36.0 39.6 42.2 41.9 41.4 39.6 











* (a) = Black-bulb in vacuum; (b) = Bright-bulb in vacuum; (c) = Maxi- 
mum air-temperature in standard screen. 


Under a steady state (when dn/dt=0), let the 
top of the “haze-layer’’ be at a height 4 above an 
arbitrary chosen base and put V/K=B. The 
boundary conditions are 


n=no at z=0 and nm=amy at z=h, where 0<aCl. 
Therefore, 


nol a—e~8*+(1—a)e~** | 


n= LS 


i-e™* 


(10) 


Integrating from z=0 to z=h in order to get the 
total number N of dust particles contained within 
a vertical cylinder of unit cross section we have 


h 
v= f ndz 


nol (a—e~8")h+(1 —a)(1 —_ BS) B| 


(11) 


When, for example, the range of horizontal 
visibility at the base is half that at the top, a=}. 
Usually K is large and in dust storms it approaches 
the order of 10° or 10’ indicating active turbulence 
and a=1, so that N=npoh. 

The thermal capacity C of N particles of dust is 
4/3xa*acCN, and on making the following substitu- 
tions for a moderate dust storm: h=3X10° cm, 
a=2 micron, and n=10*, we get C=5 units. 

It is evident from Table III that the maximum 
temperatures attainable by desert sand during the 
warm season of the region are on the average 25°C 
above those attainable by the surface air, so that 
the maximum amount of heat that may be supplied 
to the lower atmosphere by baked blown dust is 
125 cal. per cm’. 

Adopting the principle used by Gold® in calcu- 
lating the height of convection resulting from the 
addition of the above heat energy to the lower 
atmosphere, we see that the height of such con- 
vection amounts to 3300 feet. In other words, the 
energy so. supplied by dust is sufficient to create 
convection up to more than one-third of the 
assumed height of the storm. This may be taken 
as an explanation of the commonly observed 


* E. Gold, Professional Notes No. 63, 5-6 (1933). 
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feature that the warm season dust storms of the 
region have a marked tendency to attain their 
maxima in the early afternoon.® 

In Egypt accurate readings of desert temperature 
are not available, but regarding dust as inter- 
mediate between black and bright, such readings 
can be taken to lie between the values given in the 
rows marked (a) and (b) in Table III, in which 
monthly maximum actinometric averages in °C for 
lower Egypt are reproduced.'® The readings marked 
(a) are of black-bulb in vacuum, while (b) are of 
bright-bulb in vacuum. Maximum air temperatures 
in standard screen are given in (c). 

Returning now to the dusty atmosphere of 3rd 
July, 1941, in which the top was elevated by 1500 
feet through absorption of solar energy, we see that 
if such suspended dust is considered to be in 
vacuum (i.e., isolated from the air to prevent 
conduction of heat) its temperature will on the 
average be not less than 25°C higher than the 
temperature of the air in which it is suspended. In 
other words, C=1 unit in this case and corresponds 
to a concentration »=2X105 particles of radius 2 
micron per cc of air. 

Since upper layer instability may pass into the 
thunder-storm stage by convection during the day, 
the presence of sufficient dust in the atmosphere of 
subtropical Africa should introduce some modifica- 
tion to our present conception regarding this prob- 
lem of instability. However, it was suggested by 
Farquharson" that the brownish shade nearly 
always observed at the base of the high cumulus 
cloud at Khartoum is evidence that even on days of 
light winds upward drift of sand occurs by day. 


ROLE OF DUST IN RADIATION BALANCE 


As half of the total energy in the radiation from 
the sun can be taken as contained within the range 
0.4u to 0.7y, i.e., in the form of light, the average 
diameter of suspended particles is large compared 
with these radiations. Diffuse reflection takes place 
and Rayleigh’s argument” can be used. 




















Fic. 6. Wind velocity-record sine-curve. 


*M. G. El-Fandy, Quart. J. Roy. Met. Soc. 66, 323-336 
(1940). 

10 Ministry of Public Works, 1938, Climatological Normals 
for Egypt and the Sudan, Cairo, Egypt. 

u J. S. Farquharson, Quart. J. Roy. Met. Soc. 63, No. 271 
(1937). 

2 W. J. Humphreys, Physics of the Air (McGraw-Hill Book 
Company, Inc., New York, 1929), pp. 573-577. 
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On the other hand, the point of maximum 
intensity in the terrestrial radiation, with an aver- 
age temperature of 30°C, has a wave-length of 
about 10u. In other words, the greater part of the 
effective wave-length of the outgoing radiation is 
greater than 5u. Suspended dust, therefore, mainly 
scatter radiations from the earth. However, it is 
made clear by Humphreys that a veil of dust is 
more effective in shutting solar radiation out than 
it is in keeping terrestrial radiation in. 

Assuming the dust particles to be full radiators 
at their temperatures, consider a unit volume of 
unit cross-sectional area and let the maximum 
effective vertical radiation from the contained dust 
be 


2ra’n( ST"), 


where S is Stefans’ constant. Let T’’ be the effective 
temperature at which the particles are absorbing 
heat. The energy lost per unit volume will be 


2na°*n S(T" —T""), 


and the rate of cooling A7/At in °C per second 
will be given by 


(AT /At)pCp = 2xa*n S(T" —T"). (12) 


For all reasonable values of the terms in Eq. (12) 
this rate of cooling is negligible, and in northern 
Africa suspended dust is only a factor that increases 
the heat energy of the lower atmosphere. 

Bearing such factors in mind, consider a cloud of 
dust within which 80 percent of the solar radiation 
is ultimately absorbed and 10 percent of the out- 
going terrestrial radiation is also absorbed. The 
amount of solar radiation available in the lower 
atmosphere can be taken as 60 percent of the total 
amount received in the outer atmosphere as calcu- 
lated by Angot.'* Column 1 of Table IV gives the 
appropriate amount so received at latitude 30°N 
(in cm? on form 2810 tephigram, Air Ministry 
Meteorological Office, London). 80 percent of this 
amount is ultimately absorbed within the cloud 
and is represented in column 2. Columns 3 and 4 
give respectively the average effective temperatures 
at which the earth’s surface radiates and 10 percent 
of this energy in cm? on same tephigram. Column 
6 represents half the amount of heat radiated to 
space by the cloud of dust, assuming the particles 
to be 70 percent full radiators at the effective 
temperatures given in column 5. The other half, 
or even much more, may be taken to be absorbed. 
In column 7 is given the net energy gained by the 
cloud in cm? on the tephigram per day. 


MODIFICATION OF AIR MASSES BY DESERT 
TRAJECTORY 


An almost continuous process of adding dust and 
sand to the atmosphere of northern Africa takes 


'8 See reference 3, p. 112. 
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TABLE IV. Heating of air masses by clouds of dust. 








Month 18 





2b 3e 44 Se 6f 7s gb 
Jan. 13.3 96 290 3.7 280 115 2 #07 
Feb. 16.8 12.1 295 4.0 285 12.2 4 1.3 
March 20.3 14.6 300 43 290 130 6 2.0 
Apr. 23.5 16.9 305 46 295 140 7.5 2.5 
May 24.7 17.8 310 49 300 150 8 2.7 
June 26.7 19.2 315 52 300 150 9 3.0 
July 25.3 182 320 55 300 150 9 3.0 
August 23.5 16.9 315 5.2 300 150 7 2.3 
Sept. 21.3 154 310 49 300 150 5 1.7 
Oct. 17.3 125 300 43 295 140 3 1.0 
Nov. 14.2 10.2 295 40 290 130 1 03 
Dec. 123 89 290 3.7 285 12.2 0 00 








*60 percent of the solar radiation at lat. 30°N in cm? on Form 2810, 
(Air Ministry of London) tephigram. 

> 80 percent of column 1. 
+ ¢ Effective temperature of earth's radiation in °A. 

410 percent of radiation at temperatures given in column 3, represented 

cm? on tephigram. 

¢ Effective temperature of dust in °A. 

f } radiation from dust in cm? on tephigram. 

® Net energy gained by cloud of dust. 

b Rise of temperature in °C per day. 


place through dust storms or dust-raising winds 
and causes heating of the lower atmosphere. The 
net gain in energy given in column 7 of Table IV 
corresponds to an increase in temperature of the 
magnitude indicated in column 8. 

Heating is a maximum in the summer months 
and is almost ni] in December. The practical signifi- 
cance of such results lies in estimating how far 
air currents invading subtropical Africa and causing 
dusty atmospheres of any form become modified. 
The feature usually takes place on a large scale, 
and in the warm season air currents are highly 
modified by desert trajectory. For example, polar 
air that crosses Egypt takes about 4 to 5 days to 
reach the central Sudan and its temperature is then 
not less than 10°C higher than over Egypt. 


APPENDIX 


Some Characteristics of Surface Wind in 
Dust Storms 


The intensity of a dust storm, in particular the poorness of 
visibility, is observed to be governed not only by the increase 
of horizontal wind speed above a certain critical value, but 
also by the time interval that is seen to lapse in attaining a 
velocity maximum. The longer the wave on the wind velocity 
record the more intense the storm is. In a single marked 
squall that continues for an appreciable time, say, 15 minutes 
(which is the usual case with instability conditions), visibility 
falls rapidly while extensive rise of dust occurs. 

It is natural to infer to the increase in the vertical compo- 
nent of wind the variation in the intensity of the storm, or 
the variation in the density and vertical extent of suspended 
dust. Let r-seconds denote the period of a single wave on the 
wind record and assume a sine-curve for the velocity record 
as represented in Fig. 6. Also let the time be measured from 
O, so that time at O is zero. 

At any time ¢ let the horizontal wind velocity be represented 
by (u+u’), where w’ is the departure from the average value u. 
Denote the average maximum departure from u by (u’)w as 
represented in Fig. 6, and let (w, w’, (w’)w) represent the 
corresponding values for the vertical component. The increase 
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in the kinetic energy of a particle of mass m from time ¢ to 
time ¢+d is 


4md {| (u+u’)?+(w+w’)*}, 


and this is equivalent to a fraction 1/f of the energy taken in 
lifting the particle through a height (w+w’)dt, where f>1, 
since some of the energy for lifting is taken from the air, so 
that 

me(w+w’ dt = }mfd\|(u+u'?+(w+w’)*} 
or 


(w-+-w") = (f/2g)(d/dt) {(u-+u’)?+(w+w')?}. 


The average value of (w+w’) within the interval —7/4, 


+7/4 is such that 


1 T +r 
waa, (wtw')dt=a ed (utu'P+(w+w')} 
=4f/gr[u(u’)w+w(w’)av J. 
Therefore, 


(w’ wy = (gr /4f) — uu’ /w 


and to a first approximation the maximum departure from 
the usual vertical component is given by 


(w’)ay = g7/4f. (13) 
Equation (13) shows that the average maximum departure 


from the usual vertical component must increase with the 
period, and explains the observed feature. 





Some Factors Affecting the Rate of Precipitation Hardening in Cu-Be Alloys 


Paut A. BEcK 
Department of Metallurgy, University of Notre Dame, Notre Dame, Indiana 


(Received December 22, 1948) 


It was found that the lower peak hardness and faster overaging of a Co-containing Cu-Be alloy quenched 
from a low solution treating temperature is associated with a large amount of grain boundary precipitation. 
A similarly large amount of grain boundary precipitation, coupled with low peak hardness value, is obtained 
if the alloy contains a small amount of Cr impurity, even if the quenching temperature is normal. The grain 
boundary reaction resulting from the Cr impurity can be eliminated by using high solution treating tem- 


peratures. 


T has been known for a long time that the precipita- 
tion hardening process in Cu+2 percent Be is con- 
siderably affected by the addition of a few tenths of a 
percent of Co, Ni, or Fe. An outstanding effect of these 
additions is the considerable decrease in the rate of 
hardening and particularly in the rate of overaging at a 
certain temperature. This is shown in Fig. 1, where 
curve 1 represents the variation of hardness with aging 
time at 316°C for the binary Cu+2 percent Be alloy, 
and curve 2 gives similar data for an alloy with the 
same Be content, but containing in addition 0.38 per- 
cent Co. It is seen that, in comparison with the Co 
containing alloy, the binary alloy begins to harden at a 
considerably higher rate, reaches a lower peak hardness, 
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Fic. 1. Hardness vs. time at aging temperature (log scale) 


No. 2 hard strip. (1) Cu+2 percent Be alloy. (2) Cu+2 percent 
Be+0.38 percent Co alloy. 
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and then rapidly overages. Comparison of the micro- 
structures shows that in the binary alloy overaging is 
associated with the development and progress of the 
grain boundary reaction as described in a recent paper 
by Guy, Barrett, and Mehl.' As shown by these authors, 
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~ Fic. 2. Hardness vs. time at aging temperature (log scale) for 
Cu+2 percent Be+0.4 percent Co alloy after quenching from 
various temperatures (A) quenching temperature 760°C, (B) 
quenching temperature 787°C, (C) quenching temperature 
815°C, (D) quenching temperature 842°C. 


1 Guy, Barrett, and Mehl, Mechanism of Precipitation in Alloys 
of Beryllium in Copper, AIME RP. 2341 (Metals Technology, 
February 1948). 
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the area affected by the grain boundary reaction has 
lamellar structure and lower hardness than areas with 
general precipitation. During overaging this area 
gradually extends to cover practically the whole 
specimen. Under otherwise similar conditions the 
extent of the grain boundary reaction in the Co con- 
taining alloy is very much lower. The structure appears 
to be spheroidal rather than lamellar. Although even 
here the affected area increases with time to some ex- 
tent, normally it remains confined to the vicinity of the 
grain boundaries. This may explain the great stability 
of the Co containing alloy against overaging. The Co 
additions decrease the rate of both the general precipita- 
tion reaction and the grain boundary reaction. 

It was noted that by varying the quenching tempera- 
ture of a Co containing Cu-Be alloy, the rate and the 
mechanism of the precipitation reaction during subse- 
quent aging can be changed considerably. Figure 2 
gives hardness vs. aging time curves for four different 
quenching temperatures. Curve A represents harden- 
ing at 338°C after quenching from 760°C, which is con- 
sidered the normal quenching temperature for this 
alloy. Curves B, C, and D indicate the effect of increas- 
ing quenching temperatures. It is seen that the maxi- 
mum hardness is practically unaffected, but that the 
initial rate of hardening greatly increases with increas- 
ing quenching temperature. It is interesting to note 
that with increasing quenching temperature, the extent 
of the grain boundary reaction is greatly reduced. With 
quenching temperatures above about 787°C it practi- 
cally disappears. That the suppression of the grain 
boundary reaction here is not a result of an increase in 
the amount of Co entering into solid solution is indi- 
cated by the fact that the initial rate of hardening is 
at the same time increased. This is just the opposite 
of what may be expected from an increase in the Co 
content (see Fig. 1). If the quenching temperature is 
decreased below 760°C, the result is a great increase in 
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Fic. 3. Hardness vs. time at aging temperature (log scale) for 
Cu+2 percent Be+0.25 percent Co alloy quenched from (A) 
760°C, (B) 726°C. 
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Fic. 4. Cu+2 percent Be+0.25 percent Co alloy containing 
Cr impurity. Quenched from 760°C, aged for 4 hours at 316°C 
X 150. Large extent of grain boundary reaction. 


the amount of grain boundary reaction and, corre- 
spondingly, an increase in the rate of overaging. This is 
shown in Fig. 3, which indicates in addition that the 
initial rate of hardening, as well as the maximum 
hardness reached, is lowered by decreasing the quench- 
ing temperature below 760°C. Low quenching tempera- 
tures affect the rate of the grain boundary reaction and 
the extent of overaging in the same manner as a decrease 
in the Co content would, However, the initial rate of 
hardening is again affected in the opposite sense. Ap- 
parently, then, variations in the quenching temperature 
do not affect the rate and the extent of the grain bound- 
ary reaction by changing the effective amount of Co 
in solid solution. — 

For a given quenching temperature the rate of the 
grain boundary reaction in Co containing Cu-Be alloys 
was found to be greatly affected by the addition of Cr 





Fic. 5. Cu+2 percent Be+0.25 percent Co alloy containing 
Cr impurity. Quenched from 787°C, aged for 4 hours at 316°C 
X 150. Extent of grain boundary reaction decreased considerably. 
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Fic. 6. Cu+2 percent Be+0.25 percent Co alloy containing Cr 
impurity. Quenched from 815°C, aged for 4 hours at 316°C Xx 150. 
Extent of grain boundary reaction corresponds to that normally 
found in Cr-free alloy after quenching from 760°C and aging. 


in spectrographic quantities. In an alloy quenched from 
760°C, very small amounts of Cr will increase the rate 
of the grain boundary reaction to such an extent that 
the attainable maximum hardness no longer meets the 
standard specifications. However, it was found that in 
such Cr and Co containing alloys the grain boundary 
reaction may be progressively suppressed by using 
increasing quenching temperatures. This effect is similar 
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Fic. 7. Cu+2 percent Be+0.25 percent Co alloy containing Cr 
impurity. Quenched from 842°C, aged for 4 hours at 316°C X 150. 
Grain boundary reaction is almost completely suppressed. 


to that described above for the Cr-free alloy. But the 
small Cr addition increases the quenching temperature 
necessary for suppressing the grain boundary reaction 
by as much as 50°C. Figures 4, 5, 6, and 7 show the 
decrease in the amount of the grain boundary reaction 
with increasing quenching temperature in an alloy 
containing 0.3 percent Co and a spectrographic quan- 
tity of Cr. 
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Semi-Conducting Properties in Oxide Cathodes* 


N. B. Hannay, D. MacNarr, anp Appison H. Waite 
Bell Telephone Laboratories, New York, New York 


(Received September 16, 1948) 


It has been widely assumed, without adequate experi- 
mental verification, that barium-strontium oxide, as used in 
the oxide cathode, is an excess electronic semi-conductor. 
Accordingly, the electrical conductivity of (Ba, Sr)O has been 
studied as a function of temperature before and after activa- 
tion with methane, extensive precautions being taken to 
exclude spurious effects. The increase in conductivity obtained 
characterizes (Ba, Sr)O as a “reduction’’ semi-conductor, and 
hence very probably as an electronic semi-conductor whose 
conduction electrons arise from a stoichiometric excess of 
(Ba, Sr) atoms in solid solution. 

A basic prediction of the semi-conductor theory has been 


tested quantitatively with the finding that the electrical 
conductivity and the thermionic emission of a (Ba, Sr)O 
cathode are directly proportional through three orders of 
magnitude of activation; well-defined chemical and electrical 
activation and deactivation procedures were used in obtaining 
this result. It may- be concluded that activation represents an 
increase in the chemical potential of the electrons in the oxide, 
little or no change in the state of the surface occurring. It has 
also been found that deviations from the proportionality of 
conductivity and emission may be expected under conditions 
leading to inhomogeneity in the oxide, in agreement with the 
‘semi-conductor theory also. 





I. INTRODUCTION 


HE massive literature of the oxide-coated cath- 

ode shows it to be a very complex system 
which is by no means completely understood.' It 
is clear that the most nearly satisfactory theory of 
the emission of this cathode is based on the assump- 
tion that barium-strontium oxide (hereafter desig- 
nated as (Ba, Sr)O) is an electronic semi-conductor 
activated by an excess of, e.g., barium dissolved 
in it. Experiment is generally consistent with this 
semi-conductor picture, under conditions where the 
thermal equilibrium of electrons assumed for the 
theory is likely to be realized.” 

The semi-conductor theory of thermionic emission 
predicts a correlation between the thermionic emis- 
sion and the electrical conductivity if the oxide is 
an homogeneous excess electronic semi-conductor. 
If there is no change in the state of the surface 
during activation, the conductivity and emission 
should be directly proportional. A primary object 
of the present investigation has been to test this 
basic prediction by studying the emission and con- 
ductivity as a function of the state of activation 
of the ‘oxide, and thus obtain more information 
about the activation process than has hitherto 
been available. The theory underlying this work 
will be discussed briefly (Section II). 

The basic assumption that (Ba, Sr)O is an excess 
electronic semi-conductor has not, however, been 
adequately tested by experiments specifically de- 
signed for the purpose. Accordingly, the initial 


*A preliminary report of this work was made to the 
American Physical Society, Division of Electron and Ion 
Optics; Phys. Rev. 72, 153 (1947). 

‘J. P. Blewett, J. App. Phys. 10, 668, 831 (1939) ; ibid., 17, 
643 (1946). 

2 J. A. Becker and R. W. Sears, Phys. Rev. 38, 2193 (1931); 
W. Albricht, Physica 11, 146 (1931); A. L. Reimann and L. R. 
G. Treloar, Phil. Mag. 12, 1073 (1931); W. Meyer and A. 
Schmidt, Zeits. f. tech. Physik 13, 137 (1932) ; E. Nishibori and 
H. Kawamura, Proc. Phys. Math. Soc. Japan 22, 378 (1940) ; 
A. Fineman and A. Eisenstein, J. App. Phys. 17, 663 (1946). 
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phase of the investigation has been to test this 
assumption by studying the effect of high tempera- 
ture treatment with a chemical reducing agent, 
methane, on the electrical conductivity of the 
oxide. A marked increase of conductivity under 
such treatment will characterize (Ba, Sr)O as a 
‘“‘reduction”’ semi-conductor, and thus very prob- 
ably as an excess electronic semi-conductor.* 4 


II. SEMI-CONDUCTOR MODEL FOR OXIDE CATHODE 


The theory is available in the literature,’ 5* and 
so only the part pertinent to this work will be 
briefly reviewed.’ 

Without making any explicit assumptions about 
the exact model to be used, it may be shown easily 
that the conductivity and emission should be 
related. The conductivity will be approximately 
proportional to the concentration of ‘‘free’’ elec- 
trons in the conduction band of the solid, m;, at 
any fixed temperature. 


o=kin,. (1) 


The equilibrium thermionic emission current will 
be proportional to the concentration of conduction 
electrons and to e~‘*/*?), y being the work function 
for these electrons; the conductivity-emission rela- 
tionship then follows immediately. 


i = kon e~(x/#T), (2) 
=k3ae~ x!/*T), (3) 


3N.F. Mott and R. W. Gurney, Electronic Processes in Ionic 
Crystals (Clarendon Press, Oxford, 1940), Chapter 5. 

4 Methane has been reported to activate the electron emis- 
sion of (Ba, Sr)O which contains finely divided nickel. (See 
C. H. Prescott and J. Morrison, J. Am. Chem. Soc. 60, 3047, 
1938). These authors do not discuss the role of conductivity. 

5R. H. Fowler, Statistical Mechanics (University Press, 
Cambridge, 1936), second edition, pp. 401, 419; see reference 1, 

1939). 
* Frederick Seitz, Modern Theory of Solids (McGraw-Hill 
Book Company, Inc., 1940). 

7 The authors are indebted to Dr. C. Herring for many 

helpful discussions of the theory. 
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Fic. 1. Energy levels in (Ba, Sr)O. 


The direct proportionality between the conduc- 
tivity and emission expressed by Eq. (3) can be 
derived also in a more explicit form by the usual 
type of mathematical treatment, while maintaining 
the generality of the argument. In Fig. 1 we let 
the energy of an electron just outside the semi- 
conductor be —e®%,, and at the bottom of the 
conduction band ¢€o, the zero of energy being taken 
at the average electrostatic potential inside the 
crystal, and #, being the electrostatic potential in 
the vacuum just outside the semi-conductor. The 
quantity x is the work function for conduction 
electrons used above, or the “electron affinity’’ of 
the crystal, and the quantity yw is the chemical 
potential of the electrons® (the partial molar-free 
energy of the electrons), or the ‘‘Fermi level.” 
Impurity levels are shown, which may be of the 
F-center type, but their nature need not be specified 
for this discussion. The conductivity is then given 
by the expression 


o = (8e7lo/3h*)(2erm*kT) exp[(u—eo)/RkT ], (4) 


where /)= the electronic mean free path, m* =effec- 
tive mass of the electron, and the other symbols 
have their usual meanings. The thermionic emission 
current under equilibrium conditions is given by an 
expression which may be derived either from ther- 
modynamics or from statistical mechanics, and 
which may be written as 


j=((2am*kT)2kTe/h* \(1 —F) 
Xexp[(u—eo—x)/kT], (5) 
where f=a reflection coefficient, and the true 


8 Actually, the electrochemical potential, @, should be used 
throughout, where f=y—e®,, &, being the average electro- 
static potential inside the crystal. The energy zero has been 
chosen so as to make fi =u for the sake of simplicity. 
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thermionic work function is given by 


g= —[u—eo—x]. (6) 


The chemical potential, of course, will have 
definite meaning regardless of the model chosen. 
For a particular model it may be written explicitly 
in terms of the impurity centers, but this is not 
necessary for purposes of the present discussion. 
It is shown in Fig. 1 in the vicinity of the uppermost 
impurity levels, but could be drawn elsewhere with- 
out affecting the validity of the argument. 

Upon activation of the semi-conductor, the chemi- 
cal potential will increase, for example, by an in- 
crease in the number of impurity levels, and the 
increase in the chemical pdtential will increase both 
the conductivity and’ the thermionic emission, as 
shown by Eqs. (4) and (5). During this activation, 
if no quantities other than uw on the right hand 
sides of (4) and (5) change, the conductivity and 
emission will be directly proportional. If, on the 
other hand, x for example should change because 
of a change in the surface dipole, then this linear 
relationship would no longer hold. It should be 
emphasized that the validity of these predictions 
does not depend on a particular model for the semi- 
conductor. 

Results of Heinze and Wagener in electron micro- 
scope studies tend to indicate that the chemical 
potential does change during the activation process, 
with little change in the surface dipole.® !° Nishibori 
and Kawamura, Albricht, and others? have meas- 
ured the conductivity and emission of an oxide 
cathode, with the work of the first of these as 
perhaps the best experimental check, a linear rela- 
tionship being derivable from their three experi- 
mental points. The general inadequacy of the data 
made it appear desirable to investigate the con- 
ductivity-emission relationship much more exhaus- 
tively, and under activation procedures in which 
chemical and electrical processes were more clearly 
defined. 


Ill. CONDUCTIVITY MEASUREMENTS 


A condensed description of the experiments per- 
formed and the results obtained will first be pre- 
sented, and will be followed in Sections V—VI with 
a detailed discussion of experimental techniques 
and of possible spurious effects. The proper choice 
and execution of precautions designed to exclude 
such effects constitutes the major problem in work 
of this kind with impurity semi-conductors. 

Measurements of the conductivity of the active 
oxide in vacuum, as reported in the literature, have 
usually been incidental to studies of electron 


® W. Heinze and S. Wagener, Zeits. f. Physik 110, 164 (1938). 
mt. 4 Heinze and S. Wagener, Zeits. f. tech. Physik 20, 16 
1939). 
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emission. This has generally required that conduc- 
tion electrodes be separated by small! and uncertain 
amounts. The possibilities of error caused by elec- 
trode contact resistance or polarization, or by 
thermionic emission through the porous oxide have 
not all been considered. The specimens reported in 
the literature have not been prepared for measure- 
ment by well-defined chemical treatment; e.g., the 
usual procedure has been to activate the (Ba, Sr)O 
for electron emission by the conventional but 
chemically ambiguous method of drawing current 
through the specimen. 

The geometrical arrangement of the test speci- 
men used in the present work and its leads is shown 
in Fig. 2. A chemically inert insulating support for 
the oxide specimen was provided by a magnesia 
ceramic, heated internally by a helical alumina- 
coated tungsten heater, as shown. The electrode 
system consisted of four platinum leads, 0.003 
inch in diameter, mounted in appropriate grooves 
cut in the ceramic. This assembly was then coated 
by spraying with a standard cathode mixture con- 
taining BaCO; and SrCQ; in the ratio of 45 to 55 
by weight, suspended in a 1 percent solution of 
nitrocellulose binder in amyl acetate. The mean 
thickness of the porous carbonate specimen after 
evaporation of the amyl acetate was about 0.0035 
inch. This assembly was sealed to a conventional 
high vacuum mercury pump station which main- 
tained a pressure of <110-*-mm Hg throughout 
the subsequent experiments, except when their 
execution required treatment with purified gas. 

The carbonates deposited on the magnesia were 
converted in vacuum to the oxides by raising the 
pyrometer temperature of the specimen to 1270°K 
over the course of about 30 minutes. The tempera- 
ture was held at this level until the pressure dropped 
to below 10-*-mm Hg, which usually required about 
5 minutes. During this decomposition the maximum 
true temperature was from 1400 to 1475°K. This 
is higher than any temperature reached during the 
subsequent measuring or processing of the specimen. 

The (Ba, Sr)O was then heat-treated at about 
1370°K until the conductivity no longer changed 
appreciably with time. Effective conductivity was 
next measured as a function of temperature, with 
results as shown in the lower curve of Fig. 3. The 
data for one tube only are shown in the inactive 
state, for the sake of simplicity in the diagram; 
the curves in the inactive states for the other tubes, 
whose activated states only are shown, are of the 
same general nature. These data in the inactive 
state reflect partial activation of the conductivity 
due to uncontrolled factors such as traces of carbon 
introduced as nitrocellulose, and by comparison 
with each other are obviously not characteristic of 
a well-defined or reproducible state of the (Ba, Sr)O. 
The measurements were made to establish a basis 
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Fic. 2. Experimental tube for the measurement of conductivity. 


for comparison with the results after methane 
treatment. 

Unless otherwise specified, all temperatures re- 
ported are pyrometer readings corrected in accord- 
ance with the spectral emissivity of the specimen 
as described below in Section V. The primary 
measurement of effective conductivity was by a 
null method, also described below. The term “‘effec- 
tive conductivity” is employed to emphasize the 
fact that electrical measurements were converted 
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Fic. 3. Conductivity as a function of temperature in (Ba, Sr)O. 
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to conductivities without allowing for the porosity 
of the specimen. 

The activation of the (Ba, Sr)O as a semi-con- 
ductor consisted of heat treatment in an atmosphere 
of methane. This gas was admitted to the system 
at room temperature until a pressure of 1.0—2.5 
<10--mm Hg was reached. The specimen was 
then quickly heated to approximately 1320°K and 
held at this temperature for between 5 and 10 
minutes. It was next cooled to room temperature, 
and the system was again exhausted. The methane 
treatment generally increased the spectral emis- 
sivity of the specimen, for example from 0.15 to 
0.25, possibly because of deposition of invisible 
amounts of carbon by decomposition of the meth- 
ane. The oxide was then heat-treated in vacuum at 
the highest temperature of subsequent measurement 
until spectral emissivity declined to a steady value, 
which required about one hour. After this final 
treatment the (Ba, Sr)O was ready for measure- 
ment in its activated state. : 

The effective conductivity of the porous (Ba, Sr)O 
after treatment with methane is shown in the upper 
curves of Fig. 3. Clearly this conductivity is every- 
where larger after treatment, by a factor ranging 
as high as 3000. Furthermore, the increase co- 
incided with the initiation of the methane treat- 
ment. Crude resistance measurements made during 
activation showed that the entire change occurred 
during the first two minutes, although the heat 
treatment in the gas continued for 5 or 10 minutes. 
Thus it appears that the conductivity of (Ba, Sr)O 
is, in general, increased by treatment with methane, 
as would be expected if it is an excess electronic 
semi-conductor. 
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The rapidity with which steady-state values of 
the conductivity were obtained on initial activation 
is taken to be evidence that, at the temperature of 
the reaction, the concentration of impurity centers 
quickly reached a uniform value throughout the 
solid phase. This was probably due to the porosity 
of the specimen, no part of which was likely to be 
more than two or three microns from some pore 
surface. Repeated methane treatment did not in- 
crease conductivity appreciably over that observed 
after the initial activation. After pumping off the 
methane, heating at 1350°K for one hour did not 
change conductivity by more than 5 percent. The 
stability of this conductivity over the four or five 
hours required to complete a temperature run was 
repeatedly demonstrated by a final repetition of 
the initial measurement. 

The reproducibility of the results obtained from 
different activated specimens as illustrated in 
Fig. 3, shows that at no temperature do the meas- 
urements differ from their geometric mean by more 
than +50 percent. This is rather remarkable agree- 
ment for semi-conductors, particularly where the 
geometry of the specimens is so ill-defined, and 
suggests that a reasonably reproducible state of 
the system was attained. 

The most nearly comparable conductivity values 
in the literature are those of Meyer and Schmidt? 
for BaO, and those of Nishibori and Kawamura,’ 
apparently also for BaO. Although they were 
obtained from BaO which had been activated by 
the flow of current through the specimen (as well 
as by exposure to barium vapor by Nishibori and 
Kawamura), there is a striking general resemblance 
of our CH,-activation results with (Ba, Sr)O to 
those Meyer and Schmidt at high temperatures and 
to those of Nishibori and Kawamura at low tem- 
peratures. In both of these experiments measure- 
ments were made with nickel electrodes a fraction 
of a millimeter apart. 

The data recorded in Fig. 3 are in each case 
rather well represented by two straight lines with 
the lower slope at the higher temperatures. Similar 
behavior of conductivity as a function of tempera- 
ture has been observed in ZnO," and in incom- 
pletely activated BaO (Meyer and Schmidt).? 

The slopes characteristic of the respective linear 
portions of the conductivity curves in Fig. 3 are 
1.1 ev in the lower temperature range, and 0.3—0.45 
in the higher temperature range, these slopes being 
expressed as values (in electron volts) of the con- 
stant b in the equation 


o = Ae lkT), (7) 


1 W. Jander and W. Stamm, Zeits. f. anorg. allgem. Chemie 
199, 165 (1931); H. H. von Baumbach and C. Wagner, Zeits. f. 
physik. Chemie (B) 22, 199 (1933). 
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Conclusions 


These experiments justify the conclusion that 
(Ba, Sr)O is a “‘reduction”’ semi-conductor, i.e., it 
is activated as a semi-conductor by chemical treat- 
ment which tends to reduce the oxide to the free 
metal. The simplest hypothesis consistent with this 
result is that activated (Ba, Sr)O is an electronic 
semi-conductor because of a stoichiometric excess 
of (Ba, Sr) atoms in the crystal lattice. 


IV. CONDUCTIVITY-EMISSION MEASUREMENTS 


The experimental diode used in this work is 
shown in Fig. 4. The support was the same MgO 
ceramic used in the conductivity work, with the 
cathode core consisting of two parallel platinum 
wires, 0.003 inch in diameter, wound on the ceramic 
in the form of a double helix imbedded in the oxide. 
T he spacing of these in the central region was 0.01 
inch, center to center, this region 0.10 inch in length 
was that covered by the (Ba, Sr)O mixture. The 
wires were thus separated at all points by 0.007 
inch of the oxide. The anode was a bright nickel 
ribbon 0.2 inch wide surrounding the coated area 
of the cathode. The platinum core wire and the 
en used were the same as in the conductivity 
work, 


The tubes were left on the vacuum system 
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throughout the course of measurement. Starting 
with an inactive tester, the coating freshly broken 
down from the carbonates, activation was carried 
out in slow stages. Three methods were employed : 
heating alone, aging with space-charge limited d.c. 
current, and treatment with methane at 2 X 10-*-mm 
Hg pressure. The purpose of this was to compare 
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Fic. 6. Conductivity and emission showing deviations from 
linearity due to inhomogeneity in the oxide. 
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the more clearly defined methane activation with 
thermal activation, which presumably involves 
chemical reduction by traces of reducing agents 
present in contact with coating, and with the 
standard procedure usually used in the activation 
of oxide cathodes, which involves the drawing of 
d.c. emission current. In each case, the conductivity 
and emission were measured after each change in 
state of activation, the measurements being carried 
out in vacuum, with the coating temperature 700°C, 
corrected (630°C apparent pyrometer tempera- 
ture),* a temperature that was found to give no 
activation during the course of the measurements. 
Thus, all the points on the conductivity-emission 
curves, as in Fig. 5 for example, were taken under 
the same conditions, with successive points on a 
single line being obtained by the activation steps 
just described. When the maximum measurable 
activity was reached, the cathode was deactivated 
with oxygen, at 2X10-*-mm Hg pressure for 1-2 
minutes, at 700-800°C. The activation procedure 


TABLE I. Quantitative analysis of carbonates used for 
oxide specimens. 











BaCO; SrCO; 
Fe 0.0025% by wt. 0.001% by wt. 
Cl 0.006 <0.001 
Ca 0.05 0.20 
SO, <0.01 0.0024 
Na <0.001 0.02 
S (sulfide) 0.001 <0.001 
Acetic acid insol. 0.013 0.007 








*See Section V for comments on the temperature meas- 
urement. 
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could then be carried out as before. This activation- 
deactivation cycle was completely reversible, as the 
points show, no difficulty being experienced with ob- 
taining reproducible activation curves after several 
oxygen treatments. The initial activation was 
similar to subsequent activations except that the 
initial activation usually proceeded somewhat more 
slowly. The final level of activity was about the 
same in each case. In the activation by drawing 
d.c. current, a current was chosen that amounted 
to about one-half that at which departure from 
the space-charge line set in, and a temperature was 
chosen at which activation by heating alone was 
considerably slower. 

Figure 5 summarizes the results obtained in the 
comparison of conductivity and thermionic emission 
over a wide range of states of activation for one 
tube. Five tubes were tested, the data for all being 
entirely similar to those in Fig. 5. Each set of points 
corresponds to a single activation cycle, several 
such sets appearing in the figure. In no case were 
fewer than six points obtained for an activation, 
and as many as twenty were sometimes obtained. 
In all cases the points deviate from the 45° reference 
line by no more than the experimental error in the 
emission measurements, which is the principal 
cause of the scatter in the points. Two heat activa- 
tions, one d.c. aging activation, and two CH, 
activations are shown. In addition, a deactivation 
with O-2 is shown. This was obtained by starting 
with a completely activated tube, and proceeding 
with progressive deactivations by treatment with 
Oz, at 860°C, for about 2-10 seconds, at 10-*-mm- 
10-*-mm Hg pressure. It will be observed that 
excellent agreement with the other points was 
obtained. Other data for this tube are shown in 
Figs. 6 and 7. Entirely similar data were obtained 
on the other tubes. 

Comparison of the data for the various tubes 
shows that not only are the activation lines for 
each tube self-consistent with a slope experimentally 
indistinguishable from 45° in all cases, but also 
that the various sets of lines are nearly super- 
imposable, the difference between the sets for 
different tubes being no more than the scatter 
among the individual lines for one tube. This 
indicates a remarkable degree of uniformity from 
tube to tube in such matters as conditions of the 
oxide surface, as any changes would shift the lines 
along one axis. 


Deviations from Linearity 


The discussion of Section II implies that the 
oxide semi-conductor be homogeneous throughout 
activation, if a linear conductivity-emission rela- 
tionship is to be expected. Chemical treatments, 
then, must affect the interior of the oxide crystals 
in the same way as the surface layers for this to 
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hold. Herring’? has discussed the theoretical aspects 
of thermionic emission, concluding that usually the 
emission is limited by the thermodynamic state of 
the outer few thousand angstroms of the oxide, at 
most. If activation were to change only these 
surface layers, or the electron affinity due to a 
change in the surface dipole, a linear conductivity- 
emission relationship would not result. If, on the 
other hand, activation changes only the chemical 
potential » in the whole body of the oxide, leaving 
the electron affinity x unchanged, the linear rela- 
tionship would hold. 

Under certain conditions, it was found possible 
to obtain results that showed deviations from 
linearity, where this sort of inhomogeneity in the 


body of the oxide would be expected. These experi- ° 


ments are shown in Fig. 6. 

The first experiment involved the deposition of 
carbon on the oxide particle surfaces. The data in 
Fig. 6 show the initial state of activity near the 45° 
line and the final state after deposition of carbon 
by treatment with excess CH, at a relatively high 
temperature. This is a process which has changed 
the outer layers but not the interior of the oxide, 
directly contrary to the homogeneous changes in 
the experiments of the last section, and results in a 
decrease of emission without a corresponding de- 
crease in conductivity. The conduction paths for 
the electrons through the solid were essentially 
unaltered, whereas the surfaces from which therm- 
ionic emission occurred were completely altered by 
the visible deposit of carbon. 

Figure 5 showed an oxygen deactivation curve 
which followed the same 45° straight line obtained 
in the activation curves. This was obtained by 
deactivating with oxygen at 860°C. In Fig. 6 are 
shown data representing two cases of deactivation 
with oxygen at room temperature, and one de- 
activation at 640°C. The treatment was for 1-2 
minutes, at an Oz» pressure of 2X10-?-mm Hg. In 
each case of room temperature deactivation, the 
emission dropped sharply upon O, treatment from 
the initial point shown near the 45° line without a 
corresponding decrease in conductivity. As the rates 
of diffusion at room temperature are undoubtedly 
very low, it is probable that the oxygen deactivated 
only the surface layers. As the temperature was 
raised to 700°C for measurement, diffusion would 
set in, of course, but not too rapidly for a measure- 
ment of the conductivity and emission before 
appreciable changes had occurred. The lowering of 
the chemical potential near the surface would 
serve only to decrease the thermionic emission, 
without greatly changing the bulk conductivity. 
At 640°C, the rates of diffusion were still not high 





® C. Herring, private communication. 
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TABLE II. Mass spectrometer analyses of purified gases. 











CH, 97% O. 99% 
Other hydrocarbons 0.7 A 0.8 
Ne 0.66 

H. 1.38 

Oz none detected (<0.005%) 











enough to give a completely homogeneous coating, 
apparently. Por 305 3 

In the. final experiment shown in Fig. 6, a d.c. 
current was drawn from a fully activated tester 
with the plate voltage considerably in excess of 
that necessary to draw saturation current. Under 
these conditions a rapid deactivation, with regard 
to thermionic emission, set in, in line with general 
experience for these conditions, and as Fig. 6 
shows, without a corresponding decrease of con- 
ductivity. It cannot be stated definitely what 
caused the deactivation of the surface. 


Conductivity—Temperature Relationship 


The tube used for the measurements in Figs. 5 
and 6 was examined in its maximum state of 
activity over a considerable temperature range, to 
see whether the log conductivity vs. 1/T curve was 
similar to those obtained in the conductivity work, 
the data being shown in Fig. 7 along with a typical 
curve obtained in the conductivity work. In view 
of the uncertainties in the geometrical considera- 
tions in calculating conductivities from conduct- 
ances, and in the relative degree of activation, it is 
certainly demonstrated that the two curves are 
measurements of essentially the same phenomena. 
Unfortunately, the break in this curve occurs as a 
temperature at which the emission measurements 
fail because of sparking. The curve cannot be 
obtained in a lower state of activation, where the 
emission measurements would be valid, as the tube 
would activate quite rapidly at temperatures ap- 
proaching that of the break in the conductivity 
curve. Thus the possibility of a similar break in 
the emission curve could not be investigated. 


Conclusions of Conductivity—Emission Studies 
It has been found that the electrical conductivity 
and the thermionic emission of a (Ba, Sr)O cathode 


TABLE III. Spectrochemical analysis for contamination 
of (Ba, Sr)O. 








Major Minor Impurity Trace 





(>1%) (0.1-3%) (0.01-0.3%) (<0.03%) 
Magnesia ceramic Mg Si Al, Ba, Ca, Cr, Ag, Cu, Mo, 
as received Fe, Mn, Sr Na, Ni, Ti 
Carbonate spray Ba, Sr Ca Na Al, Bi, Cu 
mixture Fe, Mg, Mo, 
Pb 
(Ba, Sr)O after Ba, Sr Ca Fe, Mg, Na Al, Cu, Mn, 
conductivity Ni, Si 


tests 
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with a Pt core are directly proportional through 
three orders of magnitude of activation: the activa- 
tion was accomplished by heating without drawing 
current, by heating while drawing space-charge 
limited d.c. current, and by reduction with methane 
gas, and a deactivation curve was obtained with 
oxygen at 860°C. In all cases the same linear result 
was obtained, regardless of the method of activation 
or deactivation; furthermore, on several experi- 
mental tubes exactly the same results were found, 
no appreciable shift of the lines from tube to tube 
being observed. Thus a quantitative relationship 
involving the thermionic emission, predicted by 
the semi-conductor theory for the activation process, 
has been experimentally verified within the limits 
of error of measurement. For a cathode of this 
particular geometry, if one is assured that equi- 
librium conditions prevail, a knowledge of the 
electrical conductivity enables a prediction of the 
thremionic emission current to be made. It may 
be concluded that activation changes only the 
chemical potential, little or no change in the electron 
affinity, x, occurring. This experimental result is 
inconsistent with theories of activation involving 
changes in free barium on the surface of the oxide, 
such as have been proposed in the past. Moreover, 
this proportionality of conductivity and emission 
would not be expected if there were an appreciable 
number of surface states on the oxide. 

In addition, it has been found that certain 
processes which would be expected to change only 
the surface layers of the oxide will affect the 
emission without a corresponding change in con- 
ductivity, these processes being (1) deposition of 
carbon on the crystal surfaces with excess CH, 
treatment; (2) oxygen treatment at room tempera- 
ture, with inhomogeneity arising from the low diffu- 
sion rates at this temperature; and (3) processes 
arising from the drawing of saturation d.c. current. 
These observations are consistent with the view 
that thermionic emission is determined by the 
thermodynamic state of the outer atomic layers. 


V. DETAILS OF EXPERIMENTAL TECHNIQUES 


Some of the precautions taken to insure the 
validity of these results will be described in this 
and the next section. 


TABLE IV. Resume of spectral emissivity measurements on 
(Ba, Sr)O specimens. 








Number of 
: measure- Temperature 

Specimen ments range Ey (A =0.66n) 
C-2 before activation 11 300°-1350°K 0.17 +0.015 
after activation 8 300°-1270° 0.21 +0.025 
C-1 before activation 6 300°-1340° 0.10 +0.015 
after activation 9 300°—1320° 0.09 +0.025 
C-3 before activation 7 990°—1350° 0.14+0.025 

after activation 4 980°-1140° 


0.15 +0.015 
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Chemistry of Materials and Processes 


The principal impurities found in the original 
carbonates are shown in Table I. 

The gases used in the course of investigation, 
methane and oxygen, were purified before admission 
to the vacuum system, their purity being checked 
by mass spectrometer analyses. The methane was 
passed over hot copper and anhydrous calcium 
sulfate, and also over Ascarite in the case of the 
conductivity-emission work. The oxygen was passed 
over Ascarite and anhydrous calcium sulfate. In 
each case the liquid-nitrogen trap remained on the 
system. Typical analyses after purification are 
shown in Table II. 

The methane used in -the conductivity experi- 
ments also contained 0.25 percent COs, but this 
was removed by condensation in a liquid-nitrogen 
trap before admission to the (Ba, Sr)O specimen. 

One of the chief problems was to select for the 
insulating support of Figs. 2 and 4 a material 
which would not react while in contact with 
(Ba, Sr)O at high temperatures. Magnesium oxide 
fired in an oxidizing atmosphere was chosen because 
of its failure to form any known compounds with 
alkaline earth oxides.'* X-ray and electron diffrac- 
tion data obtained by Mr. K. H. Storks of these 
laboratories showed no evidence of reaction between 
BaO and MgO after heating the powdered mixture 
to 1000°C in hydrogen. At the conclusion of one of 
the conductivity experiments, some of the (Ba, Sr)O 
was scraped off its support and compared by quali- 
tative spectrochemical analysis with the original 
carbonate suspension. The results are recorded in 
Table III. It is apparent that the ceramic did not 
react appreciably with the (Ba, Sr)O and indeed 
contaminated it only slightly. The only new im- 
purities (0.01-0.3 percent) are iron and magnesium, 
presumably as oxidés from the ceramic. In these 
quantities they cannot affect the result of methane 
treatment except for the remote possibility that 
(Ba, Sr)O is a semi-conductor by virtue of unoxi- 
dized iron or magnesium dissolved in it. 

It was also essential to determine the effect of 
heat treatment in methane on the insulation re- 
sistance of the magnesia support. After successive 
treatments in the absence of (Ba, Sr)O at 1400°K 
in vacuum and 1300°K in methane, this resistance 
between leads 2 and 3 of Fig. 2 remained too high 
(> 10° ohmsat 1400°K) for convenient measurement. 

Platinum was chosen for the metallic leads to 
the specimen because of its chemical inertness, low 
vapor pressure, and spectroscopic purity, except 
for a trace (<0.03 percent) of copper." 


18 G. A. Rankin and H. E. Merwin, J. Am. Chem. Soc. 38, 
568 (1916). They report, for example, that MgO and CaO do 
not form compounds. 

4 This platinum was prepared by Mr. H. T. Reeve of these 
laboratories. 
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Another chemical question is the effectiveness of 


methane in reducing (Ba, Sr)O under the specified 
conditions. Such reduction is shown to be thermo- 
dynamically reasonable by estimating the Gibbs 
free energy at 1273°K for the reaction, 


BaO(solid) +CH,= Ba(gas) +CO+2H2. 


This free energy turns out to be +55,000 calories.'® 
The corresponding equilibrium constant K _ is 
4xX10-"; furthermore K=(pBa-pco- pH,”)/pcHy,. 
This together with our idealized conditions that 
pBa=pco= phe yields pBa=3X10-*- pcu,', where 
the partial pressures p are identified by the chemical 
symbols in the subscripts and are expressed in 
atmospheres. The equilibrium pressure of barium 


vapor in the presence of 0.025-mm Hg of methane _ 


is thus calculated to be 0.2 mm, with an accuracy of 
about a factor of ten. This is about one-tenth of 
the vapor pressure of solid barium as extrapolated 
to 1273°K from Van Liempt’s data.'!® Hence equi- 
librium conditions should be favorable to the 
production of gaseous barium and its solution by 
solid BaO to the extent of 0.01 or more of its 
solubility limit at this temperature. Schriel’s work 
suggests a solubility limit of about one atom percent 
of Ba in BaO at about 1420°K."* 

A mass spectrometer analysis of the gaseous 
products of the reaction showed complete recovery 
as CO of the carbon in the CH, (with small amounts 
of CO, and unreacted CH,), as expected from the 
above. However, a control tester containing no 
(Ba, Sr)O yielded a similar complete recovery of 
methane carbon as CO, showing that other hot 
oxides such as MgO, Al.O;, or WOz compete with 
(Ba, Sr)O in reacting with methane. 

Another possibility is that methane treatment 
produced a second solid phase of conducting ma- 
terial deposited on or scattered through the 
(Ba, Sr)O. The complete recovery of the carbon in 
the mass spectrometer experiment indicates that 
no considerable deposit of carbon was formed. This 
is confirmed in Table IV, which shows the spectral 
emissivity measured concurrently with conductivity 
as discussed below. Methane treatment produced 
no net change of emissivity greater than experi- 
mental error. Small quantities of carbon or other 
good conductors would be expected to affect this 
property. Hence it is concluded that no second 
conducting phase is produced by the methane 
treatment. 

The chemical evidence then supports the con- 


6 The estimate is based on data from K. K. Kelly, U. S. 
Bur. Mines Bull. No. 371 (1934); No. 383 (1935); No. 394 
(1936) ; together with the following assumptions: 


(1) aun! high-temperature molar heat capacity of BaO equals that of 
(2) The heat of sublimation of barium at 298°K is +en8 ry calories 


which is the value observed at 800°K-1050°K az. ‘. M. Van 
Liempt, Rec. Trav. Chim. de Pays-Bas 55, 468 (1936 


6 M. Schriel, Zeits. f. anorg. allgem. Chemie 231, 313 3 11937). 
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Fic. 8. Probe voltage vs. current in (Ba, Sr)O before activation, 
at 1176°K. Tube C-2. 


clusions that in these experiments reasonably pure 
(Ba, Sr)O is available for reaction with similarly 
pure methane; that the significant effect of the 
methane is to reduce small quantities of (Ba, Sr)O, 
and that after the reaction the conducting specimen 
consists of (Ba, Sr)O which should contain at least 
a fraction of any excess barium that it is capable of 
dissolving. 


Measurement of Conductivity and Temperature in 
Conductivity Experiments 


The primary measurement of effective conduc- 
tivity of the porous (Ba, Sr)O specimen was the 
measurement of the potential difference between 
probes 2 and 3:(Fig. 2) as a function of the current 
through an activated specimen, at several tem- 
peratures. The behavior of such specimens was 
invariably ohmic over the entire range of observa- 
tions, which included current densities as high as 
65 ma/cm?. The resistance determined by a slope of 
such current-voltage measurements was combined 
with the dimensions of the specimen as originally 
sprayed to calculate effective conductivity. For 
this purpose the area of cross section was deter- 
mined by the weight of the coating, together with 
its length and its bulk density. This bulk density 
was found to be 1.1 g/cm?, as measured after 
depositing a test specimen according to the same 
spraying procedure on the flat metal block. These 
dimensions take no account of over-all shrinkage of 
the specimen during decomposition of the carbon- 
ates, nor of its obvious porosity. Nevertheless the 
data are reported as effective conductivity to 
facilitate comparison among the results recorded 
here and elsewhere in the literature. 

The conductivity was observed at successively 
increasing temperatures from about 1000°K to 
1350°K, with a final check at the initial temperature 
to detect any irreversible change occurring during 
the four hours involved. After activation there was 
no such change of more than 10 percent. 
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Tas_e V. Effect of helium on effective conductivity of 
activated (Ba, Sr)O. 











Speci- o g o 

men T In vacuum In He 1 atmos. Re-evacuated 
C3 13038°K 8.1X10%ohm™'cm™ 8.2X10%ohm™"cm™ 8.210 ohm™ cm- 
C-2 1304°K 3.9X10"% 4.1xX10°% 4.7X10% 

C-2 1066°K 86x10" 7.3X10 7.3X10 

C-l 1247°K 49X10 2.5X10% 3.8X107* 





While the conductivity of the active (Ba, Sr)O 
was thus invariably ohmic, that of the inactive 
state as measured by the null method was in- 
variably non-ohmic in the sense typified by Fig. 8. 
For BaO the current-voltage curves of Meyer and 
Schmidt? exhibit the same general behavior in the 
active and inactive states respectively. This ten- 
dency towards saturation of the current in the 
inactive state may reflect thermionic emission in 
addition to conduction current, as discussed in 
Section VI. The slope (dE/dJ)1~0 was taken to be 
the resistance in determining conductivity as a 
function of temperature as shown in the lowest 
curve of Fig. 3. This procedure is justified insofar 
as the results are used to show that methane treat- 
ment increases conductivity. 

An alternative method of measuring the re- 
sistance of the specimen lying between the probes 
2 and 3 of Fig. 2 was to connect these probes to a 
d.c. Wheatstone bridge. The results were found to 
agree closely with those of the null method for the 
activated state and to disagree extensively in the 
case of only one measurement of the inactive state. 
Hence the contribution of electrode contact and 
polarization effects is negligible in this experimental 
arrangement of a two-electrode system. 

Temperatures were measured with calibrated 
optical pyrometers. In order to correct these 
pyrometer readings in accordance with the spectral 
emissivity of the specimen, a reflectometer of the 
type developed by Prescott and Morrison” was 
installed. The use of this technique with the 
translucent (Ba, Sr)O, supported by a translucent 
ceramic containing a heater, raises several questions. 
The method relies on the determination of the 
spectral emissivity E, from the measured reflec- 
tivity R, according to 


E,+R,=1. (8) 


This in turn requires that only a negligible fraction 
of the light incident on the tester should be trans- 
mitted through it.'? A rough measurement of this 
fraction for the oxide-coated magnesia yielded a 
value of <10-%, justifying the use of Eq. (8). 

The spectral emissivity was measured concur- 


17C, H. Prescott, Jr., and James Morrison, Rev. Sci. Inst- 
10, 36 (1939). 

18 American Institute of Physics, Temperature, Its Measure- 
ment and Control in Science and Industry (Reinhold Publishing 
Corporation, New York, 1941), Chapter 1. 
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rently with the conductivity throughout these 
experiments. The values obtained at the wave- 
length of 0.66u are recorded in Table IV for three 
specimens. In all cases these results were inde- 
pendent of temperature within the precision of the 
measurement; hence the table records mean values 
and their mean deviations over the indicated tem- 
perature ranges. The variation of £, among these 
three specimens may reflect differences of granu- 
larity. The mean values of Table 1V were used to 
determine the true temperature T from the pyrom- 
eter reading D in the standard way, according to 


1/T—1/D=(A/c2) Ink, (9) 


where cz was taken to be 14,380 micron-degrees. 
These corrections were so large that the uncertainty 
in T associated with that in &, is of the order 
of +25°. 

The correctness of these temperature measure- 
ments involves the assumption that the light 
employed in temperature measurement originates 
primarily in the (Ba,Sr)O rather than in the 
translucent supporting ceramic or the tungsten 
heater. This was tested by several measurements 
of the conductivity of activated (Ba, Sr)O sur- 
rounded by helium at a pressure of one atmosphere ; 
such measurements were compared with observa- 


‘tions at the same pyrometer temperature in 


vacuum. The results in Table V show remarkable 
agreement except for specimen C-1, where a par- 
tially irreversible change is ascribed to chemical 
effects. 

If the assumption concerning the origin of the 
light is valid, the presence of an atmosphere of 
helium should make no difference, but if it is not, 
the thermal conductivity of the helium should act 
to increase the true temperature, and thus the 
electrical conductivity of (Ba, Sr)O compared with 
that in vacuum for a given pyrometer reading. 
Thus in a typical case at a pyrometer temperature 
of 1190°K (T=1304°K), the heater temperature as 
estimated from its resistance was 1540°K in vacuum 
and 1350° in helium. The intensity of 0.66y light 
received from the heater and ceramic was certainly 
less in helium than in vacuum. Hence if our assump- 
tion was invalid, more light was received from the 
(Ba, Sr)O in helium than in vacuum, at this 
pyrometer reading. This would then require that 
the temperature and conductivity of (Ba, Sr)O be 
higher at a given pyrometer reading in helium than 
in vacuum. The observed agreement thus indicated 
that within the limits of error of the experiments 
and their interpretation, the primary origin of the 
light lay in the oxide, except for the possibility of 
error in the above data due to an improbable com- 
pensation of spurious effects. For example, if a large 
portion of the current were due to thermionic 
emission around the specimen or through its pores, 


JOURNAL OF APPLIED PHYSICS 














\y 


\¥ 


we 


ure 
as 
um 
ght 
nly 
np- 
the 
this 
hat 
) be 
han 
ited 
ants 
the 
y of 
om- 
urge 
onic 
res, 


SICS 








such compensation is conceivable. This latter effect 
will be considered in detail in Section VI, as it is of 
fundamental importance in all of the present work. 


Techniques in Conductivity— 
Emission Measurements 


After winding the platinum wires on the MgO 
ceramic, the structure was cleaned with solvents, 
vacuum and hydrogen-treated, and sprayed with the 
carbonate mixture to a thickness of 0.005 inch over 
the ceramic. Care was taken during the spraying 
process to pack the carbonate down well between 
the individual wires, the first part of the spraying 
as “‘dry’’ as possible, with the final layers on top 
of the wires being applied in the usual ‘‘wet”’ spray. 

The tube, when mounted on the vacuum system, 
was given a preliminary 2-hour bake-out at 400°C. 
The anode and support collars were then outgassed 
at 900°C, by high frequency induction heating. 
The carbonates were broken down thermally, at 
temperatures up to 900°C, until the pressure in 
the system dropped below 5X10-’-mm Hg. 

Temperatures were measured with an optical 
pyrometer, although once the relationship between 
pyrometer temperature and heater power was 
established for a given tube, the latter was taken 
as the more reliably reproducible indication of tem- 
perature. Corrections to true temperature were 
made by using the spectral emissivity data obtained 
in the conductivity work, the emissivity being 
taken as 0.15. Due to the difficulty of investigating 
with certainty the possibility of light originating in 
the heater-ceramic system, it cannot be stated 
definitely that the true temperature of the coating 
was not below that reported, however. 

The conductivities were measured using the two 
interlacing platinum wires as a two-electrode sys- 
tem. Justification for this lay in the results of 
Fig. 7, and in the conductivity work, where it was 
found that 4-electrode potentiometric measure- 
ments agreed with 2-electrode bridge measurements 
in every case where platinum electrodes were used, 
for that geometry. 

The conductivities were all obtained from the 
slopes of current-voltage curves taken between the 
two leads, several points being taken on both sides 
of the origin at a low enough voltage gradient to 
give ohmic results. It was observed that non-ohmic 
behavior of the type shown in Fig. 8 usually set in 
at about 0.5-1.5 volts between the leads; the slopes 
were all taken from linear current-voltage measure- 
ments in the region up to 0.05-v maximum between 
the leads in the higher activity ranges, 0.1—-0.2-v in 
the. low activity ranges. The observed conductances 
were converted into apparent conductivities for the 
porous oxide using the best available geometrical 
data for the cathode structure. 
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Thermionic emission was measured by a pulse 
technique, single 1/120-second sine-shaped con- 
tinuously variable voltage pulses being applied 
between the plate and the two core wires, which 
were connected together externally for this measure- 
ment. A signal proportional to the current through 
the tube was applied to the vertical axis of a 
cathode-ray oscilloscope, and a signal proportional 
to the three-halves power of the voltage developed 
across the tube was applied simultaneously to the 
horizontal axis. Thus, the oscilloscope trace gave a 
complete 3/2 power plot, which provided a check 
at all times as to possible fading, sparking, and 
gas effects, with a straight line oscilloscope trace as 


‘long as the emission was space-charge limited. The 


point at which the current fell 20 percent below the 
space-charge line was taken as a measure of the 
emission as shown in Fig. 9, the method being 
largely arbitrary, as are most emission measure- 
ments. The voltage was increased in each case 
above that just necessary to give 20 percent devia- 
tion, to observe the effect on the current at the 20 
percent deviation point. This usually increased 
slightly as the voltage rose, approaching a more or 
less constant value, which was the current actually 
used as a measure of the activity of the cathode. 
The crossing af the 20 percent line took place after 
about the first two milliseconds of the 1/120 second 
pulse, several single pulses being applied within 
about a minute to obtain the measurement. 


VI. MECHANISM OF CHARGE TRANSFER 


The conductivity involved in the discussion of 
Section II is the bulk electronic conductivity of the 
oxide semi-conductor, and the thermionic emission 
is assumed to be that of electrons emitted from the 
surface layers after transfer through the solid from 
the core metal to the surface layers. It is, therefore, 


CHARGE LINE 
TRAPOLATED 


/ 


OBSERVED 
Ip VS Ep372 
CHARACTERIST 


20 % DEVIATION 
POINT TAKEN AS 
MEASURE OF 
THERMIONIC 
EMISSION 


20% DEVIATION LINE 
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HORIZONTAL DEFLECTION,~ Ep 


Fic. 9. Plot of oscilloscope trace showing method used in 
measuring thermionic emission. 
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TABLE VI. 


Conductivity Emission 

(ohm~! cm~!) (ma/cm?) 
Before CH, treatment 8.8x 10-5 4.7 
After CH, treatment 8.8 10-5 <0.01 
Before O» 1.8 10-4 9.5 
After O2 (room temp. exposure) 1.2x10~¢ <0.01 
Before O2 2.910 15.0 
After O, (room temp. exposure) 1.610 0.02 


a vital point that this be the actual mechanism of 
charge transfer. The question of the method of 
charge transfer through the oxide is of fundamental 
importance to any theory of the behavior of oxide 
cathodes, and is one which is treated inadequately 
in the experimental literature. The two most 
possible sources of error in conductivity to be 
considered are (1) that there was an ionic con- 
ductivity to a sufficient extent to obscure the 
results, and (2) that it was not a bulk property of 
the solid being measured, but rather electron trans- 
fer by thermionic emission through the pores of the 
material or around the material, conductivity along 
the pore walls, high resistance contacts between 
individual crystals, or a property of the metal- 
oxide interface. ° 

Concerning the first of these, attempts to detect 
any polarization by passing a high d.c. current 
through the coating for long times, failed to show 
any detectable reverse current when the applied 
voltage was suddenly shut off. The conductivity 
behaved like excess electronic semi-conduction in- 
sofar as it was increased by reducing agents and 
decreased by oxidizing agents. The indication is 
that the ionic contribution to the conductivity was 
small. 

Concerning the possibility that the electron trans- 
fer was not really a bulk property of the solid, it is 
considered unlikely that it was a property of the 
metal-oxide interface. In the case of the conduc- 
tivity measurements, the agreement in every case 
where platinum electrodes were used between 4-elec- 
trode potentiometric measurements and 2-electrode 
bridge measurements, indicates there was no sig- 
nificant interface resistance. In a few experiments 
using Grade A nickel electrodes, on the other hand, 
it was found that these measurements no longer 
agreed after oxidation treatment, and that con- 
ductivities were not reproducible above 1000°C. 
The absence of appreciable interface resistance in 
the conductivity tester, and the agreement of the 
conductivity vs. 1/T curves for this structure with 
that for conductivity-emission structure, are con- 
sidered to justify the use of 2-electrode measure- 
ments in the latter case. 

The possibility that thermionic emission current 
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in parallel or in series with conduction current in 
the solid played an important role must be con- 
sidered in detail. The specimens studied here were 
porous, and exhibited fine cracks extending through 
at least part of the thickness of the specimen, as is 
usual with oxide coatings. Except for the experi- 
ments of Becker and Sears at much lower tempera- 
tures,” little consideration has been: given to such 
spurious effects. A number of experiments were 
carried out in the present work to investigate this 
possibility, not all of which were conclusive. The 
approach has been to subject the specimen to 
treatments which would affect thermionic emission 
current but not conduction in the solid, and to look 
for changes in the observed conductivity during 
such treatment. Three methods of accomplishing 
this have been the introduction of an atmosphere 
of helium, the application of a magnetic field, and 
the deactivation of surface layers of the oxide 
crystals by proper chemical treatment. 

In the helium experiments, vacuum measure- 
ments were compared with conductivities obtained 
in 1: atmos. of helium, the specimen being heated to 
the same pyrometer temperature. Since the mean 
free_ path of a thermal electron in helium can be 
estimated from the literature'® to be no more than 
1<10-* cm at 760 mm and 1000—1300°K, and since 
most of the volume of the bulk specimen is visibly 
occupied by voids of much larger dimensions, this 
test should detect any important spurious effect. 
The results of such experiments on conductivity 
testers in the activated state are shown in Table V. 
Before admission, the helium was purified with 
copper at 400°C, and anhydrous CaSQ,; mass spec- 
trometer analysis showed only 1.4 percent N» and 
0.1 percent Hy. The effect of the gas heating mica 
and glass surfaces resulting in transfer of con- 
taminants to the specimen appears to have been 
small, in view of the reproducibility obtained. 

It was then observed that a magnetic field of 
2000 gauss reduced the current in a conductivity 
tester in the inactive state by as much as 50 percent, 
the effect being reduced after activation to the 
range 0-20 percent. In no case in scores of tests 
made on three conductivity specimens in various 
states of activation, at temperatures ranging from 
1000—1300°K, were effects of more than 1 percent 
found to be produced by a magnetic field with an 
atmosphere of helium present. Such small changes 
could be expected from the effect of the field on 
the geometry of the tube, including the heater coil. 
Magnetic effects in vacuum of a similar magnitude 
were observed in conductivity-emission testers, 
also. 

The third type of evidence bearing on the ques- 
tion is that shown in Fig. 6 and described in Sec- 


1? R. B. Brode, Rev. Mod. Phys. 5, 257 (1933). 
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tion IV, in particular with regard to the experiments 
involving deposition of carbon, and room tempera- 
ture oxygen deactivation. Three such experiments 
on conductivity-emission testers are shown in 
Table VI. 

If the electron transfer through the solid were 
due primarily to thermionic emission currents in 
the pores and cracks, rather than electronic semi- 
conduction through the solid with subsequent emis- 
sion from the outer oxide surfaces in the case of 
thermionic emission current to the anode, the 
conductivity should have been decreased in the 
same way as the emission current. In the case of 
the oxygen deactivations, the small change in the 
conductivity may be ascribed to the fact that as 
the temperature was raised to 700°C for the meas- 
urements, diffusion would become more rapid, and 
the effect of the deactivated surface layers would 
begin to appear throughout the body of the oxide. 

It appears established from these data that in 
the activated state conductivities, as measured in 
these specimens, reflect primarily the flow of elec- 
trons through the solid; in the inactive state 
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spurious thermionic currents are of more impor- 
tance. However, it is concluded that with an error 
no greater than that of the emission measurements, 
the conductivity measurements reflect solid con- 
duction, in view of the following: 


(1) In the activated state this current was not affected 
by the presence of helium at one atmosphere 
pressure. 

(2) Processes designed to affect the surface layers through- 
out the porous specimens greatly reduced the 
emission to the anode but did not affect the con- 
ductivity. 

(3) Conductivity measurements made in specimens of 


widely variant geometry agree within reasonable 
error. 


(4) The current-voltage curves in the conductivity- 
emission testers were ohmic up to a voltage 
gradient of the order of 50 v/cm over three orders 
of magnitude in the conductivity. 
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To obtain suitable conductivity for direct heating, tungsten powder is added to sintered thoria 
cathodes. The thermionic properties of a cathode, sintered from a mixture of 67 percent thoria and 33 
percent tungsten, were investigated. Emission constants were determined and the emission was found 
to be somewhat lower than that of cathodes sintered from pure thoria. 

Change of cathode activity with temperature was observed and studied in detail. The cathode 
showed high activity at temperatures which would bring about deactivation of ordinary thoriated 
tungsten and decreased its activity when the temperature was lowered to the range where thoriated 
tungsten would be activated. A theory is proposed for the mechanism of the change in activity, in- 
volving the production of free thorium due to reduction of thoria by tungsten, diffusion and evapora- 
tion of thorium from thoria. Evaporation of thorium as well as thoria was proved by the use of a 
tungsten filament probe which detected evaporation of materials from the cathode by the change of its 
own emission. The energy of evaporation of free thorium from thoria is estiniated from the experiments 
to be 46600 cal./g-atom. The energy of evaporation of thoria is found to be 184000 cal./g-molecule. 


Thermionic Emission from Sintered Cathode of Thoria and Tungsten Mixture* 


H. Y. Fan** 
Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received November 15, 1948) 





The behavior of the cathode is satisfactorily explained by the proposed theory. 





INTRODUCTION 


INTERED thoria cathodes, developed in recent 

years, possess a number of desirable properties, ! 
of which the two most important are: (1) The 
ability to withstand severe punishment caused by 
sparking, bombardment and other destructive phe- 
nomena. (2) Operation at higher temperatures, as 
compared with ordinary oxide-coated cathodes, 
with consequent ability to dissipate extra heating 
power, at the same time quite a large density of 
electron emission is obtainable. These properties 
make this type of cathode especially suitable for 
applications such as high powered magnetrons. 
However, due to the extremely high resistivity of 
thoria it is necessary to provide such cathodes with 
a separate heater which has to be run at very high 
temperatures. The design of such heaters is a 
problem and burn out of the heater is frequently a 
source of trouble. 

To simplify the problem of heating, cathodes are 
now being developed which are sintered from thoria 
mixed with tungsten powder.? By varying the per- 
centage of added tungsten the resistivity of the 
cathode can be adjusted to a value suitable for 
direct heating. The cathode studied here consists 
of 67 percent thoria and 33 percent tungsten by 
weight. It is made in the form of a hollow cylinder 
pressed from such mixture and fired in a hydrogen 
furnace at 1800°C for one half hour. The dimensions 
of the cathode are: 1” in length, 0.265” outside 
diameter and 0.192” inside diameter. Molybdenum 


* This work has been supported in part by the Signal Corps, 
the Air Materiel Command, and ONR 


** On leave of absence from National Tsing Hua University, 


Peiping, China. Now at Department of Physics, Purdue Uni- 
versity, Lafayette, Indiana. 

1M. A. Pomerantz, N.D.R.C. Division 14, Report 517 
(1945). 

2 This development seems to have been originated by the 
Raytheon Manufacturing Company. 
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sleeves platinum brazed to the two ends serve as 
current leads. Figure 1 shows the resistivity p and 
the required heating power W as function of tem- 
perature. The resistivity of tungsten p», multiplied 
by a factor 845 to adjust it to the same value as 
the resistivity of the cathode at room temperature, 
is shown for comparison. During the operation of 
over five hundred hours at temperatures varying 
from 1100°K to 2300°K the resistance of the 
cathode, as measured between the cathode leads 
of the tube in which the cathode is mounted, in- 
creased by approximately 10 percent. The results 
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Fic. 1. Heating power and resistivity of the cathode as 
functions of temperature. 
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Fic. 2. Diagram of tube structure and measurement circuit. 


of the study of the properties of such a cathode are 
reported here. 


EXPERIMENTAL ARRANGEMENT 


A conventional diode structure is used for the 
study of thermionic emission (Fig. 2). The anode is 
a cylinder, concentric with the cathode, 1” in 
diameter and ;%”’ long. It is provided with guard 
rings, each of 1” j in diameter and 3” long. For the 
measurement of true temperature a hole, 10 mil in 
diameter, is drilled through the wall in the middle 
of the cathode. A larger hole is made in the anode 
so that the hole in the cathode wall can be viewed 
conveniently with an optical pyrometer. Before 
sealing off the pump the vacuum as read by an 
ionization gauge was below 10-7 mm Hg with the 
cathode at a high operating temperature and with 
a d.c. emission giving an anode dissipation several 
times the highest value ever reached during the 
operation after seal-off. The tube is also gettered. 
The probe arrangement shown by the dotted lines 
was added in a second tube only. It will be de- 
scribed later. 

D.c. emission was measured for low temperatures 
of the cathode. For higher temperatures a pulse 
modulator was used which gives pulse lengths from 
0.5 to 10 usec. The pulse current and voltage were 
measured on the screen of a synchroscope. The 
circuit is shown schematically in Fig. 2. The cathode 
was heated from the a.c. line through a Sola con- 
stant voltage transformer, a variac and a high 
insulation transformer in series. The last mentioned 
was used because the cathode was raised to high 
negative voltages during the pulses. 


EMISSION IN ACCELERATING FIELDS 


Figure 3 shows the curves of J as function of V?. 
It is seen that for sufficiently high fields the curves 
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approach the theoretical Schottky slopes. Since the 
surface of the cathode consisting of thoria and 
tungsten mixtures is nonuniform, deviation from 
the theoretical slopes at low fields is expected on 
the basis of patch theory.* The fact that the curves 
approach the theoretical slopes more readily for 
higher temperatures is also consistent with the 
predictions of the patch theory. The sharp bending 
away from the theoretical line in the curve for 
2040°K is the transition into space charge limited 
region; below the bend the curve follows closely 
the calculated $ power curve. 

Since the cathode is directly heated its zero field 
emission cannot be measured readily. As the emis- 


_sion approaches Schottky lines for high voltages we 


may well use the emission values at a definite, 
sufficiently high voltage for the determination of 
emission constants. All emissions were accordingly 
measured at 2500 v. A still higher voltage would be 
preferable, as the emissions at low temperatures 
have not quite reached the theoretical slopes at 
this voltage. But we were limited in this respect by 
the pulse modulator. It can be seen from Fig. 3 
that a higher voltage will not make a great dif- 
ference. 

The emission under applied d.c. voltage at low 
temperatures and under pulse voltage at high tem- 
peratures maintained a definite value and showed 
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Fic. 3. Schottky ~ of emission current. Current scales for 
different curves: 1=2.0 amp. (2040°K); 1=0.126 amp. 
(1770°K); 1=10 Te (1660°K); 1=10 ma (1650°K); 1=0.20 
ma (1420°K); d.c. and pulsed emission curves are not for the 
same state of the cathode. 


3 J. A. Becker, Rev. Mod. Phys. 7, 95 (1935). 
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Fic. 4. Richardson plot for equilibrium state of the cathode at 
two different temperatures. 


no sign of decay. For temperatures where both d.c. 
and pulsed emissions could be measured con- 
veniently the values were the same. This result 
shows that there is no fundamental difference be- 
tween emissions under d.c. and pulsed conditions. 
Furthermore the tube was opened to air and re- 
built twice. The emission properties of the cathode 
showed no marked change. This is indeed one of the 
desirable properties of this type of cathode. If 
however, the temperature was changed and main- 
tained at a different value the emission underwent 
a time change irrespective of applied voltage. This 
indicates a change of the state of the cathode with 
temperature. This effect will be discussed more fully 
later on. For determining the emission constants the 
cathode was maintained at a definite temperature 
T> and the emissions at different temperature 7;’s 
were taken to be the initial values when the tem- 
perature was brought from J>) to T;. Thus all 
emissions corresponded to a definite state of the 
cathode at 7». Thoria being a semi-conductor, its 
emission equation 

4=ATe*o/kT (1) 


may have different possible values of m according 
to’ theory, depending upon the type of semi-con- 
ductor it is. Without knowing what type of semi- 
conductor the thoria in the cathode is, there is no 
basis for taking any particular value for m. We shall 
use the conventional value »=2 for convenient 
comparison of the emission constants A and ¢g with 
data given by other authors. Figure 4 shows curves 
of log(t/T?) versus 1/T for two values of JT): 2040°K 
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and 2320°K. Both g and A are different for the 
cathode states at these two temperatures. ‘The 
state at 2320°K gives higher emission on account of 
higher value of A. 

Wright’ gives for thoria coated on tantalum 
g=2.5(4) ev, A=2.5 amp./cm? deg.” for d.c. and 
¢g =2.6(2) ev, A=7.5 amp./cm? deg.? under pulsed 
conditions. He gives also g=2.5(5) ev, A=3.5 
amp./cm? deg.? for the pulsed emission from an 
extruded thoria sleeve heated by a tungsten spiral. 
Hanley® gives for cataphoretically coated thoria 
cathodes average values: g=2.67 ev, A=2.63 
amp./cm? deg.? under d.c. conditions and g=2.55 
ev, A=5.62 amp./cm? deg.? under pulsed condi- 
tions. Our values of ¢ are higher than those given 
by these authors. This result may be due to the 
fact that a part of the surface of our cathode is 
presented by tungsten which has a higher work 
function. Dependence of emission constants on the 
cathode temperature was also observed by Wein- 
reich® for thoria coated on tungsten. Contrary to 
our results the constant A remains unchanged while 


the work function is lower for higher temperature 
of the cathode. 


TEMPERATURE DEPENDENCE OF ACTIVITY 


It has been mentioned that the cathode activity 
shows temperature dependence. In fact the emission 
constants have different values for different condi- 
tioning temperatures of the cathode. The observa- 
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Fic. 5. Reversible change of activity with cathode temperature. 


4D. A. Wright, Nature 160, 129 (1947). 
5 T. E. Hanley, J. App. Phys. 19, 583 (1948). 
6 M. O. Weinreich, Rev. Gen. Elec. 54, 243 (1945). 
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tions on this effect will now be given in detail. If the 
cathode was kept at a temperature 7) until its 
emission had reached a steady value then brought 
down to a lower value T2, the emission tested inter- 
mittently at T, changed with time. If the tempera- 
ture was again restored to 7), the emission would 
build up gradually to its original equilibrium value. 
The process could then be repeated again. This 
reversible change is shown in Fig. 5 for 7; = 2040°K 
and T,=1670°K. It must be emphasized that the 
anode voltage was applied only when a point was 
taken and only long enough for the measurement 
to be made. This is, therefore, an entirely different 
phenomenon from the decay of emission under 
applied voltage observed for oxide-coated cathodes. 
Figure 6 shows the good reproducibility of such 
curves taken at T,~1670°K for two values of 7): 
2320°K and 2040°K. Curves taken many days 
apart, during which time the cathode was subjected 
to different operating conditions, would check each 
other if only the cathode was initially brought to 
equilibrium conditions at the same 7). This fact 
indicates that it was not the effect of residual gases 


in the tube. In a gettered tube gas effect cannot be’ 


expected to show such good reproducibility. Fur- 
thermore, it was found that such a curve could be 
interrupted for a long time, overnight, with the 
cathode cold; the emission would start again at 
the same value where the curve was interrupted 
and keep on changing in the same way as if there 
had been no interruption. Curves were taken at 
different JT, values ranging from 1200°K to 1920°K 
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Fic. 6. Variation of emission with time after a decrease in 
cathode temperature. 
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Fic. 7. Variation of emission with time after a decrease in 
cathode temperature. 


for the same two values of T;. Figures 7 and 8 give 
the curves for two other values of T;. They show 
the same behavior. Comparing Figs. 6, 7 and 8 
we see that with decreasing temperature the 
rate of change of emission decreases. All curves 
for 7,=2040°K show initial rise, whereas for 
T, = 2340°K they start off rather flat. Generally for 
the test temperature range 1200°K to 1920°K 
curves are similar in character. At very low testing 
temperatures, 7;=1210°K, the initial rise in curves 
for T7;=2040°K could still be observed. The emis- 
sion rose in 6 or 7 hours to a rather steady value, 
thereafter changed little in the course of 24 hours. 
The drop in emission observed at higher 72 had 
either disappeared or proceeded at a very slow rate. 

To follow the change in emission constants data 
for curves of log(i/T?) versus 1/T were taken 
quickly at four points A, B, C and D of the run 
shown by curve 2, Fig. 8. These curves and the 
values of emission constants they yield are given 
in Fig. 9. They show that both g and A underwent 
changes. 

In a cathode with thoria and tungsten present to- 
gether such reversible and reproducible changes of 
activity can naturally be ascribed to the production, 
diffusion and evaporation of free thorium. The 
effect of free thorium on the activity can be either 
a bulk effect or a surface effect. When dispersed 
inside thoria crystals free thorium may change ¢ 
and A of thoria as impurities in a semi-conductor. 
It may also change ¢g and A by being absorbed on 
the surface of the cathode. It may be argued that 
free thorium must be produced during the operation 
of the cathode, since if we have only a fixed amount 
of thorium as produced in the cathode originally 
during its fabrication, then without evaporation an 
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equilibrium condition should soon be established 
and there would then be no changes of activity. 
Evaporation should cause irreversible changes in 
the activity which are not observed. 

The results of the following experiment, shown 
in Fig. 10, are significant for understanding the 
mechanism of the change of activity. After the 
emission reached a steady value at 7,=2040°K, 
the temperature of the cathode was lowered to 
T:=1670°K and changes of emission with time 
were observed. After an interval of time the cathode 
temperature was increased to 7; for a very short 
time and brought back to 72 again. In one of the 
curves the short time elevation to T; was made at 
the time when the emission at 1679°K had reached 
its peak. In the second curve this was made when 
the emission at 1670°K had passed beyond the peak. 
In both cases immediately after the short time at 
the higher temperature the emission was greatly 
reduced. It then partially recovered. 

In the light of all the experiments we have 
arrived at the following picture of the mechanism 
of the cathode: (1) The variation of the activity is 
caused mainly by the variation of the amount of 
free thorium adsorbed on the surface. (2) Free 
thorium is produced through reduction of thoria by 
tungsten throughout the cathode. (3) The heat of 
reduction of thoria by tungsten is higher than the 
heat of evaporation of thorium from the cathode, 
which is in turn higher than the heat of diffusion 
of thorium in the cathode. 

According to this hypothesis, when the tempera- 
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Fic. 8. Variation of emission with time after a decrease in 
cathode temperature. 


686 





*3.00ev A96 
*2.9lev A008 
*3.26ev A616 
*3.40ev A*69.) 


oo @p 
+t fF 


\AA 











nN 1%10 he's aad a 






































os 
1000 
=r 


Fic. 9. Richardson plots for different states of the cathode. 


ture of the cathode is suddenly raised the amount of 
absorbed thorium will at first decrease on account 
of more rapid increase of evaporation than diffusion. 


However, the increased rate of thorium production. 


more than compensates the increased rate of evapo- 
ration. Therefore, if the cathode is held at the higher 
temperature, eventually a higher thorium concen- 
tration will be built up, which in equilibrium sup- 
ports a larger number of thorium atoms on the 
surface. The experimental facts: (1) Increasing the 
cathode temperature for a short time decreases the 
emission (Fig. 10). (2) After a decrease in cathode 
temperature there is an initial rise of emission with 
time (curves 2, Figs. 6, 7 and 8), and (3) Equi- 
librium emission activity is higher for higher tem- 
perature (compare initial emission values of curves 1 
and 2, Figs. 6, 7 and 8), can thus be all explained. 
The fact that curves 1, Figs. 6, 7 and 8, do not show 
initial rise in emission, we assume to indicate that 
at 7,=2320°K the amount of adsorbed thorium is 
close to a full monolayer. 

Although in thoriated tungsten filaments con- 
siderable production of thorium takes place at 
temperatures higher than the temperatures reached 
by our cathode, the production should be much 
more rapid in our cathode than in a thoriated 
tungsten filament at the same temperature. Where- 
as a thoriated tungsten filament contains about 
1 percent thoria our cathode is a mixture of 33 
percent tungsten and 77 percent thoria. 

The evaporation of thorium from tungsten is 
very slow at the temperatures at which changes of 
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Fic. 10. Effect of short-time-elevation of the 
cathode temperature. 


activity of our cathode were observed as shown by 
Figs. 6, 7, and 8. Also the heat of evaporation of 
thorium from tungsten is higher than the heat of 
reduction of thoria. But in our case thorium may be 
evaporated from surfaces of thoria rather than from 
tungsten. The data on thoriated tungsten do not 
necessarily contradict our hypothesis. However to 
substantiate our hypothesis an independent proof 
of evaporation of thorium during the change of 
cathode activity is desirable. The following experi- 
ment was made for this purpose. 


PROBE MEASUREMENT 


The tube was rebuilt with the same cathode. 
The structure is the same as described previously 
except for the addition of a probe. Its arrangement 
is shown by the dotted lines in Fig. 2. It is a hairpin 
filament of 10 mil tungsten wire. Its temperature is 
calculated from the heating current and its emis- 
sion measured at 1390°K, 1550°K and 1760°K 
checks with the calculated values. It serves to 
detect evaporation of materials from the cathode, 
which condensing on it changes its emission. Its 
emission could always be brought back to normal 
values by flashing it at 2460°K for a sufficient 
period of time. A few seconds of flashing was 
enough to bring the emission almost back to normal 
values. It was found necessary sometimes to keep 
the probe at the flashing temperature for many 
minutes before its emission at the low testing tem- 
perature, 1390°K, was quite restored to its normal 
value. This result is probably due to the fact that 


VOLUME 20, JULY, 1949 


the top and the base of the hairpin are not far away 
enough from the cathode. The small amounts of 
material deposited on the tips of the lead wires, 
the spring, and the short lengths of filament near 
the contacts with the wires and the spring were not 
easily driven off and kept to diffuse on to the fila- 
ment. In testing the probe emission, the cathode, 
the anode and the guard rings were tied together 
to serve as the anode. The emission was measured 
with a galvanometer the sensitivity of which is 
4.4-10-” amp./mm. 


EVAPORATION OF THORIUM 


In the following experiments the probe was first 
cleaned by being kept at 2460°K for several minutes 
(with cathode at room temperature) until its 
emission at the testing temperature 7=1250°K 
gave the normal value of the clean surface. It was 
found that after glowing the cathode for six hours 
at 1770°K the probe emission at the testing tem- 
perature would not change appreciably, if the 
cathode had been previously brought to equi- 
librium at 1770°K. However, if the cathode had 
been previously heated at 2220°K, glowing the 
cathode for two hours at 1770°K increased the 
probe emission 20,000 times. These experiments 
prove conclusively that the increase of probe emis- 
sion in the second case must be due to free thorium 
evaporated from the cathode. The only other material 
that could come from the cathode to enhance the 
probe emission is thoria itself. But the evaporation 
of thoria depends upon the cathode temperature 
only and should not show such decided dependence 
on the previous treatment of the cathode. The 
conclusion is that considerable reduction of thoria 
takes place in the cathode at 2220°K building up a 
concentration of free thorium which keeps on 
evaporating at 1770°K. This is then a direct evi- 
dence of the correctness of our picture of the 
mechanism in the cathode. 
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Fic. 11. Probe emission versus time of glowing the cathode, 
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Fic. 12. [d(logi)/dt] versus time, after the cathode temperature 
is brought down from a high value. 


Figure 11 shows the probe emission at 1250°K 
as function of the time interval of glowing the 
cathode. Points taken for four different glowing 
temperatures of the cathode can be fitted by a 
single curve which is shifted horizontally a suitable 
distance for each cathode temperature. From these 
curves we get relative time intervals ¢ 9 necessary to 
glow the cathode at the different temperatures to 
obtain a certain value of the probe emission. Thus 
1/to is the relative rate of thorium evaporation. 

The rate of thorium evaporation can also be 
estimated from other experiments. Assuming the 
rate of evaporation to be proportional to the 
number of atoms on the surface we have 


dN/dt=DG—EN (2) 


‘ where N is the number of thorium atoms per cm? of 


the surface, G is the concentration gradient of free 
thorium at the surface, D is the diffusion constant, 
and £ is the constant for evaporation. It can be 
assumed : 

D=Dyer(#alk?, 

E= Eye-(2elk?), (3) 


In the case of thorium adsorbed on tungsten, it has 
been shown that? 


log(t/to) = L(N/N-), (4) 


where JN, is the number of adsorbed thorium atoms 
per cm? when the surface is completely covered, 


7W. H. Brattain and J. A. Becker, Phys. Rev. 43, 428 
(1933): 
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Fic. 13. Determination of the heat of evaporation of 
free thorium. 


i and i, are the emission currents corresponding to NV 
and N, respectively. The factor L depends on N but 
may be taken as constant for limited ranges of 
N/N.. By analogy we shall assume 


logit « (N+ const.). (5) 


When the cathode is brought down from a higher 
temperature to a temperature at which the reduc- 
tion of thoria is negligible, the concentration of free 
thorium becomes depleted with time. Eventually G 
becomes small and the evaporation term in Eq. (2) 
becomes predominant over the diffusion term: 


dN/dt~—EN. (6) 
In view of Eq. (5) we have 
d(logi) dN 





= ae-*, (7) 
dt dt 
d(logz) 
- log r —=const.— Et. (8) 
t 


Figure 12 shows curves of [d(logi)/dt] versus ¢ in 
semilogarithmic plot. They are calculated from the 
experimental curves of Fig. 7. When the experi- 
mental curves begin to fall quickly (according to 
our interpretation, thorium concentration has been 
depleted enough for the evaporation to predomi- 
nate) these plots do become straight lines. Accord- 
ing to our theory the slope of such straight line 
being equal to E should be determined by 72 
independent of 7}. 

Figure 13 gives logE and log1/ty (from Fig. 11) 
plotted against 1/7. The two sets of poiats logE 
and logi/t can be shifted vertically with respect to 
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‘ each other. All the points should fall on a straight 


line according to Eq. (3). The straight line drawn to 
fit these points is given by 


log io = const. + (10200/7). (10) 


The slope of this line multiplied by 4.575 gives the 
heat of evaporation to be 46600 calories per g-atom. 
Although this slope is not determined very accu- 
rately due to the scattering of the points, it serves 
to indicate the order of magnitude. In view of the 
fact that the points were determined from experi- 
ments made many days apart and the probe 
measurements were made with a different tube, 
the scattering may be considered not too bad. The 
fact that logi/t: and logE derived from different 
kinds of experiments can be made to fall around the 
same line lends support to our theory. That the 
heat of evaporation obtained is much lower than 
the value for evaporation of thorium from tungsten 
(4.575 X 44500 cal./g-atom)® indicates that thorium 
is evaporated from the surface of thoria. 


EVAPORATION OF THORIA 


At cathode temperatures above 1900°K the fol- 
lowing phenomena were observed in the probe 
emission. After the cathode had been glowed for a 
certain period of time the probe emission measured 
at 1390°K was above its normal value. When the 
probe temperature was raised to 1550°K the emis- 
sion decayed fairly quickly and might even drop 
below the normal emission. But flashing for a few 
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Fic. 14. Variation of probe emission after a period of heating 
the cathode. Emission measured at 1390°K, 1550°K, and at 
1390°K while being kept at 1970°K. 


‘I. Langmuir, Phys. Rev. 22, 357 (1923). 
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seconds at 1970°K made the emission tested at 
1390°K rise again. With successive flashing the 
emission went through a peak then fell approaching 
its normal value. This is shown in Fig. 14. The 
measurements were taken after the cathode was 
glowed at 2130°K for 5 minutes with the probe 
unheated. 

The fall of emission at 1550°K is apparently due 
to rearrangement by diffusion of materials deposited 
on the probe. The rise of the emission during the 
flashing of the probe we assume to be due to 
reduction of thoria deposited on the probe from 
the cathode. Figure 15 shows the probe emission, 
measured at 1390°K while the probe was flashed at 
1970°K, after the cathodes had been glowed for 
different lengths of time at 2030°K. The bottom line 
corresponds to about ten times the normal emission. 
The time required for the emission to fall to this 
value is seen to be approximately proportional to 
the time for which the cathode had been glowed. 
The same result was observed at 2230°K, 2130°K 
and 1940°K. We take this time as a measure of 
the amount of thoria deposited on the probe. 
Below the value of about ten times the normal 
emission the curves became irregular. Materials 
deposited on the ends of the filament may be the 
cause of this irregularity. With this measure of 
the amount of deposited thoria we can compute 
from data taken for different temperatures T of the 
cathode (curves similar to those shown in Fig. 15) 
the time, ¢., required for glowing the cathode to 
deposit a certain amount of thoria on the probe. 
Figure 16 shows the straight line relation between 
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Fic. 15. Time variation of probe emission at 1970°K (testing 


temperature 1390°K) after glowing the cathode at 2030°K for 
different intervals of time. 
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Fic. 16. Determination of the heat of evaporation of thoria. 


logiote and 1/7. The equation of this line is: 
log of. = const. — (39,400/T). (11) 


The heat of evaporation of thoria is thus 184,000 
calories per g-molecule. 


DIFFUSION OF THORIUM 
Equations (2) and (5) give: 
[d(logi)/dt]« DG—EN. 


By allowing the cathode to come to equilibrium 
at a fixed temperature 7, we established the same 
values of Gand N. The initial values, | d(logi) /dt| .~o, 
when the cathode is brought down from 7; to 
different temperatures should then be a combination 
of two exponential functions of 1/T determined by 
D and E according to (3). Figure 17 gives the curve 
log | d(logt)/dt| .-o versus 1/T for T; = 2040°K. These 
data are taken from curves similar to curves 2 of 
Figs. 6-8. For temperatures below 1700°K the curve 
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Fic. 17. Determination of the heat of diffusion of free thorium. 
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is straight down to 1100°K. At the lower tempera- 
ture the evaporation will be very slow. The slope 
of this line corresponds then to the heat of diffusion. 
The equation of this line is 


log 10| d(logz) /dt| :~0 = const. — (6840/T) 


giving the heat of diffusion to be 31,900 calories 
per g-atom. At 1700°K the curve begins to bend 
due to the increased predominance of evaporation. 


CONCLUDING REMARKS 


The addition of tungsten powder for the purpose 
of simplifying the heating problem seems to retain 
all advantages of the sintered thoria cathode, 
except for a somewhat higher work function. Thus 
at 2040°K the maximum space charge limited 
emission at the bend of the curve in Fig. 3 is 0.9 
amp./cm?, whereas Hanley’s data for thoria coated 
on tungsten give about 4 amp./cm?. This is due to 
the fact that part of the surface is occupied by 
tungsten which has a high work function. 

Our theory seems to give satisfactory explana- 
tions to the observed changes of activity with tem- 
perature. The probe measurements support our 
theory in confirming the evaporation of thorium. 
The estimated values of the heat of evaporation, 
46,000 cal./g-atom, is higher than the estimated 
value of heat of diffusion, 31,900 cal./g-atom, and 
is lower than the known value of the heat of 
reduction of thoria, 138,000 cal./g-atom.*® This is in 
accordance with the assumption of our theory. 

It appears from our results that all thoria cath- 
odes, where hot tungsten is in cdntact with thoria, 
should show changes of activity with temperature 
variation of the cathode. Indeed, Weinreich® ob- 
served that the emission of thoria coated on 
tungsten at 1800°K starting at a high value decayed 
with time, if the cathode had been previously 
flashed at a temperature above 2600°K. Thus in his 
case change of activity does occur, but only by 
flashing the cathode at temperatures above 2600°K. 
It is significant, however, that, similar to our case, 
the emission right after flashing was much higher 
than the equilibrium value at the lower tempera- 
ture, 1800°K. A thoriated tungsten filament would 
be completely deactivated at the flashing tempera- 
ture (above 2600°K) and would begin to activate at 
the lower temperature. As discussed above, this is, 
according to our theory, because the rate of thorium 
evaporation from thoria rises with temperature 
more slowly as compared with its evaporation from 
tungsten. 
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High Temperature Furnace for Electron Diffraction Studies of Thin Films 


EILEEN I. ALESSANDRINI 
Research Laboratory, General Electric Company, Schenectady, New York 


(Received January 18, 1949) 


A furnace has been designed which enables one to carry out electron diffraction studies of thin films at 
high temperatures. Electron diffraction patterns are obtainable at temperatures up to 1100°C. It is possible 
to bring the furnace up to 1500°C, but at temperatures over 1100°C the photographic plates are fogged by 
light emitted from the furnace. With the use of this furnace, recrystallizations and phase changes in thin 
metal and alloy films can be studied at the temperatures where these phenomena occur. 





GENERAL DESCRIPTION 


HIGH temperature furnace has been designed for 

use in making electron diffraction studies of thin 
metallic films of metals and alloys. It has been used 
with a laboratory model, but can be easily adapted to 
any electron diffraction instrument. The furnace is 
constructed so that it is possible to obtain transmission 
patterns of these films at temperatures up to 1100°C. 
The furnace can be mounted in place of the sample 
holder in an electron diffraction instrument, and vertical 
and horizontal motions of the thin film mount are pro- 
vided, so that any chosen portion of the sample can be 
brought into the electron beam. 

The furnace mounting consists of a vacuum seal 
through which it is possible to move the specimen, heat 
it and observe temperatures by means of a thermo- 
couple. The films are heated by a helically wound 
tungsten wire, which surrounds the outer edge of the 
circular film holder. A direct current is applied to the 
tungsten wire, and the heat is radiated to the film. 
A bank of four storage batteries is sufficient to obtain 
temperatures up to 1100°C. In order to obtain still 
higher temperatures, electron bombardment of the 
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Fic. 1. Furnace manipulator. 


VOLUME 20, JULY, 1949 


specimen can be accomplished. By connecting the 
positive terminal of one or more 90-volt batteries to the 
specimen holder and the negative terminal to one leg 
of the tungsten wire, the specimen can be bombarded 
with electrons and heated to temperatures over 1500°C. 
A radiation shield of tantalum is used to keep heat 
losses down to a minimum. Any amount of gases evolved 
when the furnace is run at high temperatures is readily 
taken care of by the pumps of the diffraction instru- 
ment. 

Specimens for use in the high temperature furnace are 
prepared by evaporation. A metal or metals can be 
evaporated onto a plastic film. This film of plastic plus 
metal, in turn, is supported by a fine meshed screen 
which is placed in the furnace. Before films are mounted 
in the furnace for high temperature studies, the plastic 
material is removed by a solvent, and therefore no 
contamination of the metal film occurs. 
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Fic. 2. High temperature furnace. 
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Fic. 3. Detail of furnace 
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SPECIFICATIONS 


The furnace mounting is arranged to pass through a 
vacuum seal, making it possible to move the specimen 
from outside. This arrangement is illustrated in Fig. 1. 
A stainless steel tube (1) has been made to fit through a 
heavy brass disk (2). The vacuum is maintained by 
means of a Wilson seal (3). The steel tube is used to give 
support to the specimen holder and also to make it 
possible to move the specimen vertically and _ hori- 
zontally in the electron beam. Additional support is 
given the tube by a }-inch track (4). A slide and knurled 
nut (5) are used to secure horizontal positioning of the 
specimen, and a split chuck and knurled nut (6) are used 
to hold the specimen in any desired vertical position. 
The two filament terminals and thermocouple leads are 


brought out through the heavy brass disk (not shown in 
drawing). Glass beaded seals are used to insulate them. 

On the end of the stainless steel tube which extends 
into the vacuum (Fig. 2), a stainless steel plug (1) is 
copper brazed. A slit and screw are provided which 
hold a 20 mil, $-inch wide molybdenum lath (2). This 
ribbon lath terminates in a disk 4 inch in diameter with 
a $-inch aperture in its center. This disk has a lip which 
is 3 inch in width, as shown in Fig. 3. 

A 10 mil helically wound tungesten wire (3) encircles 
the lip of this disk without touching it. This wire serves 
as the heating element of the furnace. This heating 
element is supported by two 80 mil molybdenum wires 
(4). Copper braided wires ‘extend from these supports 
to the terminals (8). The thermocouple wires (5), of 10 
mil Chromel and Alumel, are spot welded to the molyb- 
denum disk at the edge of the aperture. The terminals 
(9) for these wires extend to the outside of the furnace 
through glass beaded seals. All of these leads are held 
rigidly in place along the molybdenum lath by means 
of a mica lined stainless steel block (6), which consists 
of two pieces screwed down tightly against the lath. 

In order to mount specimens in this furnace, a molyb- 
denum disk $$ inch in diameter and having an aperture 
$ inch in diameter was made (Fig. 4). This disk is 20 
mils in thickness and has a } inch ribbon extending 1} 
inches from an edge (Fig. 4 (1)). At the edge of the ribbon 
a piece of 5 mil molybdenum was spot welded to form 
tabs (Fig. 4 (2)) for securing the disk to the furnace. 
This specimen holder fits flush against the molybdenum 
disk of the furnace, and the two }-inch apertures are 
concentric. Across the aperture of the specimen holder 
a 150-mesh platinum-rhodium screen (Fig. 4 (3)) was 
spot welded. This screen serves as the film support. 

A radiation shield of tantalum is used to cut down 
heat loss to the instrument. This shield is a piece of 
tantalum which is spot welded to form a cylinder 
(Fig. 2(7)) ? inch in diameter. A disk of tantalum is spot 





Fic. 5. Germanium evaporated on a high temperature 
furnace mesh. 
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Fic. 6. Germanium heated to 350°C showing 
recrystallization. 
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welded across the bottom and three tabs were cut out 
at the top. This shield fits around the furnace and is 
secured by means of the tabs to the stainless steel, 
mica lined, supporting block (Fig. 2(6)). Two 3:-inch 
apertures were made for the entrance and exit of the 
electron beam. 

A film for diffraction studies is mounted in the fur- 
nace, the filament of which is heated with the use of 
storage batteries. Twelve cells are connected in series, 
and the amount of current through the filament is 
controlled by a variable resistor. In this way, specimens 
can be brought to temperature gradually. For tempera- 
tures over 1000°C, the specimens are bombarded with 
electrons. Dry batteries are connected across the speci- 
men holder and one filament terminal. If the specimen 
is made positive, electrons from the filament bombard 
the molybdenum disk and cause an increase in tem- 
perature. 
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RESULTS 


Several preliminary experiments have been carried 
out using this high temperature furnace. Germanium 
was evaporated on Formvar. The Formvar was removed 
by placing the film in a solution of ethylene dichloride. 
The metal film was then mounted in the high tempera- 
ture furnace and a diffraction pattern was taken (Fig. 5). 
This pattern shows four diffraction rings. The germa- 
nium was heated to 350°C, and recrystallization oc- 
curred, as shown by the increased number and sharpness 
of the rings (Fig. 6). 
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A Method of Equivalent Linearization for Non-Linear Oscillatory 
Systems with Large Non-Linearity 


Cari. A. LUDEKE 
University of Cincinnati, Cincinnati, Ohio 
(Received October 15, 1948) 


This paper introduces a method, based partially on existing non-linear theory and partially on ex- 
perimental results, for transforming a non-linear oscillatory system into an equivalent linear one, 
regardless of the amount of non-linearity in the original system. 


INTRODUCTION 


HE applied physicist, concerned with non- 

linear forced damped oscillations, is interested 
in the wave form, the frequency-amplitude relation- 
ship, and the dependence of the maximum ampli- 
tude on damping. In most practical problems, the 
damping force and the applied force are so small 
compared to the restoring force that if questions 
concerning these three points of interest can be 
satisfactorily answered, the complete behavior of 
the system is essentially established. Although 
methods are available for answering such questions 
if the non-linearity is small, this author believes that 
more should be said about the application of these 
methods to systems with large non-linearity. 

Of the methods available two have been chosen as 
characteristic of all others, and the difficulties in- 
volved in applying them to systems with large non- 
linearity are cited. In an attempt to remove some of 
these difficulties, the author uses experimental re- 
sults obtained by means of an apparatus previously 
reported.! 


THEORY 


Before examining these two methods, let us con- 
sider the differential equation, 


16+k,0+k;6=0, (1) 
and assume it has a solution 
6=A sinwt+ pA sin3ot, (2) 


in which p is so small that terms of the order of p? 
may be neglected. Under what conditions is this 
assumption approximately correct? 

Upon substituting Eq. (2) into Eq. (1) and neg- 
lecting terms of the order of p? we have 
[ (ki — Iw?) A +3k3A*/4 | sinwt 
~ +[ (ki —9Iw*) pA + 3pk3A*/2 —k3A*/4] sin3wt 

— (3pk3A*/2) cos2wt sin3wt=0. (3) 

Obviously Eq. (3) can be true for all values of ¢ only 
if p=0, in which case k; =0, and w?=k,/J; this is the 
well known linear case. However, with p#0, we can 
still find an expression for w? which will satisfy 


1C. A. Ludeke, J. App. Phys. 20, 600 (1949). 
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Eq. (3) for a particular value of ¢. When restricted 
to one value of ¢, the usual choice is t=27/w, in 
which case we have the relationship, 


Tu? =[(1—p)/(1—9) Jer 
+[(1—3p)/(1—9p)(1—p) Than? (4) 


in which @,, is the maximum value of @. 

If the solution represented by Eq. (2) exists with 
small values of p regardless of the magnitude of the 
terms in Eq. (1), then it is possible to write Eq. (4) 
in the approximate form, 


Ta? —~k 1 +h 3On?—~k1 + k3A?. (5) 


If, however, a small value of » demands a small 
value of ks, then Eq. (3) reduces to 


[ —AwI+k,A + 3k3A'/4 | sinwt =0, (6) 


which can be satisfied for all values of ¢, by the 
expression 
Tw? = ki +3k3A?/4. (7) 


A consideration of Eqs. (5) and (7) indicates (a) 
that the form of the expression relating frequency 
and amplitude might remain the same for all values 
of non-linearity, and (b) if the solution remains 
essentially sinusoidal for large non-linearities, then 
the coefficient of the term 3A? is 3/4 when the non- 
linearity is small and approximately 1 when the 
non-linearity is large. 

If the right side of Eq. (1) is replaced by Pp sinat, 
Eq. (1) becomes the so-called Duffing? equation, and 
Eqs. (5) and (7) are replaced by 


Te? ~ki +h 38m? — Po/Om~kitk3:A?—Po/A, (8) 


and 
Iw? =k, +3k3A?/4—Po/A (9) 


respectively. 

in order to determine whether or not the solutions 
are approximately sinusoidal when the non-linearity 
is large, and in order to determine the best coeffi- 
cient for the term k3A? let us consider a well-known 
method? available for solving such problems when 
the non-linearity is small. In reference 3, p. 135, this 


2G. Duffing, “Sammlung Vieweg,”’ Heft 41/42, 1918. : 
3N. Minorsky, Introduction to Non-Linear Mechanics (J. W. 
Edwards, Ann Arbor, Michigan, 1947). 
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Fic. 1. Frequency-amplitude data for a typical non-linear 
free oscillation. 


restriction of small non-linearity is explicitly stated. 
In article 71, p. 222, Eq. (1) is considered, and from 
Eqs. (71.5) and (71.6) on p. 223 it is evident that 
these equations which give the solution and the 
frequency-amplitude relationship are only true 
when the non-linearity is small. Thus the applied 
physicist would not find the method discussed in 
reference 3 suitable for problems with large non- 
linearity even though it is a method capable of 
yielding approximations of higher orders. 
As a specific example, consider Eq. (1) written in 
the form 
d*6/d7r?+ 6+ »é=0, (10) 


in which 7 = (k,/I)'t =wot and p=k3/Iwo?=k3/k1 


TABLE I. Experimentally determined values of C, the coeffi- 


cient of the term 8@,,2, for systems having a restoring torque 
ki0+k36°. 











Spring 

identi- ki ks 

fication c (g cm/rad.) (g cm/rad.*) 
13;0-0 1.01 1.27107 15.2 «10° 
13;0-1 1.01 1.17 14.0 
13;0-2 0.93 1.14 14.0 
13;0-3 0.94 1.10 11.9 
13;0-4 0.89 1.06 11.9 
13;0-5 0.84 1.03 51.3 
13;0-6 0.99 0.98 9.8 
13;0-7 0.95 0.94 10.3 
13;0-8 0.95 0.92 9.7 
20;0-3 0.88 0.11 9.0 
20;0-4 0.90 0.02 7.7 
20;0-2 0.94 0.00 8.1 
20;0-8.A 0.83 0.04 6.9 
20;0-1 0.93 0.07 9.3 
20;0-8B 0.87 0.49 7.0 
20;0-0 0.86 0.10 10.2 
20;0-2A 0.82 0.15 9.7 
20;0-2B 0.86 0.68 8.6 
21;0-0 0.87 1.01 16.3 
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Fic. 2. Frequency-amplitude data for several typical non- 
linear forces oscillations. 


= 9200. With this large value of u, the higher power 
terms in Eq. (71.6) of reference 3 would predomi- 
nate and hence this equation could not be used. The 
actual experimental results are shown in Fig. 1, 
from which we see that it is only necessary to retain 
the term in the square of the amplitude. Thus, as 
clearly stated in reference 3 this method should not 
be used for systems with large non-linearity. 

Consider now the work of Friedrichs and Stoker,* 
p. 36, in which we encounter an iteration method 
which does not depend explicitly upon the magni- 
tude of the non-linearity but which does start with 
the free linear oscillation as the initial solution. 
Indication that the magnitude of the non-linearity 
is involved is given by the author of reference 4 on 
p. 39 where he states that the frequency-amplitude 
relationship can be expected to be accurate for not 
too large values of the term containing the square of 
the amplitude. Since the method under discussion 
can be applied repeatedly, let us consider the results 
of two such applications on the Duffing equation 
when it is written in the form, 


6+ W070 + BO = F coswt, (11) 
in which wo? =k;/I, B=k3/I and F=P,/I. 


The first application yields as a solution for 
Eq. (11) the expression 


6,=A coswt 
+ [BA*/32(wo?+38A?/4— F/A)]cos3wt. (12) 


The first approximation to the frequency-amplitude 
relationship is 


wi” = wo’ + 38A?/4—F/A. (13) 


*K. O. Friedrichs and J. J. Stoker, “‘Introduction to Non- 
Linear Mechanics,”” New York University Lecture Notes. 
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TABLE II. Results of the harmonic analysis of three typical! 
wave forms. (Note that the third harmonic is the only higher 
harmonic present to any appreciable amount, and ‘that its 
emgutade € decreases as the non- rlimearity ¢ decreases. i 











@=a0+ SE an cosnut + SB bn sinnut 
a7! ant 








Non- Equation 
Spring constants n an bn linearity (approximately) 
1 0.00 —2.34 8.67 
k, =0.13 107 2 0.00 0.02 6= —2.34 sinwl 
ks=11.95X10® 3 0.00 0.12 +0.12 sin3wt 
4 —0.01 0.00 
5 0.01 0.00 
6 0.01 0.01 
7 0.02 0.01 
8 0.01 0.01 
9 0.01 0.01 
10 0.01 0.01 
1 —0.11 2.79 1.97 
k, = 0.60107 2 —0.05 0.00 6=2.79 sinwl 
ks=11.95X10® 3 0.01 —0.07 —0.07 sin3wt 
4 0.01 01 
5 0.00 02 
6 0.01 01 
7 0.03 —0.01 
8 0.00 01 
9 0.00 00 
10 0.00 01 
0.06 0.98 
ky, =1.48X 10? —0. 6 = 3.98 sinwt 


k3=11.95 X 10° 
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* Note: (1) an and b» are measured in inches with the wave-length equal 
to 6 inches. (2) The smallest reading on the analyzer was 0.01 inch. (3) 
From theoretical considerations there should be no.even harmonics. 


A second application of the process yields as a 
solution 


62= A coswt — (1/8w*) {[ (wo?) /32u1?— 1/4 8A? 
— 38°A*/64w,? — 3BA®/131072} cos3wt 
+ (1/802) [BA ®/128w32-+8°A7/4096w3"] cosSwt 
+ (1/16w?)[8*A7/4096w;* ] cos7wt 


— (1/320w?) [B*A*®/32768w1°] cos9wt (14) 


while the frequency-amplitude equation becomes, 
wo” = wo?+ (38A?/4)[1+8A?/32w:t]—F/A. (15) 


Let us examine Eqs. (12), (13), (14), and (15) 
when the non-linearity is so large that Eq. (13) can 
be written 


wyr~3BA?/4, (16) 
Equation (12) then becomes 
0:~A coswt+ (A /24) cos3wt, (17) 


which consists predominantly of the term A coswt, 
and which does not contain 8 explicitly. Similarly, 
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Eqs. (14) and (15) become, 


~A coswt+ (3/64)A cos3wt 
+(1/600)A cosSwt+ (1/25248)A cos7wt 
+ (4.43X10-*)A cos9wt, (18) 


_and 


= wo?+ (3/4) (397/384)BA? 
— F/A~(3/4)(397/384)BA?. (19) 


From Eqs. (18) and (19) we see that even for large 
non-linearities the solution is still dominated by the 
term, A coswt, while the frequency-amplitude re- 
lationship is altered by a slight increase in the 
coefficient of the term 8A*. The discussion leading to 
Eq. (8) indicates that the eet limit of this coeffi- 
cient is one. 

If only the first two terms of Eq. (18) are retained, 
and if the maximum value of @ is designated with a 
subscript m, then Eq. (19) can be written 


= wo? +0.71802m? — F/ 02m. (20) 
Under similar circumstances Eq. (13) can be written 
@i 2 = w,_?+0. 69861m? — F/@im. (21) 


By comparing Eqs. (20) and (21) we see that as the 
iteration process is repeated the coefficient of the 
term 86,,? does increase slightly. The applied physi- 
cist is therefore faced with the task of repeating the 
process until this coefficient is determined accu- 
rately enough to give satisfactory results for the 
particular non-linearity involved. It seemed to this 
author much more appealing to determine the 
coefficient of BA? from experimental results. 

The results of many experiments indicates that 
for systems having non-linearities between 0.1 and 
o, the value of 0.9 is a good approximation for the 
coefficient of the term §6,,?. For specific results see 
Table I. In Fig. 2 we see the actual experimental 
values of amplitude and frequency obtained for 
three typical non-linear systems with different 
amounts of non-linearity. In these three examples 
the term F/6,, was negligible. In Table II we see the 


TABLE III. Typical experimental results showing that the 
relationship between the amplitude @,, and the frequency », 
is essentially independent of the small damping torque 








4m in millirad. ~ 6.4 12.9 19.3 25.8 32.1 38.6 45.1 





Constants c * w in radians per second 

1.05 10.2 12.9 16.4 20.9 25.3 29.8 33.8 

2.39 10.2 13.2 16.9 21.2 25.8 30.2 34.2 

ki =0.15 X10? 2.64 10.5 13.3 16.9 21.3 25.8 29.8 33.9 
k3 =9.7 X10° 3.04 10.7 13.3 17.3 21.3 26.8 30.2 34.2 
3.46 10.7 13.8 17.3 21.3 26.0 30.2 34.2 

k =32.6 

ki =0.10 X10? 1.05 12.4 16.9 21.8 26.3 31. 35.3 
ks =10.2 X10° 1.95 9.8 12.9 16.9 21.8 26.5 31.1 35.4 
k =32.6 2.39 9.8 12.9 16.9 21.8 26.7 31.1 35.6 
2.64 9.8 12.9 17.6 21.8 26.4 31.1 35.6 

3.04 9.8 13.2 17.4 22.2 26.7 31.1 35.6 

3.46 10.2 13.3 17.8 22.2 26.9 31. 35.7 








* The magnitude of the applied torque is kw. 
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TABLE IV. Typical phase angle data for non-linear systems.* 








Amplitude Frequency Phase angle Non- 





Spring constants millirad. rad./sec. degrees linearity 
3.2 24.9 17 0.01 
k, =1.48 X10? 6.4 27.9 20 0.03 
k3=11.08 X 10° 14.2 29.2 27 0.15 
19.3 33.1 34 0.28 
25.8 36.3 35 0.50 
32.8 39.2 38 0.80 
45.0 44.8 52 1.51 
x R x R 90 
3.9 15.1 35 0.03 
k, =0.60 X 107 7.1 17.8 50 0.09 
k;= 11.08 x 10° 12.9 20.1 49 0.30 
19.7 23.5 40 0.71 
26.7 27.6 46 1.38 
32.2 30.8 50 1.92 
35.4 32.9 53 2.31 
x R x R 
20.3 17.8 78 3.51 
k, =0.13 X10? 26.1 24.5 70 5.80 
k;=11.08 X 10° 30.0 25.1 72 7.67 
33.2 27.0 69 9.40 
35.4 29.0 69 10.70 
38.0 $8.2 72 12.31 
41.8 33.5 82 14.88 
43.8 34.5 82 16.35 
45.1R 35.2R 90 17.30 








* Note: (1) R indicates resonance. (2) % indicates off scale. 


results of the harmonic analysis of the wave forms of 
several typical solutions. We nete that even for 
large non-linearity, the fundamental is the domi- 
nating term in the solution. 

Consider next the effect of introducing a small 
amount of viscous damping into the system under 
consideration. An examination of reference 3 shows 
that the methods therein considered are suitable for 
systems with small non-linearities. For small 
damping the methods yield a frequency-amplitude 
relationship essentially the same as Eq. (9) except at 
resonance. At resonance the methods tell us that the 
phase angle must be approximately 90° and the fre- 
quency and amplitude must, in addition to Eq. (9), 
satisfy the relationship 


A,=F/cw, (22) 


in which A, is the resonance amplitude, F is the 
magnitude of the applied force, c is the coefficient 
(small) of viscous damping, and w is the frequency. 
We note that Eq. (22) is the usual approximation 
for the linear case when the damping is small; how- 
ever, w is now a function of A. See Eq. (100.30), 
reference 3. 

Turning now to reference 4, p. 42, we find that 
although the procedure is different the first ap- 
proximation yields the same result. For example, 
consider a system represented by the differential 
equation, 


6+ 6+? +86' =F coswt—G sinwt. (23) 
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Upon assuming an initial solution, 6)=A coswt, and 
imposing the conditions of periodicity, we find 


and 


cAw=G, 


(wo? —w*)A +38A*/4= F. 


(24) 


(25) 


Equations (24) and (25) can be combined in the 


form 


S?(w, A)+CA*%w? = H?, 


(26) 


TaBLE V. A comparison of the experimental resonance 
amplitude @,, and the theoretical amplitude F/cw, for various 


non-linearities. 








Spring 





F/cw Experimental Non- 
identification millirad. 6, (millirad.) L =¢w6,/F linearity 
13;0-0 20.0 28.2 1.41 0.95 

15.8 20.7 1.31 0.51 

12.5 14.2 1.14 0.24 

10.6 11.5 1.09 0.16 
13;0-1 18.1 24.8 1.37 0.74 

14.6 18.3 1.25 0.40 

11.8 13.3 1.13 0.21 

10.1 - 10.7 1.06 0.14 
13;0-2 17.1 23.1 1.35 0.66 

14.0 17.0 1.21 0.36 

11.5 12.7 1.10 0.20 

9.9 10.4 1.05 0.13 
13;0-3 16.1 21.7 1.35 0.51 
13.5 16.0 1.19 0.28 
11.2 12.4 1.11 0.17 
9.7 10.0 1.03 0.11 
13;0-4 15.5 20.6 1.33 0.48 
13.1 15.3 1.17 0.26 
11.0 12.0 1.09 0.16 
9.5 9.8 1.03 0.11 
13;0-5 14.6 19.3 1.32 0.40 
12.7 14.6 1.15 0.23 
10.7 11.7 1.09 0.15 
9.3 9.6 1.03 0.10 
13;0-6 14.4 17.7 1.23 0.31 
12.3 13.8 1.12 0.19 
10.4 11.3 1.09 0.13 
9.0 9.3 1.03 0.09 
13;0-7 13.9 16.9 1.22 0.31 
12.1 13.4 1.11 0.19 
10.2 11.0 1.08 0.13 
8.8 9.1 1.03 0.09 
13;0-8 13.6 16.1 1.18 0.27 
11.8 12.9 1.09 0.17 
10.1 10.8 1.07 0.12 
8.8 9.0 1.02 0.08 
20;0-2B94 23.7 45.1 1.90 13.4 

21.6 38.5 1.78 9.8 

19.0 he Ie 1.70 6.8 

16.4 25.8 1.57 4.4 

15.1 19.3 1.28 2.5 
20;0-2B93 23.0 36.4 1.58 2.1 

19.7 30.3 1.54 1.5 

16.3 24.2 1.48 0.9 

13.6 17.7 1.30 0.5 

9.3 10.3 1.11 0.2 
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Fic. 3. The factor L as a function of the non-linearity. 


in which 


Si(w, A) = (wo? —w?)A +3/48A?, (27) 
and 
I? = F?+G?*. (28) 


At resonance S;(w, A) =c?A/2, H=F, and for small 
values of c, Eq. (26) becomes 


A,=F/cw, (29) 


which is the same as Eq. (22). It is perhaps well to 
mention at this time, that in non-linear systems the 
resonance curve is discontinuous and exhibits the so- 
called ‘“‘jump” phenomena. These discontinuities or 
jumps occur at different frequencies depending upon 
whether the frequency is being increased or de- 
creased. In the system described by Eq. (23) with 
B-positive, the larger resonance amplitude occurs 
when the frequency is being increased. Throughout 
this paper this is the resonance amplitude referred 
to. 
If we repeat the process assuming a solution 


6,:=A coswt+ (BA*/32w,") coswt, (30) 
Eq. (24) remains unchanged, and Eq. (25) becomes 


(wo? — w*) A + (38.4*/4) 
X [1+BA?/320+6°A'/5120:4]=F (31) 


in which w,;’ is given by Eq. (13). If we consider now 
systems with non-linearity so large that Eq. (16) can 
replace Eq. (13) then Eq. (31) becomes 

(wo? — w*)A + (38/4)(397/384)A*=F. (32) 


Upon squaring and adding Eqs. (24) and (32) we 
have 
‘ S2?(w, A)+CA*w’ = H’, (33) 
in which 

S2(w, A) = (wo? — w*) A + (38/4) (397/384)A*. (34) 


Proceeding as before, S2(w, A) =c?A/2 at resonance, 
and for small values of c, Eq. (33) becomes, 


A,=F/cw. (35) 
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A comparison of Eqs. (29) and (35) shows that no 
matter how many times we repeat this process, the 
amplitude of the fundamental at resonance will 
always be given by the same equation utilizing how- 
ever higher and higher approximations for the fre- 
quency. In terms of resonance amplitude of the 
solution rather than of the fundamental, Eqs. (29) 
and (35) become 


61,= F/cw, (36) 


TaBLe VI. A comparison of the ratio w/wo with L, the 
corrective factor in the equation 0,=LF/cw. 








Spring Non- 
identification w/wo L linearity 
13;0-0 1.40 1.41 0.95 

1.22 1.31 0.51 
1.11 1.14 0.24 
1.07 
13;0-1 1.32 1.37 0.74 
1.18 1.25 0.40 
1.10 1.13 0.21 
1.07 1.06 0.14 
13;0—2 1.27 1.35 0.66 
1.15 1.21 0.36 
1.08 1.10 0.20 
1.06 1.05 0.13 
13;0-3 Be 1.35 0.51 
1.13 1.19 0.28 
1.08 1.11 0.17 
1.06 1.03 0.11 
13;0-4 1.20 1.33 0.48 
1.11 1.17 0.26 
1.07 1.09 0.16 
1.05 1.03 0.11 
13;0-5 1.14 1.32 0.40 
-1.09 ‘3 0.23 
1.06 1.09 0.15 
1.04 1.03 0.10 
13;0-6 1.15 1.23 0.31 
1.09 1.12 0.19 
1.06 1.09 0.13 
1.04 1.03 0.09 
13;0-7 1.14 1.22 0.31 
1.09 1.11 0.19 
1.06 1.08 0.13 
1.04 1.03 0.09 
13;0-8 1.12 1.18 0.27 
1.08 1.09 0.17 
1.06 1.07 0.12 
1.04 1.02 0.08 
20;0-2 B94 3.47 1.90 13.4 
3.07 1.78 9.8 
2.66 1.70 6.8 
2.20 1.57 4.4 
1.76 1.28 2.5 
1.37 1.09 1.2 
20;0-2B93 1.86 1.58 23 
1.67 1.54 2.5 
1.48 1.48 0.9 
1.32 1.30 0.5 
1.14 a.43 0.2 
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and 
62, = (25/24) F/cw. (37) 


A comparison of Eqs. (36) and (37) reveals that for 
systems with large non-linearities the amplitude at 
resonance will be larger than that predicted by the 
theory of the first approximation. 

Thus for large non-linearity we have the following 
conclusions concerning periodic solutions of Eq. (23) 
which have a frequency the same as the frequency 
of the forcing function: 

(1) For most practical purposes, the solution can 
be considered sinusoidal in form. 

(2) The frequency-amplitude relationship for all 
large non-linearities is approximately 


w? = wo? +0.986,,? — F/Om, (38) 


and is essentially independent of damping except at 
resonance. 

(3) At resonance not only must Eq. (38) be 
satisfied but also the relationship 


6,=LH/cw, L>1 (39) 


must hold true. In Eq. (39) the value of L depends 
upon the value of the non-linearity; Z approaches 1 
as the non-linearity approaches zero. 

As with the coefficient of the term £8,,”, this 
author decided to determine the relationship be- 
tween Z and the non-linearity from experimental 
results. 


EXPERIMENTAL RESULTS 


Experimental results have already been used to 
verify our first two conclusions; and some additional 
data is displayed in Table II] to show that the 
frequency-amplitude relationship is essentially inde- 
pendent of small amounts of viscous damping except 
at resonance. 

For the case of resonance we note from Table IV 
that, as expected, the phase angle between the 
applied force and the solution is approximately 90° 
at resonance if the damping is small. The depend- 
ence of the factor Z on the non-linearity k36,7/k, is 
shown in Table V. We note that L>1, L—1 as the 
non-linearity approaches zero, and L increases as 
the non-linearity increases. These are just the re- 
sults which can be inferred from Eqs. (36) and (37). 
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A plot of Z versus non-linearity is shown in Fig. 3, 
and supplies the information needed to make Eq. 
(39) valuable in the actual solution of problems. 
It is interesting to note at this point that Eq. 
(100.23) p. 316 of reference 3 indicates that for not 
too large non-linearity the value of L should be 
w/wo if the forcing function is supplied by a rotating 
unbalance. A comparison of L and w/w is made in 
Table VI and we note that they are approximately 


equal provided the non-linearity remains less than 
one. 


EQUIVALENT LINEARIZATION 


Returning to Fig. 3, we note that the results of 
this figure together with Eqs. (38) and (39) are 
sufficient to determine 6,. The procedure is simply 
to choose a value of 6,, compute the non-linearity 
k30,2/k,, and with this value of non-linearity find the 
value of Z from Fig. 3. Using this value of Z and the 
chosen value of 6,, solve Eq. (39) for w. This value of 
w and the original choice of @, must then satisfy 
Eq. (38). When this is true, the choice of 6, is correct. 

It is obvious that this procedure, which resulted 
from a study of Eq. (23) can be generalized to 
include systems with small viscous damping and 
any kind of a non-linear restoring force. It should be 
noted, however, that the periodic solutions are re- 
stricted to have the same frequency as the forcing 
function. These are not the only periodic solutions 
to equations of this type and the author will report 
on his investigations on subharmonic solutions in a 
subsequent paper. 

Since for most practical purposes, the solutions 
examined in this paper can be considered sinusoidal, 
and since Eq. (39) differs from the resonance ampli- 
tude equation for a linear system only by the factor 
L, the procedure outlined above is essentially a 
method of equivalent linearization for systems with 
large non-linearity; the equivalent linearization of w 
being defined by Eq. (38). 
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The Illuminating System of the Electron Microscope 


James HILuier Anp S. G. ELtis 
RCA Laboratories, Princeton, New Jersey 


(Received December 2, 1948) 


Using a self-biased electron gun the current arriving at the specimen, the current density at the 
specimen, and the area of illumination of the specimen have been determined for different values of 
the filament height. The variation of current density at the specimen with angular aperture of illumi- 
nation has been determined and the results compared with those for the zero-biased gun. The per- 
formance of the self-biased gun is assessed and some of the factors involved in its operation and 


maintenance are discussed. 





INTRODUCTION 


HE purpose of this paper is to provide data 

on the illuminating system of an electron 
microscope when it is operated with an electron 
gun biased with a cathode resistor and operated at 
saturated beam current. For simplicity in termi- 
nology such an electron gun operated at saturated 
beam current will be called a self-biased gun. The 
paper also contains a discussion of some of the 
factors to be taken into consideration in the opera- 
tion of such an illuminating system. 

Hillier and Baker! have described the adjustment 
and performance of an illuminating system em- 
ploying a fixed-bias gun with zero bias. Ellis? has 
described the optics of three-electrode electron guns 
and has also described qualitatively the behavior of 
the self-biased gun. References to the earlier litera- 
ture will be found in the above papers. 

For efficient operation the illuminating system 
should satisfy the following requirements. First, it 
should illuminate only that part of the specimen 
which is under examination, for illumination which 
reaches other parts of the specimen serves only to 
heat the specimen. The temperature of the region 
under examination and therefore the tendency to 
specimen destruction is increased without any gain 
of image intensity.* * The general heating also tends 
to cause thermal drifting of the specimen. The 
illumination of parts of the specimen that are not 
under observation has the further undesirable effects 
that these parts of the specimen are contaminated, 
and also that scattered electrons from them may 
enter the objective lens and decrease the contrast 
in the final image. The electron gun should there- 
fore provide a small and preferably circular source. 
Secondly, the illuminating system should provide 
sufficient intensity of illumination at the desired 
angular aperture that the microscope may be 
readily focused. Thirdly, the angular aperture of 
the illumination should, for most purposes, be as 
small as possible, consistent with the foregoing 

1 J. Hillier and R. F. Baker, J p. Phys. 16, 469 (1945). 

25. G. Ellis, J. App. Phys. a a5 9 (1947). 


3M. von Ardenne, Kolloid Zeits. 108, 195 (1944). 
‘V. K. Zworykin and J. Hillier, J. App. Phys. 14, 658 (1943). 
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requirements. Small angular aperture of the illumi- 
nation is desirable for two reasons: if the objective 
lens has any asymmetries it is better to use the 
more paraxial portions of the lens. This minimizes 
the effects of the asymmetries which become more 
apparent as the rays traversing the objective do so 
further from the axis. The other reason for requiring 
small angular aperture of illumination is that it will 
tend to sharpen and make more distinct the con- 
tour and phase phenomena in the final image and 
hence facilitate the detection of exact focus. It is 
now possible and indeed desirable to sacrifice in- 
tensity in order to achieve small angular aperture 
of illumination since the focusing can be performed 
at lower intensity of illumination if a high magnifi- 
cation viewer is employed either with or without 
a deflection focusing device of the type described 
by le Poole® and modified by the authors. The total 
beam current reaching the specimen should be kept 
as small as possible, consistent with the foregoing 
requirements, since for all the condenser strengths 
employed in practice this beam current, apart from 
the fraction hitting the specimen supporting screen, 
passes into the objective lens. If an objective dia- 
phragm is employed, that part of the beam current 
which does not go through the diaphragm deposits 
contamination upon it. When the beam current is 
large the rate of contamination is large and more- 
over the static charge maintained by the beam 
current on-the diaphragm is considerable so that 
the fraction of the beam current which proceeds 
through the diaphragm to form the intermediate 
image suffers irregular deflections as a result of the 
electrostatic field produced by this charge. Even if 
no objective diaphragm is used, the rate of con- 
tamination of the objective lens will increase as 
the current reaching the specimen is made larger. 
The foregoing optical requirements cannot be met 
independently and are, in fact, to a certain extent 
contradictory. In practice a compromise is reached 
by producing a small single intense source in the 
electron gun and by choosing the best condenser 


5 J. B. le Poole, Philips Tech. Rev. 9, 33 (1947). 
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diaphragm radius and condenser lens strength for 
the immediate purpose in view. 

In addition to the optical requirements, certain 
mechanical requirements have to be met. These are 
briefly that the illuminating system shall be as 
simple in construction as possible, that it shall not 
be too difficult to align it to the objective lens, and 
that the illuminating system shall be readily de- 
mountable for the renewal of the filament and the 
cleaning of the whole system. Finally, the operator 
may in special cases, wish to employ unusually high 
current densities of bombardment of the specimen 
or other unusual conditions of operation. It is, 
therefore, convenient if possible to avoid making 
the design so rigid in concept as to preclude such 
special departures from normal practice. 

The measurements reported below were made on 
the RCA model EMU electron microscope in use in 
this Laboratory. The symbols employed in this 
paper are those previously introduced in references 1 
and 2. They are listed below, together with the 
parameters of the illuminating system of this 
microscope. 


uz=spacing of filament above lower surface of the gun 
focusing aperture. 


ue=height of effective source above the condenser lens 
center. 


v-=distance from condenser lens center to the specimen 
= 150 mm. 
2r,'’ =diameter of focusing aperture = 3.75 mm (0.15 in.). 
2r.'' =diameter of condenser aperture =0.5 mm (0.020 in. ). 
2r,’=diameter of region in specimen plane illuminated when 
the intensity is a maximum, (i.e., focused illumi- 
nation). 
f-=focal length of condenser lens. 
I.=condenser lens current (in test meter scale divisions 
unless otherwise stated). 
Ig=beam current (in beam current meter reading minus 
leakage current shown when the filament is off). 
ig=current reaching the specimen plane. 
o=current density in beam reaching the specimen plane. 
oo=maximum value of o. 


n’ =o0/Ig the illuminating efficiency of the illuminating 
system. 


a«-=aperture of the illumination at the specimen. 


The shield to anode potential difference was 50 
kilovolts. The filament diameter employed was 4 
mils (=10-? cm), the focusing aperture to anode 
spacing was 15 mm (approx.), and the anode 


aperture 18 mm. The biasing resistor was 0.5 
megohm. 


EXPERIMENTAL PROCEDURE 


With the condenser lens strength adjusted to 
provide maximum illumination, the following quan- 
tities were measured: the total current arriving at 
the specimen, the distribution of current density in 
the illuminating beam at the specimen, and the 
size and shape of the area of illumination at the 
specimen. These measurements were obtained for 
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different values of the filament-aperture spacing 
using the self-biased gun. For comparison purposes 
similar measurements were made on the zero bias 
gun by shorting out the cathode biasing resistor and 
adjusting the filament temperature to produce the 
same total beam current leaving the filament. 

The apparatus used in making these measure- 
ments may be described with the aid of Fig. 1. 
A collector plate E about 1 square centimeter in 
extent was mounted on the axis of the microscope, 
midway between the projector lens and the final 
viewing screen. It was maintained at a potential of 
one hundred volts positive with respect to the body 
of the instrument by means of the battery B and 
the current reaching it was measured by a sensitive 


‘electronic microammeter M. Parallel plates D,, D2 


were mounted between the objective and condenser 
lenses so that the illuminating beam could be de- 
flected through small angles. 

The total beam current was measured by re- 
moving the projector polepieces and focusing the 
microscope so that the current passing through the 
specimen plane was completely collected by the 
plate E. In measuring the distribution of current 
density at the specimen plane the microscope was 
set at normal focus and operated at a magnification 
of approximately 6000. The microscope was aligned 
so that the microammeter gave a maximum reading 
and then voltages were applied to the plates D; 
and D: so that the collector E received, successively, 
currents from areas spaced out radially from the 
center of the illuminated region. The filament was 
oriented in the electron gun so that the deflection 
produced by the plates D;, Dz was parallel to the 
short diameter of the area of illumination. The col- 
lector E was conjugate to an area of approximately 
3X3 microns in the normal specimen plane. The 
measurements of distribution of current density 
were, therefore, somewhat crude. To obtain the size 
and shape of the area of illumination of the speci- 
men, the microscope magnification was changed to 
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Fic. 2. Variation of total beam current with filament height. 
Full line experimental self-biased gun used in these experiments 
(rx’’ =0.150 inch). Dashed line standard RCA self-biased gun 
for comparison (r;’’ =0.125 inch). 


approximately 1000 and the spot of illumination 
photographed with exposures of 1, 2, 4, 8, and 16 
seconds. Traces from these plates were obtained on 
a Leeds and Northrup recording microphotometer 
and from these traces, it was possible to determine 
the distribution in intensity and also the half-width 
of the area of illumination along the long and short 
diameters. In order to be able to determine the in- 
tensities of illumination under the conditions in 
which photography is normally carried out, the 
relation between condenser lens current and lens 
strength was determined in the following way. The 
radius ro’ of the area of illumination at the specimen 
was measured for different condenser lens currents 
I, and the condenser lens strength 1/f. deduced 
from the approximate relation 


1/fe=1/v-[(2ro’/re’’) +1]—1/ue, 


where u.=distance from condenser lens to the 
source in the electron gun = 160 mm with the self- 
biased gun. The condenser lens current J, was 
measured on the panel test meter which has 100 
scale divisions corresponding to 25-ma full scale. 
The relations and graphs given below are in terms 
of scale divisions in order to simplify their use. 

It was found that the condenser lens strength 
was given very accurately by Busch’s formula 
where 


1/fe=2-51X10-5-I2, 


and f, is in mm, J, in test meter scale divisions. 
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EXPERIMENTAL RESULTS 


In presenting the experimental results some con- 
sideration must be given to the choice of inde- 
pendent variable. The filament to focusing aperture 
distance uz is the quantity which is preset and 
controls the behavior of both the self-biased and the 
zero bias guns. When operating the self-biased gun 
the total beam current leaving the gun is the 


quantity most readily observed in the laboratory. , 


The relation between these quantities is shown in 
Fig. 2. In operating the zero bias gun the optimum 
filament height is ascertained by observations of 
the nature of the illumination at the specimen.! 
With this gun the filament height controls the 
distance of the effective source from the condenser 
lens, and consequently affects the condenser current 
required to focus an image of the source at the 
specimen, Fig. 3a. It will be seen from Fig. 3b that 
within the range of filament heights normally em- 
ployed the self-biased gun produces an effective 
source at a nearly constant distance from the con- 
denser lens. 

In the following sections the filament height will 
be treated as the independent variable. The varia- 
tion of the beam current 7, arriving at the specimen 
with filament height is given in Fig. 4. In Fig. 5 the 
current density, oo at the center of the focused 
illumination is shown as a function of filament 
height. It may be mentioned in passing, that the 
focusing is quite critical. Moreover the maximum 
current density attainable is dependent upon the 
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Fic. 3a. Condenser current required to focus illumination for 
different filament heights. Full line zero bias gun, dashed line 
self-biased gun. 
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shape of the filament being larger for filaments 
which are sharply bent at their apexes. 

The illuminating efficiency defined as oo/Ip (see 
Fig. 6) is of interest because a high efficiency indi- 
cates that good illumination can be obtained with 
reasonably low values of Jz, which in turn implies 
that a lower filament temperature, giving a longer 
filament life, can be used than if the efficiency were 
poor. This argument is not complete. The self- 
biased gun is operating closer to cut-off so that 
somewhat higher filament temperatures are re- 
quired to yield the same total beam current Jz. 
The filament of the self-biased gun is usually run 
at a higher temperature than that of the zero bias 
gun, but the difference in temperature and life of 


the filament is not large enough to be objectionable * 


in view of the higher intensity of illumination 
obtained. It is here assumed that the filament tem- 
perature is not raised above the value required to 
produce saturation of the beam current. 

Since for visual focusing it is desirable to obtain 
good illumination of that part of the specimen which 
is being observed, and at the same time have small 
total illumination of the specimen to minimize con- 
tamination the variation of the quantity oo/iz with 
filament height is of interest. This is shown in 
Fig. 7. It may be concluded that it is desirable to 
operate the self-biased gun with a filament height 
of 4.8 mm and choose the condenser diaphragm 
radius to provide the proper level of maximum 
illumination for the work in hand. For the standard 
RCA biased gun the optimum filament height will 
be slightly less. 

The quantity oo/t, is obviously connected with 
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Fic. 3b. Effective source positions relative to condenser lens. 
Full line zero bias gun, dashed line self-biased gun. 
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Fic. 4. Current reaching specimen for different filament 
heights. Full line self-biased gun, dashed line zero bias gun. 


the area of the specimen which is illuminated by 
the focused beam. The relation between oo/ig and 
ro’ the radius at which ¢=a/2, is also dependent 
upon the distribution of intensity in the focused 
illumination. This distribution is critically depend- 
ent on the accuracy of focusing of ‘the illumination. 
Finally the illuminated area is usually elongated 
with a length approximately 50 percent greater 
than the breadth in the best cases when the filament 
tip is sharply pointed. Figure 8 shows the variation 
of ‘spot radius’ r.’, with filament height. The radius 
recorded here is the minimum radius of the ‘‘spot.”’ 
The points at which the intensity in the illumination 
fell to half the maximum value were found from 
photo-densitometer traces. The visual observation 
of the spot at low magnification on the final viewing 
screen is not a satisfactory method of determining 
the spot size since the light output of the fluorescent 
screen is not a linear function of the incident current 
density but shows saturation at higher current 
densities. Visual observation therefore leads to an 
over estimate of the half-width. 

The factors discussed above pertain primarily to 
the electron gun. With the results there obtained it 
is a straightforward matter to calculate by the 
formulae of von Borries and Ruska® the dependence 
of the aperture of illumination a, and the intensity 


6 B. von Borries and E. Ruska, Zeits. F. Tech. Physik 20, 225 
(1939). 
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of illumination ¢ upon the condenser lens current J,. 
These are presented graphically in Figs. 9 and 10. 
In calculating the values of a, the minimum spot 
radius r,’ has been used. Again in Fig. 10 the curves 
are not shown flat topped. Flat topped curves would 
only be obtained if the source were circular and 
had uniform current output across its cross section. 
The maxima of the curves in Fig. 10 have been 
adjusted to coincide with the previously observed 
values of oo. These are in both cases somewhat 
greater than would have been calculated from the 
formulas of von Borries and Ruska. This is again 
due to the way the intensity is distributed in the 
focused illumination. 


DISCUSSION 


In the graphs above, the normal operating points 
have been circled. It will be seen that the self- 
biased gun produces a source which is approxi- 
mately one-half the diameter (or one-quarter the 
area) of that produced by the zero biased gun. At 
the same time it provides from ten to fifty times 
the intensity depending on the filament height and 
the shape of the filament apex. This means that at 
the same intensity of the final image, focusing and 
photography can be done at an angular aperture 
of the illumination which is less, by a factor between 
five and seven, than that required with the zero 
bias gun. The total bombardment of the specimen 
is some twenty times greater with the self-biased 
gun. When the illumination is focused, the small 
size of the source tends to minimize the temperature 
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Fic. 5. Current density at specimen. Full line self-biased gun, 
dashed line zero bias gun. 
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Fic. 6. Illuminating efficiency (¢0/Jg). Full line self-biased gun, 
dashed line zero bias gun. 


rise in the specimen, the higher intensity is accom- 
panied however by a higher rate of contamination. 
When the angular aperture is reduced, the heating 
takes place over a larger area, much of it now in 
the mesh. The same is true of the contamination, 
the rate of deposition in the field of view is reduced, 
but the deposition becomes more widespread. 

The above results, together with experience 
gained in its use, make it possible to assess the 
general performance of the self-biased gun. Its most 
notable feature is that it permits focusing and pho- 
tography at low angular aperture with a consequent 
sharpening of the contour and phase contrast 
phenomena. Moreover the high available intensity 
and inherent crispness of the images facilitate the 
detection of defects which may occur in the imaging 
system of the microscope. The next most attractive 
feature is that the setting of the filament height, 
filament centering, and gun tilt are not critical. 
The source is definitely single so that double images 
are not produced. The single source has great 
advantages in electron diffraction studies. 

For a given filament temperature the optical 
behavior of the gun will be the same whether the 
bias is produced by a cathode resistor or a battery. 
The use of a resistor is obviously simpler, and has 
the further advantage that saturated beam current 
operation is attainable. The latter provides a stable 
electron current through the instrument which re- 
sults in a stable distribution of charge throughout 
the imaging system, and a consequent freedom of 
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the final image from shifts which might otherwise 
be caused by varying electrostatic deflection of the 
imaging beam. 

The major disadvantages to be cited against the 
self-biased gun are the heating and contamination 
caused by the higher current reaching the specimen. 
The heating effects will be found of negligible im- 
portance, however, at the intensities and angular 
apertures which are proper for focusing. The rate 
of contamination of the objective lens and its dia- 
phragm can be reduced by using a smaller con- 
denser diaphragm. Finally the higher operating 
temperature of the filament results in a slightly 
shorter filament life, the decrease is not great enough 
to be serious. These disadvantages are far out- 
weighed by the advantages cited above, moreover 
as indicated their effects can be minimized by the 
proper handling of the electron microscope. 


THE OPERATION OF A SELF-BIASED GUN 


To obtain the best results when using the self- 
biased gun attention must be given to several 
factors involved in its operation. The filament tip 
should be bent as sharply as possible (see Fig. 18 of 
reference 1). The tip of the filament should be 
centered with respect to the shield when the fila- 
ment is installed. Any bending of the filament must 
be done before it is heated since heating makes the 
tungsten wire brittle. The behavior of the self- 
biased gun is not critically dependent on filament 
centering so that it is sufficient to center the fila- 
ment by eye. 

For stable operation the illuminating system 
must be clean and must be designed so that insu- 
lating surfaces which may charge up cannot deflect 
the illuminating beam. In the model EMU elec- 
tron microscope the chief insulating surfaces re- 
quiring attention are the fluorescent screen below 
the anode and the window through which it is 
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viewed. Since in using the self-biased gun this 
fluorescent screen serves no useful purpose, it is 
better to clean it down to the metal with steel wool. 
The wire screen over the viewing window should 
not be in contact with the window and should be 
renewed whenever it becomes blackened by con- 
tamination. Alternatively the window may be re- 
placed by a metal plate. 

With the exception of the condenser polepieces all 
metal surfaces in the illuminating system should be 
cleaned with the finest grade steel wool until they 
shine. This is particularly true of the inside of the 
focusing shield near the aperture, and those other 
parts of the system such as the condenser aperture 
insert which are in the direct path of the beam. 


‘Strands of steel wool left behind after the cleaning 


should be carefully removed with a bar magnet. 
In some places this cleaning will eventually remove 
the plating. This will not affect the performance of 
the gun provided that no partially attached pieces 
of plating are left in the gun. The sharp edges of 
such plating might lead to field emission during 
operation. 

The insulating surfaces in the electron gun should 
not be struck by high voltage electrons for in this 
case they become partially conducting and the gun 
eventually breaks down. With the self-biased gun 
it is also desirable to widen the hole through the 
anode so that electrons do not strike any metallic 
surface there. If this precaution is not observed a 
conducting layer may build up on the glass of the 
electron gun and eventually lead to break down. 

The external cleanliness of the electron gun and 
other surfaces which are at high voltages is im- 
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Fic. 8. Spot size variation with filament height. Full line self- 
biased gun, dashed line zero bias gun. 
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Fic. 9. Variation of angular aperture of illumination with 
condenser current. Full lines self-biased gun, dashed lines zero 
bias gun. 


portant since dust will cause corona discharges and 
sometimes spark breakdown. In cleaning the glass 
of the electron gun care should be taken to avoid 
damage to the glaze of the external surface. If there 
is much dust on it, it should be wiped circum- 
ferentially rather than longitudinally. Longitudinal 
scratches in the glass tend to accumulate continuous 
lines of moisture which serve as a track for break- 
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Fic. 10. Variation of intensity of illumination of specimen 
with condenser current. Full lines self-biased gun, dashed line 
zero bias gun. 
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down over the surface. Deposits on the inner face 
of the glass will usually prove intractable if treated 
with organic solvents. They may be removed by 
boiling the glass (without the end rings) in a strong 
caustic solution. 

The adjustment of the zero bias electron gun for 
optimum performance has been described by Hillier 
and Baker. In the case of the self-biased gun it 
will be seen from Fig. 5 that the most favorable 
ratio of intensity of illumination to total current 
reaching the specimen is obtained with the filament 
set rather high with respect to the shield. If the 
filament is set to give a total beam current of 250 
microamps at saturation (—u,=4.5 mm) then the 
maximum intensity of illumination with a 20-mil 
diameter aperture in the condenser is of the order 
0.1 amp./cm? and the area of illumination approxi- 
mately 10X16 microns. It is worth repeating that 
the attainment of high maximum intensity and 
small source size is dependent on using a filament 
with a sharp bend at the apex. Residual asym- 
metries (‘‘wings’’) on the illumination can usually 
be removed with the tilt adjustment on the elec- 
tron gun. 

For some samples an intensity of illumination of 
0.1 amp./cm? will prove destructive. The operator is 
not, of course, obliged to peak the illumination. 
However, if much work is to be done with such 
samples it is convenient to employ a condenser 
aperture of smaller diameter so that the maximum 
illumination will be less intense. When employing a 
smaller condenser aperture the alignment of the 
microscope becomes a slower procedure. It is there- 
fore a simplification if good alignment has been 
previously achieved using a larger aperture and if 
the condenser lens is not moved when the smaller 
aperture is inserted. In all cases the condenser 
aperture should be accurately coaxial in the con- 
denser lens. 

If a very large condenser aperture is used, or if 
the condenser aperture assembly is removed and if 
the microscope is then focused at low magnifica- 
tion, a relatively large fraction of the total beam 
current will reach the final viewing screen. Under 
these conditions the output of x-rays from the final 
viewing screen may rise above the generally ac- 
cepted tolerance level. For this reason it is strongly 
advised that the use of condenser apertures greater 
than the standard aperture supplied with the instru- 
ment be avoided. If the nature of the work involves 
the use of a larger aperture, the x-ray level should 
be measured before proceeding with the experi- 
ments. 
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Point-to-Plane Corona Onsets* 


WruiaM N. ENGLisH** AND LEONARD B. LOEB 
Department of Physics, University of California, Berkeley, California 
(Received December 20, 1948) 


The effect of point material and point radius of curvature on positive and negative intermittent 
corona onset potentials has been studied with a point-to-plane gap in air at atmospheric pressure. The 
negative Trichel pulse onset strangely is independent of point material but does depend on point 
history and radius. This surprising result is shown to come from the circumstance that the Trichel 
pulse onset depends on current densities needed to condition the point surface and yield a higher value 
of the second Townsend coefficient to give the increased currents. Trichel pulse onset thus-does not 
mark the onset of a self-sustaining discharge. The self-sustaining discharge initiates at lower potentials 
and leads to currents of a low order until the cathode spot cleans up. 





N a survey article “Recent developments in 
analysis of the mechanisms of positive and nega- 
tive coronas in air,’’! Loeb*** discusses the rather 
surprising equality of positive and negative corona 
onset potentials in air, as contrasted with a much 
lower negative onset in pure He and Ne. He con- 
cludes that this equality must be due to an increase 
in the work function of the point, corresponding to a 
decrease in ‘“‘“gamma,”’ the secondary electron emis- 
sion coefficient, which happens to raise the negative 
onset to just about the same value as the positive. 
With this dependence on the ‘‘gamma”’ of the point, 
one would at first expect that materials of different 
work function would show different negative corona 
onset potentials. Indeed, experiments with a water 
drop point? had revealed that here was a substance 
with such a low “gamma”’ that apparently negative 
corona could not take place at all! The purpose of 
the present work was to look for any solid substances 
which might exhibit such an effect, and this involved 
also an investigation into other factors affecting the 
onsets. A discussion of the nature of the onsets, and 
of the types of pulses observed, will be found in a 
previous paper.® 


APPARATUS AND TECHNIQUE 


The arrangement of the point-to-plane gap is 
shown in Fig. 1. 

The intermittent corona onset potential, Vg, was 
measured by means of a potentiometer calibrated 
against a standard cell. For positive point, Vg+ was 
the lowest potential at which burst pulses could be 
seen on the Dumont 241 oscilloscope, with full gain 
(100). For negative point, Vg- was the lowest 
potential at which Trichel pulses could be obtained. 
The observed onset potential for both polarities is 
considerably affected by the amount of “‘triggering”’ 





* This work was done under contract with the ONR. 

** Now at the National Research Council, Chalk River, 
Canada. 

1 L. B. Loeb, J. App. Phys. 19, 882 (1948). 

*** In a recent article, J. D. Craggs and J. M. Meek, Proc. 
Phys. Soc. London 61, 327 (1948). 

*W. N. English, Phys. Rev. 74, 179 (1948). 

*W. N. English, Phys. Rev. 74, 170 (1948). 
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ionization which is introduced into the gap in order 
to start the electron avalanches. A “‘natural’’ onset 
will be defined as that obtained with the ionization 
from cosmic rays and the usual slightly radioactive 
surroundings. 

Except where otherwise indicated, the gap geome- 
try was as recommended by Loeb,‘ with a ratio of 
point diameter to gap length of 1/80, and plane 
radius greater than gap length. The work was done 
in room air. 


EXPERIMENTAL RESULTS 
Triggering Ionization 


The dependence of positive intermittent onset 
potential on the amount of triggering ionization is 
marked, since nearly all the negative ions produced 
in the gap drift into the high field region, and the 
electrons formed undergo multiplication by collision 
with gas molecules. Thus, a radioactive source must 
be used with caution. The practice followed has been 
to place a weak 8-y-ray source (50 grams of thoria- 
nite ore) at such a distance as to give a 3- to 5-fold in- 
crease in the frequency of burst pulses at onset. This 
greatly facilitates rapid and accurate determination 
of Vgt without lowering it appreciably. Much more 
triggering than this causes the burst pulses to merge 
into one another on the oscilloscope screen, and 
produces a lowering of the onset through cumulative 
avalanches. 

With negative point, unfortunately, the ‘‘natural’”’ 
onset cannot be determined very accurately with 


Hvt 





Fic. 1. Diagram of the 
corona gap, showing the point 
of insertion of the oscilloscope 
leads. 
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4 J. M. Parker and L. B. Loeb Phys. Rev. 74, 1557A (1948). 
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Fic. 2. Characteristic results from the measurements of 
negative and positive corona onset thresholds, using points 
made of the materials indicated. The numbers and pointer lines 
represent different numbered metal points of the series used in 
the measurement. 


the oscilloscope, owing to the long and randomly 
varying time between successive Trichel pulses. 
This arises from statistical fluctuations in the rather 
complex ionizing sequence, from the low and per- 
haps variable secondary electron emission coeffi- 
cient of the point, and from the inherently minute 
sensitive volume for multiplication of electrons, 
since here the negative carriers move out of the high 
field region, as well as from the limited sweep time 
of the oscilloscope. 

The same weak §-y-ray source mentioned earlier 

produced an increase in the number of Trichel 
pulses at onset (and a lowering of the onset) for 
large negative points (diam. 2 mm), although the 
effect was of course much less than for the corre- 
sponding positive point. It was somewhat surprising 
to find an appreciable effect also with even the 
smallest point (diam. 0.03 mm). This was probably 
a “‘counting”’ action due to positive ions from elec- 
tron tracks. Although it is doubtful if there is much 
danger of falsifying the negative onset through a 
cumulative ionization effect, only the minimum 
amount of triggering necessary to give a reasonably 
sharp and reproducible onset was used for onset 
measurements. 
’ An attempt was made to trigger the negative 
corona by illuminating brass, platinum, and zinc 
points with ultraviolet light from an Hg arc. Reli- 
able triggering was not obtained owing to the 
susceptibility of the surfaces to change in work 
function and photo-emission connected with illumi- 
nation in air. Sometimes illumination caused a 
marked increase in the number of pulses at onset ; at 
other times there was little effect, or even a notice- 
able decrease. Similar results have been reported by 
Loeb et al. 





® Loeb, Kip, Hudson, and Bennett, Phys. Rev. 60, 714 
(1941). 
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Point Material 


‘Kip and Bennett’ found no measurable difference 
in negative onset for points of platinum, copper, 
aluminum, iron, and various alloys. In view of the 
experimental uncertainties in the determination of 
negative onset, and of the all-important role played 
by the secondary electron emission coefficient of the 
point, it was essential to confirm this result, using 
the positive onset, which is independent of point 
material, as a check. Accordingly Vgt and Vg~ were 
determined (for ‘‘natural” ionization), using 1-mm 
diameter points of carbon, zinc, copper, aluminum, 


_ and platinum. 


Points were ground hemispherical and carefully 
polished, but no special effort was made to get 
identical radii of curvature since the difference be- 
tween Vgt and Vg~ gives the information needed. 
For each determination of the quantity Vg~— Vgt, 
four measurements of onset (+ — + —) were taken. 
Some of the results are shown in Fig. 2. 

As found in all previous experiments, positive 
onset was reproducible and apparently independent 
of point history or material.t The increase in po- 
tential, from burst pulse onset Vg*t to the first 
appearance of streamers, was essentially constant 
and averaged 270 volts. : 

Negative onset as determined was higher than the 
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Fic. 3. The onset of positive burst pulse, Vg*, and of nega- 
tive Trichel pulse corona, Vg~, with and without triggering 
radioactive sources as a function of point diameter in mm. The 
decrease in the negative potential relative to the positive onset 
potential with decreasing radius is clearly seen. 


t+ For example, an onset determination made with the 
platinum point 9 months earlier checked exactly, although the 
point had been used considerably in the meantime. 
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positive by from 125 to 1000 volts.f It was notably 
dependent on point ‘‘condition” varying from one 
measurement to the next, and even during a de- 
termination. This was especially true for Pt. Tests 
showed that the effect was not due to changes in the 
gas filling of the chamber. Taking into consideration 
the uncertainty in the negative onset, for a given 
point, there was no significant difference in behavior 
attributable to point material. 
The following side experiments were done: 

(1) a A copper point was oxidized black in a bunsen flame. 


b Point was etched with dilute nitric acid. 
c Point was oxidized again in bunsen flame. 


There was no change in positive onset after a, b, 
or c, and no significant difference in the negative 
onset, except that the etched point gave more pulses 
and a cleaner onset. 

(2) a A platinum point was bombarded with positive hydro- 
gen ions at 1.5 cm Hg pressure of hydrogen, current 10— 

20 wa for 65 minutes. 
b Hz was pumped out, air admitted, and negative and 
positive onsets immediately determined. There was no 


significant difference from the values previously ob- 
tained. 


Point Size 


Some early measurements with points of FeS, 
CuO, and As had given Vg~ considerably below Vgt. 
These were necessarily rough points because of the 
nature of the material. The results with the 1-mm 
points at once suggested that this lowering of Vg- 
was due to point shape. Onsets were then determined 
for a fine steel needle point, a fine hemispherical 
platinum point (diam. 0.03 mm), and one of the 
1-mm carbon points with the tip broken off to give 
sharp edges. For such fine points the negative onset 
was below the positive onset by several hundred toa 
thousand volts! 

The next step was to study the variation of onset 
potential with point radius of curvature. A series of 
12 polished hemispherically tipped platinum points 
(diam. 0.03 to 2.4 mm) was made and measured 
under 100 to 4400 power magnification. The ratio of 
point diameter to gap length was kept constant at 
1/80, as previously indicated. At least three meas- 
urements were made for each point, Vg-, Vgt, Vg-, 
negative onset being determined both for ‘‘natural”’ 
ionization, and with the radioactive source at the 
maximum distance which gave a marked increase in 
pulses near onset. Results are shown in Fig. 3. The 
consistent character of the data indicates a real and 
smooth lowering of the Trichel pulse onset relative 
to the positive onset as point radius decreases. 





} Similar results have been observed with concentric 
cylinders by C. G. Miller (see reference 8) in air down to 90 mm 
pressure. 
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Fic. 4. The various thresholds for positive point corona as a 
function of the gap length. Vg* is the beginning of intermittent 
burst pulse onset and the beginning of the positive point 
Geiger regime. Vstr.* represents the appearance of what are 
called pre-onset streamers. Vo* represents the disappearance of 
streamers and the onset of the continuous burst pulse type 
corona, marking the end of the Geiger counter regime for the 
point. It is noted that the counting plateau, which is pre- 
sumably somewhat shorter than the vertical distance between 
Vg* and Vo", increases with gap length, and that the various 
thresholds are more clearly delineated the longer the gap for 
the same diameter point. 


Gap Distance 


Although not of fundamental importance, a plot 
of the three positive onsets versus gap length is given 
in Fig. 4 to show the type of variation. The benefit 
obtained in using relatively large gaps, such as 
recommended by Loeb, so that the onsets will be 
well separated, is evident. Negative onset would be 
very close to Vg* for a point of this size, as can be 
seen from Fig. 3. 


DISCUSSION OF RESULTS 


By specifying the level of triggering ionization, in 
terms of ‘“‘natural’’ ionization or otherwise, it is 
possible to obtain reproducible positive onsets from 
the oscilloscope, although even here the amplifica- 
tion used has some effect on apparent onset. For 
negative onset, a large difference between apparent 
oscilloscope onset and the “break’’ in the current- 
voltage characteristic may be found. Thus, in one 
case, apparent oscilloscope onset was 900 volts 
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above the change in slope of the current voltage 
curve. However, on lowering the voltage to corre- 
spond to the break in the curve, and waiting ten 
minutes, a single Trichel pulse was obtained on the 
oscilloscope. The significance of the abrupt current 
increase and of Trichel pulse onset is discussed later. 
A more satisfactory method of triggering a negative 
point would greatly facilitate the study of negative 
corona. The lack of success with ultraviolet light in 
this regard is to be expected from gas-covered and 
oxide-coated surfaces and is behavior familiar to all 
workers on photoelectric effects in air. 

In theory, the onset of a self-sustaining discharge 
from a negative point is given by y expQ’ adx =1, 
where a is Townsend's first coefficient, dx is the 
distance from the center of the point, r its radius, a 
is the distance from the point at which Townsend's 
coefficient loses its efficiency in the declining field. 
y is the chance of emission of electrons from the 
cathode per ion created in the gap, and could be 
ascribed to a photo-effect at the surface, or more 
probably to electron emission by positive ion bom- 
bardment. Thus y¥ is critically dependent on condi- 
tion of the surface and on the work function, which 
must be less than the ionization potential of the 
positive ions. From the nature of the current- 
potential curves for the negative point, it had been 
assumed that the Trichel pulse appearance set the 
threshold defined above. 

For the positive burst pulse threshold, Loeb® has 
shown that #f exp JZ” adx = 1, where f is the fraction 
of the collisions producing ions that give photons 
capable of ionizing in the gas, and 8 is a geometrical 
factor dealing with the chance of absorption of such 
a photon beyond a in the gas. Here onset depends 
entirely on the properties of the gas and not on the 
surface. In glow discharges at lower pressures, and 
in studies of Townsend's coefficients at low pres- 
sures, y is of the order of 10-? to 10-*. It has never 
been evaluated at atmospheric pressure in gases, but 
owing to back diffusion’ of emitted electrons alone, 
it could be smaller at 760 mm. Another factor that 
cannot be neglected is that y is dependent on the 
energy of impact of positive ions on the surface. No 
data below positive ion energies of 20 volts exist. 
‘For air, 8f has the value of about 10-*. Thus, it is 
clear that if the Trichel pulse onset evaluates y in 
air at atmospheric pressure, y must be of the order 
of 10~*, while in pure nitrogen and otherwise, it can 
_ well be 10- to 10-*. 

However, this is not the worst difficulty the 
present observations present. They indicate that the 
Trichel pulse onset determined value of y is inde- 
pendent of the metal or the metal oxide surface. 


°L. B. Loeb, Phys. Rev. 73, 798 (1948). 

7L. B. Loeb, Fundamental Processes of Electrical Discharge in 
re eae Wiley & Sons, Inc., New York 1939), pp. 315 
an ; 
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This independence of y for the widely different 
metals or their oxides was tentatively ascribed to 
the presence of some common gaseous film on all 
metal surfaces, involving oxygen or nitrogen oxides. 
That some sort of an adsorbed film does form in the 
presence of oxygen which is common in properties, 
and behavior on all surfaces cannot be denied, and 
copious evidence for the existence of such films 
exists. However, it is hardly probable that such a 
film of gas could alter the value of y to the extent of 
making it completely independent of the metal, as 
different metals might have films of different y. 

As long as the work was confined to studies of the 
corona with point-to-plane geometry, the character 
of the observed onset curves was such that one could 
hardly avoid concluding that the Trichel pulse onset 
was the actual threshold of a self-sustaining negative 
discharge, and that it fixed the value of y. The 
difficulty in a more discriminating decision here, lies 
in the instrumental difficulties of observing the 
exact point of Trichel pulse onset relative to the 
very sharp increase in current which seems to mark 
the appearance of the pulses. The recent work of C. 
G. Miller,® using concentric cylinder geometry with 
a fine central wire, and studying gases ranging from 
pure nitrogen through various mixtures of oxygen 
and nitrogen to pure oxygen, has disclosed a set of 
circumstances which now modify the picture. With 
point-to-plane corona and relatively small points, 
pre-onset currents remain at about 10~ ampere 
until a potential is reached at which there is an 
abrupt change of current to about 10-7 ampere with 
the appearance of Trichel pulses in the steeply 
rising portion. With concentric cylinders, depending 
on the character of the gas, Miller observed that, in 
general, the rise of current with potential up to the 
point where Trichel pulses occurred, was far more 
gradual than with the point. Thus, the jump of cur- 
rent on the appearance of the Trichel pulses was less 
sudden. The Trichel pulses also first appeared well 
up on the steep rising portion of the current voltage 
curve. His results with guard ring cylinders indi- 
cated that there was a definite and clearly defined 
pre-onset current along the whole wire in contrast to 
localized pulses, which in pure oxygen at 760 mm 
was half as large as the total current when the 
Trichel pulses appeared. In other mixtures, while 
the rise was more abrupt, it was less abrupt than any 
observed with point-to-plane corona. This indicates 
that well below the threshold of Trichel pulse onset, the 
ordinary field intensified gaseous discharge has been 
replaced by a self-sustaining Townsend discharge of 
low current value along the whole length of the wire. 
Since this Townsend discharge is below an intensity 
where a glow is visible, and since it develops in a 


* C. G. Miller, Ph.D. thesis, University of California (1948); 
Phys. Rev. 75, 1460 (1949). 
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current range which makes it inaccessible to oscillo- 
scopic observation, the transition from the inter- 
mittent Townsend discharge to the steady self- 
sustaining Townsend discharge evaluating y cannot 
be indicated even with Miller’s set-up. It is, how- 
ever, considerably below Trichel pulse onset in 
potential as judged from the magnitudes of Miller’s 
currents. Thus, the onset of a self-sustaining dis- 
charge. which defines the value of y, lies at po- 
tentials far below Trichel pulse onset. Hence, the 
onset and thus, y may quite definitely be a function 
of electrode material as well as of gas coating. Until 
the true onset can be observed, this point will not be 
decided. Trichel pulse onset, therefore, does not repre- 
sent a threshold defining a y. With a point-to-plane 
geometry, the Townsend discharge current, owing 
to the limited area of the point, would be so small 
that its contribution to the current would hardly be 
noted. Thus, the abrupt increase of current near 
Trichel pulse onset was naturally but mistakenly 
assumed to represent the threshold for self-sus- 
taining discharge. 

This conclusion permits immediate interpretation 
of the present observations. The onset of the self- 
sustaining Townsend discharge, fixed by the true y 
of the surface, leads to a type of glow discharge 
whose value is limited by the relatively low value of 
y. The currents in this limited discharge increase in 
a regular and consistent fashion with increase in 
potential, since the first Townsend coefficient is in- 
creasing irrespective of any small increase in y. It 
appears, however, that with the low y from the gas- 
coated surface, the secondary electron emission is so 
low that the positive ion space charge created by the 
secondary electrons coming from the point is not 
sufficient to initiate the autocatalytic build-up of 
ionization described by Loeb.! However, as the 
potential is raised, the surface of the point or wire is 
being bombarded by greater densities of positive 
ions and possibly with positive ions of increasing 
energy. At higher positive ion energies, such bom- 
bardment leads to sputtering of the point material ; 
here it leads to removal of the gas and/or oxide 
films. There is, then, at each current, an equilibrium 
process set up in which there is competition between 
the removal of the gas film by positive ion bombard- 
ment and the tendency of the gas film to re-form by 
the diffusion of oxygen atoms and molecules to the 
point. Eventually a current density must be reached 
at which the surface of the metal becomes denuded 
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of its gas and perhaps even of its oxide layers, i.e., 
the point is cleaned up or “‘conditioned.”’ At this 
juncture the work function of the surface is very 
much decreased, and the value of 7 may increase by 
an order of magnitude, i.e., from 10-* to 10-*. 

The value of the threshold current density for 
clean-up will depend primarily on the character of 
the gas, the percentage of active molecules, and the 
energy of the incoming positive ions. Therefore, the 
higher the field, the more concentrated the dis- 
charge is by the geometry, or the lower the concen- 
tration of oxygen, the lower will be the current 
density at which the film will be removed. 

With the film suddenly dispersed at a given po- 
tential, the value of y is such that the potential is 


‘far above that for the threshold value of any sort 


of a discharge of this character, i.e., the point is 
highly overvolted. Thus there will be a secondary 
liberation of electrons from the overvolted surface 
which will produce the autocatalytic increase in 
current which was described by Loeb.' This leads to 
a burst of current which can either go to arc dis- 
charge as in pure ‘nitrogen with no limiting resistor, 
or which can, in a gas like air, create so many 
electrons in an explosive fashion that negative ion 
space charges choke off the discharge producing the 
Trichel pulse. It is seen that this mechanism com- 
pletely accounts for the Trichel pulse onset. Through 
the differentiation of Townsend discharge onset and 
Trichel pulse onset, this picture disassociates the 
latter from the evaluation of y and makes possible 
a variation of y with point material. The mechanism 
proposed accounts for higher pre-onset currents 
with other geometries. It accounts for the relative 
lowering of Trichel pulse onset with point radius 
through the increased current densities and higher 
tip fields. It accounts for the aging effects with Pt 
and other surfaces, and for the more gradual in- 
crease of pre-onset current with potential before 
Trichel pulse onset in pure O2 where the O, reforms 
the surface layer.6 In pure Ne, and Ne with small 
amounts of Os, it explains the sudden current in- 
creases and change of discharge form from a diffuse 
to a concentrated discharge when points are initially 
subjected to constant current for some time at the 
same potential.® It finally indicates the need for 
much more work on the location of the true 
Townsend threshold and the evaluation and study 
of y. 


°G. L. Weissler, Phys. Rev. 63, 96 (1943). 
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Magnetic Properties of Iron Compacts in Relation to Sintering Temperature* 7 
di: 
ROBERT STEINITZ 
American Electro Metal Corporation, Yonkers, New York sir 
(Received January 3, 1949) va 
wi 
Previous experiments have indicated that the magnetic permeability of sintered iron compacts is to a . 
large extent determined by the final density of the compact. The permeability of five different iron powders P 
has proved to be independent of the origin of the powder, if the compacts were sintered at the same tempera- 
ture. The results were compared with the theory of Polder and Van Santen and agreed very well with their ga 
predicted curve, if the residual pores were assumed to be flat disks. It could be expected that at higher TI 
sintering temperatures, the disk-shaped pores would change into a more equiaxed or spherical shape. 12 
Experiments to check this theory have been made, and the results are in good agreement. Compared for . 
identical densities, the permeability of different iron powders is appreciably higher if sintering is done at co 
1250-1350°C and for 24 hours, instead of at 1150°C for one hour. Different powders, however, now show 
quite different permeabilities, which could be explained by a different inclination of the powder to form mi 
spheroidal pores. We 
It is shown that coining after high temperature sintering brings the permeability ¢alues back to the = 
curve of the disk-shaped flat pores. The effect of the increase in permeability at higher temperature cannot, 
therefore, be due to a purification process, but must be attributed to a change in the pore shape. wi 
a 
su 
. a tae , " , , oe of 
ECAUSE of the limitations of commercial furnaces, pact is a variable, very well suited to follow this sinter- Al 
iron parts for magnetic applications are usually ing behavior. Polder and Van Santen! have published a i 
sintered at a temperature between 1000° and 1150°C. theory on the permeability of porous compacts with 
Furthermore, it is possible to hold dimensions to close special reference to the influence of the shape of the 
tolerances at this sintering temperature, as shrinkage is pores. It was shown in a previous paper’ that iron com- 
only insignificant. It is well known, however, that, at pacts, sintered at around 1150°C for about 1 hr., have a 
only 100°C higher temperature, dimensional changes relationship between maximum permeability and den- om 
and densification already become much more pro- sity corresponding very closely to the one predicted by tel 
nounced. It can be expected that this increased sintering the Dutch authors for flat holes. Such a flat disk-like for 
activity at higher temperatures manifests itself not only pore shape is to be expected in compressed parts, if fre 
in an increased and accelerated shrinkage, but also has__ sintering temperature and time are not high enough to © = 
an influence on the physical properties of the compact. change it. However, results published elsewhere indi- bata 
This should be especially true for properties which de- cated that an increase in sintering temperature results _ 
pend on the shape and distribution of the pores. The in a considerable improvement of the permeability, ae 
modern theory of sintering postulates that pores change though no systematic studies of this effect have been pe 
to a spherical or almost spherical shape before they published. th 
disappear completely. A comparison of properties ob- The purpose of this paper is to investigate the po 
tained after different sintering temperatures should permeability-density relationship for higher sintering = 
indicate whether such a spheroidization of pores actually temperatures. It is, therefore, a continuation of work im 
takes place, provided the comparison 1s made for identi- previously published,? but the present results can still me 
cal densities. This would eliminate the influence of the be considered only as preliminary, due to a lack of ~ 
amount of porosity or the pore volume, and make the sufficient data. Further work in this line is in progress. - 
pore shape the determining variable. ' The results obtained, however, indicate satisfactory of 
The maximum magnetic permeability of an iron com- agreement with theory. sit 
TABLE I. Characteristics of i rd educed ders). 
Pisincnutead aracteris mae iron powders (r uced powders) ; MATERIALS AND PROCED S . 
No. 1 No. 2 No. 3 No. 4 The raw materials used were four reduced iron 
Apparent density tt 2s th 28 -. powders from different manufacturers. One powder, 
arbon content, % 0.1 1 Ri ' . * ‘ : 
+ ping sf er weve Re © 073 0.43 0.40 oss No. 3, was subjected by the manufacturer to a special 
“Screen analysis process which resulted in almost perfectly spherical po 
+100 0 0 6.5 0 ; isti iven 
~100 +150 175 970 19.0 75 particles. The characteristics of the powders are gi as 
—150 +200 34.5 25.0 16.5 235 in Table I. Experiments with electrolytic powders, pe 
== t= oe = Ry a ‘which were included in the previous investigation,” are pe 
— 325 +e 20:0 3.0 37.5 379 ‘not yet concluded, so that the results are not reported an 
—_—_—— th 
cane 1D. Polder and J. H. Van Santen, Physica 12, 257 (1946). Nx 
* This paper was presented at the International Powder Metal- 2 R. Steinitz, Trans. A.I.M.E., Tech. Pap. No. 2335 (February, to 
lurgy Conference at Graz, Austria, July 12-17, 1948. 1948). 
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here. The powders are all quite similar in particle size 
distribution and other characteristics. 

In order to obtain different densities at identical 
sintering temperatures, the powders were pressed with 
various pressures, ranging from 10 to 50 t.s.i. A ring 
with an average diameter of about 2 in. was directly 
pressed, so that no machining was necessary afterwards. 

Rings proved to be the only type of sample which 
gave reproducible results in magnetic measurements. 
These rings were sintered in hydrogen for 24 hr. at 
1250°C and 1350°C, respectively, and cooled in the 
cooling chamber of the furnace. 

After sintering, the density of the sample was deter- 
mined by weighing it in air and water, after the sample 
was coated with a lacquer to prevent water from enter- 
ing the pores. The rings were then wound with two 
windings, a magnetizing winding of about 120 turns and 
a measuring secondary winding of 25 turns. This re- 
sulted, with a current of 5 amp., in a magnetizing field 
of about 50 oersteds, which was considered sufficient. 
All measurements were done with a ballistic galva- 
nometer. 


RESULTS 


The results of the magnetic measurements of all 
samples are given in Table II, together with sintering 
temperature and time and resulting density. The results 
for a sintering temperature of 1150°C were either taken 
from the previous investigation? or from new measure- 
ments, where new powders were used. In order to get 
well reproducible results, it proved necessary to demag- 
netize with special care. Due to the effect of eddy cur- 
rents, especially in dense materials, the maximum am- 
perage of the demagnetizing a.c. must be higher than 
the maximum d.c. used, and should be decreased con- 
tinuously to a very low value before being switched off. 
The measurement of the maximum permeability on the 
initial magnetization curve is very sensitive to any 
remaining magnetism. But even if all precautions are 
taken, results still fluctuate and stray values with con- 
siderable deviations occur. Another source of error is 
the density determination, as well as the determination 
of length and cross section of the ring sample. As the 
sintering temperature has a pronounced effect on the 
permeability, a uniform temperature in the sintering 
furnace is required. 


DISCUSSION OF RESULTS 


In Fig. 1, the maximum permeability for all four 
powders after sintering for 1 hr. at 1150°C is reproduced, 
as function of the density. (The density is given in 
percent, 7.85 g/cc being taken as the value for 100 
percent dense iron.) The values for powders Nos. 1, 2, 
and 3 are all well on the curve which was taken from 
the previous publication. The one value for powder 
No. 4, however, is somewhat higher ; this is an exception 
to the observations on all other powders, the perme- 
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TABLE II. Magnetic properties of compacts made 
from reduced powders. 











Maxi- 
Coer- mum 
Residual cive perme- Induc- 
Density magnetism force ability tion 
Powder Sintering % Br He Bmax Bu 

No. 1 1150°C 93 11,900 3.5 1590 14,300 

1 hr. 89 11,600 3.5 1600 13,400 

81 7000 4.5 700 9800 

78 6400 4.3 700 8800 

75 6500 4.5 730 8200 

1250°C 90 10,200 1.9 2500 12,800 

24 hr. 88.5 11,400 18 3100 13,800 

83.5 9300 2 2100 11,300 

80.5 8300 2.1 1900 10,600 

79 8500 2 2100 10,200 

1350°C 92 10,800 1.4 3000 13,500 

24 hr. 82 9200 1.8 2000 ‘11,000 

77 7300 1.8 1800 9400 

No. 2 1150°C 90.5 9600 2.8 1470 13,000 

1 hr. 74 4900 3a 680 7000 

1250°C 89 8100 3.3 1350 13,000 

24 hr. 74.5 5500 3.6 850 8700 

No. 3 1150°C 71 5500 3.7 750 7800 
1 hr. 

1250°C 88.5 9500 2.1 2000 13,000 

24 hr. 86.5 8500 3 2000 12,400 

79.5 7500 3.2 1500 10,100 

1350°C 93.5 10,400 1 2600 14,900 

24 hr. 90 8300 2 2100 12,500 

73.5 6150 2 1500 8800 

No. 4 1150°C 80 7500 2.5 1400 11,000 
1 hr. 

1250°C 93 12 250 1.6 3300 15,000 

24 -hr. 86 11,000 1.5 3100 13,200 

79.5 9300 1.9 2350 11,500 

1350°C 95 10,000 1.6 3200 15,500 

24 hr. 86 10,100 1.5 3150 13,000 

79.5 7800 1.8 2500 10,400 








ability values for which are all well represented by this 
curve. A reason for the behavior of this powder is not 
yet known; however, it might be noted that the shrink- 
age in this case is already appreciably higher than with 
the ‘‘normal” powders. It is possible that the sintering 
temperature of 1150°C is already quite high for powder 
No. 4. 

Figure 2 gives the values of the maximum perme- 
ability after sintering for 24 hr. at 1250°C. The points 
for the 4 powders are no longer on one curve; on the 
contrary, they differ appreciably, powder No. 4 having 
about 3 times the permeability of powder No. 2. 
Despite considerable fluctuation, the general trend of 
the curve for each powder can be easily recognized. 
For powder No. 1, this curve gives permeability values 
about 2% times higher than after 1150°C sintering, 
which is a very considerable improvement. 

The sintering at 1250°C was extended for 24 hr., 
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Fic. 1. Permeability-density relationship for sintering 
for 1 hr. at 1150°C, 


while, at 1150°C, the sintering time was only 1 hr. This 
was done to make more pronounced any effect the 
higher temperature might have. However, a conse- 
quence of this procedure is that there is little difference 
between 1250°C and 1350°C sintering. Values for both 
temperatures are practically on the same curve. Some 
kind of an equilibrium appears to have been reached. 
The change of the permeability with time for the same 
’ sintering temperatures is still being investigated. 

The induction value for a magnetizing field of 40 
oersteds, By, was taken as a measure for the saturation. 
It is interesting to note that these induction values are 
not influenced by the sintering temperature. They all lie 
well within the expected fluctuations on one curve, if 
plotted against density, as is true for the residual mag- 
netism B,. No trend with temperature is discernible. 

It was claimed that the improvement in permeability 
is due to the purification by loss of carbon at the higher 
sintering temperatures. That this is not the essential 
effect can be seen from the following: A ring was 
sintered at about 1350°C, then coined and annealed at 
760°C. This resulted in a density of 87.5 percent and a 
maximum permeability of 1530, which is close to the 
curve for 1150°C sintering. (See reference 2, Table II, 
Specimen No. 12.) After the high sintering tempera- 
ture, a much higher permeability would be expected for 
this powder. Any purification would, of course, still be 
effective after coining, but it is believed that the coining 
process flattens the pores again and thereby removes 
any spheroidization effect of the high sintering tem- 
perature. 


714 





5000 





Pon es Ph 
3000 —_ Za a A 
—" / 
> */ 

7 

4 









































bal 
2000 © —14 Z 
enn Wana curve 
: ra Aa 
wen am 
1000 = 
x _——— 
° 
70 75 80 65 90 95 100% 
DENSITY Fe 
° PowoeR”: 
1250 °C ; x POWDER “2 
24h. + POWDER *3 
& POWDER *4 


Fic. 2. Permeability-density relationship for sintering 
for 24 hr. at 1250°C. 


SUMMARY AND CONCLUSIONS 


Permeabilities of iron compacts prepared in the form 
of rings were measured after sintering at 1250° and 
1350°C for 24 hr., and compared with values of samples 
sintered at lower temperatures. A considerable increase 
of the permeability for identical densities was found 
for most of the powders, but values for different 
powders no longer lay on the same curve. 

An explanation for the increase in permeability is 
given by the theory of Polder and Van Santen,! if it is 
assumed that the pores change their shape from flat 
disks to approximately spheroids at the high sintering 
temperature. An attempt to follow this change in pore 
shape under the microscope was not successful. The 
pores are probably too small and too easily changed 
in the polishing procedure. The difference in the be- 
havior of the permeability of the various powders cannot 
yet be completely explained. It is probably not due to a 
difference in chemical composition, but to a greater or 
lesser sintering activity of the powder which is also 
shown by an increased shrinkage. Measurements of 
relative shrinkage and behavior of mechanical as well 
as magnetic properties at different sintering tempera- 
tures and times, always compared for equal densities, 
will help to explain these phenomena and show whether 
the sintering process and the change in pore shape can 
be followed by these experiments. The improvement of 
the magnetic properties by the higher sintering tem- 
perature will also be of commercial importance. 
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The Vaporization of the Cathode in the Electric Arc 


RAGNAR Horm 
Stack pole Carbon Company, St. Marys, Pennsylvania 


(Received March 7, 1949) 


Calculations are carried out which show that, because of the smallness of the cathode spot, the cathode 
is not able to dissipate the heat generated with heavy currents by conduction. Therefore, a much higher 
vaporization of cathode material per coulomb occurs for heavy currents than for low currents. 





I. INTRODUCTION 


HE following computations provide an explana- 
tion for the fact that the vaporization per 
coulomb from the cathode of an arc is greater at heavy 
currents than at low currents. The treatment is based 
on observations and on the assumptions that the tem- 
perature in the cathode spot is of the order of the 
boiling point of the cathode material. All of the obser- 
vations on the arc, which have been used for the 
calculations, are of an approximate nature, and do not 
motivate exact calculations. Therefore, several simplifi- 
cations have been used. One of the aims of this paper is 
to call the attention to a problem where new and more 
accurate observations would make possible important 
conclusions. 
For information about some details, references are 
made to the book, Electric Contacts! where further 
literature is listed. 


Il. THE ENERGY BALANCE AT THE CATHODE SPOT 


We neglect the relatively small radiation from the 
cathode spot, and assume that the total energy pro- 
duced in the cathode drop region is transferred to the 
cathode. 

The gain of energy, W, in the spot is given by the 
product of the current, 7, and a voltage which is the 
difference between the cathode drop, V., and the heat 
of vaporization of the electrons given by the work 
function, g. Hence approximately 


and 


for metals W=(V.— ¢)J= “ 
; (1) 


for carbon W= 15] 


For simplicity, the calculations are carried out for 
metals only in the following derivation. We regard the 
cathode spot as a circle with the radius a and designate 
the current density by 7. Thus 


I=7a’}. (2) 


It is important to notice that 7 is fairly uniform across 
the cathode spot because of the uniformity of the 
temperature. Consequently, 7 in Eq. (2) may be re- 
garded as a constant for the concerned metal and 
current intensity. 


1 Ragnar Holm, Electric Contacts (H. Geber, Stockholm, 1946). 
This book is denoted E.C.H. in the following. 
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The losses of energy per second are the heat W, con- 
ducted away from the spot into the cathode, and the 
heat W, used for vaporization of cathode material. 

The heat W, is equal to the temperature # of the 


- cathode spot times the thermal conductance away from 


it. We regard the cathode as a semi-infinite body. The 
conductance? is then 4aX, where ) is the heat conductiv- 
ity and 


W.=4and. (3) 


The latent heat of vaporization may be H joule/cm* 
and the volume evaporated per second v cm*/sec. Thus 


W.=H2. (4) 
The energy balance requires that 
W=W+W, 
or 
107 = 107a?j= 4ad0+ Hv. (5) 


We see that the radius a of the spot enters with 
different powers in two terms of Eq. (5). This fact is 
the basis for important conclusions. It means that, with 
increasing current J and radius a, the term 107a’j 
increases quicker than the term 4a\#. Finally, at large 
currents, the term 4a\0 becomes negligible, and all the 
heat is used for vaporization. 

We shall follow the calculations through in some 
detail. First we introduce two variables: 


5=v/I, the evaporated volume per coulomb removed 
from the cathode by the arc, and 
x=the fraction of the total energy wh'ch is used for 
vaporization ; 
thus 
sW=W,=Hov 
or 
5= 10x/H. (6) 


We see that the volume v may be replaced by /6 in 














TABLE I. 
Materials Cu Ww ¢ 
0 = 2300 5000 5000 gs 
» = 3.8 1.5 0.04 watt/cm deg 
H = 5 9.2 21 10* joule/cm‘ 











2 See reference 1, p. 23. 
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Fic. 1. Volume 6 of the cathode material vaporized per coulomb 
. through the arc plotted against the current intensity. The number 
written above a plotted point gives the current density 7 amp./cm*; 
the number below the point gives x, the fraction of the energy 
developed in the cathode drop region, which is available for 
vaporization. 


Eq. (5). Then Eqs. (2), (5) and (6) have the variables 
a, j, x, 6, d, H, 3, I. 


The equations determine three of these variables if the 
others are known. Values of J may be assumed. The 
values of 3, \ and H which have been used are given in 
Table I. 

As for a, x, 5, and 7 we must regard one of them to be 
known. For some points plotted in Fig. 1 6 has been 
measured.* Thus, for these points the radius a and j 
and x can be calculated. For the largest currents, for 
instance, /=5X10* amps. we regard 4a\# as negligible 
and put x=1. Then the equations are able to deter- 
mine 6, but 7 and @ remain undetermined, since under 
this condition Eqs. (2) and (6) do not tell more than 
Eq. (5). 

In the case of carbon j is known,*‘ and a, 6 and x can 
be calculated. 


Ill. NUMERICAL RESULTS 


The results of the calculations are illustrated in 
Fig. 1. The number written above a plotted point 
gives j7; the number under the point gives x. 

All curves begin at small currents with a horizontal 
branch and end at very large currents in another hori- 


3 Some of these measurements are new and some are taken from 
E.C.H. Table (57.09) after conversion from cathode loss to cathode 
evaporation according to E.C.H. Eq. (57.03). 

“4 W. Finkelnburg, Hochstromkohlebogen (Springer-Verlag, Heidel- 
berg, 1948), p. 18; or “The high current carbon arc,” FIAT Final 
Report 1052, PB-81644, Office of Technical Services, Department 
of Commerce, Washington, D. C., 1948, p. 11. 
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zontal branch. It may be noticed that a horizontal line 
here indicates that the volume evaporated per second is 
proportional to the current J. 

The lower horizontal branch is an interesting singu- 
larity. It represents measurements cited in E.C.H.5 
Here only a small fraction of the energy which is carried 
by the ions to the cathode is used for vaporization. The 
metal atoms obtain the energy for their escape directly 
from the impinging positive ions. I would say, they are 
kicked out, and this effect is independent of the tem- 
perature of the cathode spot. In fact, it seems® that the 
fraction of the energy of the ions, which is used for 
the vaporization, or call it the disintegration or sputter- 
ing of the cathode, is nearly the same at the cool 
cathode spot of a glow discharge and the hot spot of an 
arc, at a small current. A. von Hippel’ has discussed 
this independence of the temperature as far as the glow 
discharge is concerned, but it is valid also for arcs. 

The kicking out of atoms will probably give the clue 
to the explanation of the high speed jets of metal vapor 
that have been observed at arc cathodes by Tanberg 
and Haynes.* The atoms in these jets have taken over 
more energy from the ions than is needed for the escape. 

The horizontal branch at large currents is a con- 
sequence of the condition that x= 1. Here the vaporiza- 
tion is primarily a normal thermal one. Since the term 
4ad8 in Eq. (5) is negligible, the conductivity \ does 
not affect this branch. The metal property which is 
decisive is H, the heat of vaporization, and it is not 
able to cause very great differences for different ma- 
terials. 

This high vaporization from the cathode is, of course, 
a nuisance for the constructors of switches, but it is 
inevitable if our assumptions are right. The vaporization 
from the anode and splashing of molten metal are other 
disturbing phenomena which shall not be discussed 
here. 

Finally, attention may be called to the high values of 
the current density 7 at small currents, which have been 
indicated in Fig. 1. They are consequences of the 
assumptions made. Recent measurements? indicate the 
existence of such current densities in transient arcs. 


5See E.C.H. paragraph 57, particularly Tables (57.13) and 
(57.14) and Eq. (57.03), also Appendices I and IV. 

® See reference 1, paragraph 57 E. 

7 The concerned papers by A. v. Hippel are cited and criticized 
by N. D. Morgulis, M. P. Bernardiner, and A. M. Patiocha. 
Physik. Zeits. Sowjetunion 9, 302 (1936); see also E.C.H. p. 316. 

8 See J. R. Haynes, Phys. Rev. 73, 891 (1948). 

®See H. Roth, Dissertation, Dresden 1939; L. H. Germer and 
collaborators, Lectures given at the gas discharge conference, 
Brookhaven National Laboratory, Oct. 28, 1948; J. D. Cobine 
and C. J. Gallacher, Phys. Rev. 74, 1524 (1948). 
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A Study of the e.m.f. Method* 


. C. T. Tat 
Cruft Laboratory, Harvard University, Cambridge, Massachusetts 


(Received September 15, 1948) 


The present work contains an investigation of the e.m.f. method to determine the input impedance of a 
thin biconical antenna as first suggested by Schelkunoff. By considering the various components of the cur- 
rents flowing on the antenna, it has been shown that the sinusoidal part is, indeed, the dominating compo- 
nent among all these currents. An investigation of the total current at the end of the antenna, where the 
lateral surface of the cone and the spherical cap meet, shows that the current is not identically zero at that 
point, but vanishes as the reciprocal of the square of the characteristic impedance of the cone. The two 
approximate expressions of the complex conjugate power used in the e.m.f. method are derived, as well as 
the exact expression. Computation of some integrals occurring in the formulation is also treated in detail to 


show the various approximations involved. 





HE e.m.f. method was first suggested by Bril- 
louin'? to compute the power radiated by an 
antenna, or by a system of antennas. It is an alternative 
to the standard Poynting vector method. The difference 
between the two is that the latter deals with the inte- 
gration of the Poynting vector over a big surface 
enclosing the system, while the former is based on 
direct integration of the Poynting vector along the 
conductor. As was shown by Bechmann,’ both methods 
are equivalent and one can be converted into the other 
by means of Poynting’s theorem. In this investigation, 
we will consider the applicability of the e.m.f. method 
to the determination not only of the power radiated 
by antenna but also its input impedance. This applica- 
tion was first studied by Schelkunoff‘ in his theory of 
biconical antennas, where he derived an expression for 
the effective load admittance of a thin biconical antenna 
based upon two apparently different methods, of which 
one is the e.m.f. method. A recent investigation’ of the 
same problem shows that Schelkunoff’s expression can 
be derived using a more rigorous method based upon 
Smith’s formulation of the problem,® and involving 
fewer assumptions than Schelkunoff’s original work. 
We have thus verified a result obtained by the e.m.f. 
method. What remains to be done, therefore and the 
object of the present work, is to give some interpreta- 
tions and justifications of the various assumptions in- 
volved in that method. 


* The research reported in this document was made possible 
through support extended Cruft Laboratory, Harvard University, 
jointly by the Navy Department (ONR) and the Signal Corps, 
U. S. Army, under Contract NSori-76, T. O. 1. 

!L. Brillouin, “Sur l’origine de la resistance de rayonnement,” 
Radioelectricite 3, 147 (1922). 

2A. A. Pistolkors, “The Radiation Resistance of Beam An- 
tennas,” Proc. I.R.E. 17, 562 (1929). 

3R. Bechmann, “On the Calculation of Radiation Resistance 
—- and Antenna Combinations,” Proc. I.R.E. 19, 1471 

1). 

*S. A. Schelkunoff, “Theory of Antennas of Arbitrary Size 
and Shape,” Proc. I.R.E. 29, 493 (1941). 

5 C. T. Tai, “On the Theory of Biconical Antennas,” Cruft Lab- 
oratory Technical Report No. 47, June 10, 1948; J. App. Phys. 
19, 1155 (1948). 

*P. D. P. Smith, “On the Conical Dipole of Wide Angle,” 
J. App. Phys. 19, 11 (1948). 
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‘FIELDS ASSOCIATED WITH A BICONICAL ANTENNA 


Without going into the details of the derivation, we 
will simply write down most of the equations derived 
previously.*~* The principal wave existing in region J 
of Fig. 1 is given by 


—jlo 
Es= fat KY ,)e~ #80 
4nr sind 
+(1-KY Jem], (ae 
—jlo 
Hy=——"| KY) 
4rr sind 


— (t= Ky jee» } (2) 


where J, as will be shown later, is the amplitude of the 
principal current on the antenna; Y;, is the effective 
admittance acting as a load upon the finite conical 
transmission line formed by the lateral surfaces of the 
cones; and K is the characteristic impedance of the 
cones, being given by 

















K=20/x (log cot@,/2); (3) 
where 
20> (u0/€o)*= 120x ohms. (4) 
The complementary waves existing in region J are 
J an Sn(Br) 
E,= — L,(8), (S) 
weer? » 2x S,(6l) 
P —jzo __ An 1 S,(Br) dL,(0) 
gn ‘S Qe n(n+1) Sa(Bl) 20 
—1 a, 1 S,(8r) dL,(0) 
i,-— ; (7) 


r ‘x Qn n(n+1) Sa(Bl) 90 


** The arbitrary constant V occurring in Eqs. (1) and (2) of 
reference 5 has been replaced by —j/o/4 in this paper, where J 
is a new constant having the dimension of current. The change is 
necessary in order to make our notation of the principal current 
the same as that of Schelkunoff. 
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Fic. 1. The biconical antenna. 


where 
S,.(8r) = (8r)*J n44(Br) (8) 
S,'(8r) =[dS,,(Br)/d(Br) }. (9) 
The complementary waves existing in region // are 
given by 
—j by, Ry(Br) 


E,=——- } ———P,(6), (10) 
weor k =1,3---2r R, (Bl) 





—jZo b, 1 R,. (Br) OP,(8) 
E,= — + —__——_— ———_ ——., (Il) 
r k=1,3---29 R(R+1) R,(Bl) 00 
a,=— — — (12) 

r k=1,3---29 R(R+1) R,(6l) 00 
where 

Ry (Br) = (8r)'H+4 (Br), (13) 
Ry’ (6r) = (dR, (6r)/d(Gr) }. (14) 


It was shown by Schelkunoff* that the problem of de- 
termining the input impedance is essentially that of 
determining Y;. Once this quantity is found the input 
impedance of the antenna is given by 





KY, sinBl—j cosél 
Z; -| lx. (15) 


sin6l—j KY, cospl 


By matching the H,- and E,-components of the in- 
terior and exterior fields across the boundary sphere 
and at the end-surfaces, the following set of equations 
can be derived :5 ® 


20 by; 1 
Woe (-) Pius), (16) 
wK? ¢=1,3--- \Io/ k(R+1) 


b, 
~jP.(us)= (—)w etre 1)P-(u1) 


0 





by P,.(11) Sn 
) g=1,3---, (17) 


, * 2 Se (r,n)(kyn). 
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where 











i= Cos6;, (18) 
_2n+1 * dn S,' (Bl) 
S,=j— (1—p)— —, (19) 
n(n+1) du; S,(6l) 
(r,.n)=r+r—n?—n, (20) 
: j Rl) 
W =—_ ——_. (21) 
2r+1 R,(sl) 


Equations (16) and (17) contain the formal solution 
for Y, for biconical antennas of any angle. For small- 
angle cones, it has been shown’ that Eq. (17) yields a 
rather simple solution for the coefficient (b,/Jo), which 
is given by . 


» WT (2) 
—= —j—(2r+1)81 4;(BD A 44(6)); (22) 
Te 2 


hence the corresponding value of Y; is 


zoGl 2k+1 @) 
Y.=— DY ——Jiy(6)Aeu(6l). (23) 
2K2x=1.3---k(k+1) 


For convenience, we may write 


K*Y ,=Z,, (24) 
with 
Z,.= RitjXa 
Z 2k+1 


“0 yA 


=-s1 © —Jyus3(81)Hiss(81). (25) 
2 «=1,3---k(R+1) 





The above formulation forms the basis of the ensuring 
discussion of the e.m.f. method. 


GENERAL EXPRESSIONS OF THE VOLTAGE 
AND THE CURRENT 


It is obvious that in order to apply the e.m.f. method 
to finding the power radiated from the antenna or to 
determining its input impedance, we need the expres- 
sions of both the voltage and the current acting upon 
the antenna. As shown by Schelkunoff,* if we define 
the voltage as the line integral of the E»-component 
along a meridian, then the interior complementary 
waves contribute nothing to the voltage thus defined, 
while the principal wave gives 


V(r)=KI,LKY;, sin@(l—r)—j cosB(I—r) ]. (26) 


The current flowing in the antenna can best be de- 
scribed by separating it into several parts. On the lateral 
surface of the antenna, the current associated with the 
principal wave is given by 


I(r) =IoLsin@(l—r) — 7 KY, cosB(/—r) ]. (27) 


Equation (27) is obtained by putting @=@, in Eq. (2) 
and then multiplying the resultant expression of H, by 
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2rrsin6;. The current associated with the interior 
complementary wave is obtained by using Eq. (7) 
instead of (2). That gives 


an S,(8r) dL ,(41) 
sin6; 
n n(n-+1) S,(6l) 00; 
an S,(8r) OL (U1) 


=> =p, 
n n(n+1) S,(Bl) Our 


I (r)= on 














In order to investigate the relative magnitudes of the 
currents defined by Eqs. (27) and (28), it is necessary 
to transform Eq. (28) into a form containing the co- 
efficients 5, instead of a,. It was shown by Smith® that 
as a result of matching the interior and exterior fields 
across the boundary sphere, the following relation 











exists: 
=f = ha byl nk (29) 
k=1,3--- 
where 
i—pr OL» (us) 
IT uae= P,(u1) " (30) 
k,n) Our 
1—p? OL n(u1) aL nas) 
lan = (31) 


2n+1 On Our 
By making use of the relation 


OL (m1) 
On dur Ou 


OL »(u1) dn 








(32) 


and Eqs. (29)-(31), there is no difficulty in changing 
Eq. (28) into the following form: 

















2n+1 1 
—— a k=1,3-- GAD (k,n) 
dn S,(8r) 
xi~- Kr a (gl) 
> 5 be fi ‘. 
n k=1,3- R(k+D Ln n—k n+1 
- -|a- ut) Pad) 56”) (33) 

n+k+ Sl) 


On the spherical end-surfaces of the cones, the current 
is equal to 2ml sindH,y, where Hy, is the magnetic field 
defined in region JJ and given by Eq. (12). Thus, we 
have 











1 dP, (8) 
I(o=- dX & sin@ 
k=1,3--- R(R+1) 06 
1 Pi(u) 
~ ———(1-)—, 34) 
k=1,3--- k(R+1) Ou 
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where u4=cos$. Equations (27), (33) and (34) are the 
general expressions of the currents flowing on the lateral 
surface and on the end-surfaces of a biconical antenna. 
Their specific values depend upon the amplitudes of the 
exterior complementary waves }, which satisfy an 
infinite set of linear equations as given by Eq. (17). 


CURRENTS ON SMALL-ANGLE CONES 


For small-angle cones we have shown that the value 
of 5, is given by Eq. (22) and that of Y; by Eqs. (24)- 
(25). The principal current defined by Eq. (27) can 
then be written as 


I(r) = Io sin8(—r) — j(Z./K)cosB(/—r)]. (35) 


‘For convenience, the following notations will be used: 


I,(r) =I sing(/—r) (36) 
I2(r)= — j(Zalo/K)cosB(l—r). (37) 


For J.(r), it has been shown® that when y:—>1, the 
following relations exist: 


dn 
lim (1—:?)——-= — (n—k)? 


mimi Mi 


1 0 
-- =-—, (38) 
[log(2/0,)? Ke 


lim (n—k) =[1/log(2/6:) ]=20/rK, (39) 


wil 





lim Pela) = ¥ (40) 


wil 


As a result, J-(r) can be split into two parts; the first 
corresponds to the sum of those terms of the double 
series that contain the factor 1/(n—k) for values of k 
nearly equal to n, while the remainder constitute the 


second part. Denoted respectively by I(r) and J;(r), 
they are 


_— a Ss 
i 2 ee. a 
WK k=1,3---k(k+1) S,(6l) 
Zo" b, 1 S,(8r) 


T4( )=—— oe 
ORR -R(k-+1) (n—B) Sa(B0) 








n=1,3-++k=1,3-- 
(n#¥k) 


by rs 
5. K+ Dn n+1 


1 " 
TR 
n+k+115,(8l) 








Substituting the value of ), as given by Eq. (22) into 
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(41), we have 


2olo 2k+1 
if)-j— 2 = 
2K &=1,3---R(R+1) 


| anuasian | @r 41.00% | (43) 


A similar expression can be derived for J,, which is then 
proportional to J)/K*. Thus, if J; is designated as the 
zeroth order distribution of current on the antenna, both 
I; and J; then have a weighting factor that is propor- 
tional to 1/K, while the weighting factor of J, is pro- 
portional to 1/K*. At the feeding point r=0, both J;(0) 
and J,(0) are identically zero. The input current is 
therefore contributed by /,(0) and /.(0) only. At r=, 
1,(2) is zero but J2(1), Js(/), and J,(/) are, in general, 
nonvanishing. From Eqs. (37), (43) and (25) it is seen 
that 


I.(l) = — j(Zalo/K), (44) 
and 
T3(1) = j(Zalo/K). (45) 


For J,(1), by putting r=/ in Eq. (42) and summing the 
series with respect to m, the following expression can‘ be 


2a—| fe 


| 
| 
he 


ie 
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(A) (B) 


Fic. 2. Small-angle biconical antenna showing 
the various surfaces. 
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derived: 
Ze by —jBlzcTo 
I.oO=— SY ——X(C+H(k+1)]J=— -— 
wK? «=1,3---k(kR+1) 24K? 
(2k+1)[C+y(k+1) ] 


ne i ol 
k(k+1) 


&=1,3°°: 


where C is the Euler’s constant and ¥(k+1) is the 
logarithmic derivative of '(k+1). Since the sum of 
I.(l) and J;(1) is zero, we conclude that the total axial 
current at r=/ vanishes as 1/K*, provided the series 
occurring in Eq. (46) is convergent. The latter state- 
ment will be postulated here pending a mathematical 
proof. In passing we may observe that in Schelkunoff’s 
method of partial matching,*® the vanishing of the 
total current at the end of the antenna is assumed as 
one intermediate step in his analysis. The actual picture 
now acquires a greater clarity from what has been dis- 
cussed here. Knowing the relative magnitudes of the 
various currents associated with the principal and 
complementary waves, we can proceed to a discussion 
of the parts taken by these currents in the e.m.f. 
method. 


POWER RELATIONS 


The essential idea in applying the e.m.f. method to 
find the input impedance of a small-angle biconical 
antenna is based upon a comparison of two different 
expressions of the complex power fed into the system. 
These power relations will be discussed under separate 
headings. 


A. Exact Expression for the Complex Conjugate 
Power Fed into the Antenna 


By means of Eqs. (15), (26) and (27) we find that the 
complex conjugate power fed into the antenna is given 
by 


P=3V(0)-1*(0)=4|1(0)|*Zin 


: = —K sin2£l 
=) hi] RG) —_ 





— 





2 sin2pl 
— K*B, cos26l+ K?* | Y,\? )} 47) 
2 
where G, and B, are respectively the conductance and 
susceptance of Y;; ie., 
Y ,=Git+ jBi. ‘ (48) 


As will be shown later, in making use of the e.m.f. 
method to find the input impedance of the antenna, 


instead of this exact expression it is the approximate — 


expression of the complex conjugate power which will be 
effective. 
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B. Approximate Expression for the Complex 
Conjugate Power Fed into the Antenna 


By approximate expression for the complex conjugate 
power we mean that expression of P in which an ap- 
proximate expression of J(0) is used but a more exact 
expression of Z;, is retained. Thus for 7(0) in 37(0)?Z;n 
we use Josinfl, but neglect the term —jJoKY; cosfl 
in Eq. (27); and for Z;, we expand the denominator of 
Eq. (15) using Taylor’s theorem and neglect the high- 
order terms. This gives 








if —K sin£l 7 
Ziun® - RG +i(— +B.) . (49) 
sinBlL 2 7 


Hence, 


1 —K singl : 
Papp=}|Io\?* KG+j(———+«8.) |. (50) 


L 








Comparing Eqs. (47) and (50) it will be noted that in 
addition to the term (K*|Y,|?sin26l)/2, which is 
negligibly small, the coefficients associated with B; are 
quite different for the two cases. This is of considerable 
consequence in determining Y, be means of the e.m.f. 
method, as will be seen later. 


C. Expression for the Complex Conjugate Power 
Derived from Poynting’s Theorem 


The foregoing expressions for the complex conjugate 
power are derived by regarding the antenna simply as 
a circuit element. On the other hand, an equivalent 
expression for the complex conjugate power can be de- 
rived by applying the Poynting theorem and integrat- 
ing the tetal power along the surface of the conductor. 
Referring to Fig. 2(A), we observe that the total com- 
plex conjugate power radiated from the antenna is 
given by 


p=3f E- H*do (51) 
71,092,093 


where a; is an infinitesimal spherical surface surround- 
ing the generator, o2 is the lateral surface of the cones 
and o3 their top surface. If we use the exact expressions 
of E and H to compute P, then, since the tangential 
component of E. is zero along o2 and o3, the only con- 
tribution is due to the integration along o:. It can be 
shown without any difficulty that by using (1) and (2) 
an expression for P is derived that is identically the 
same as Eq. (47). There is, however, another approach 
to the problem that enables us to find an approximate 
expression for the complex conjugate power based upon 
Poynting’s theorem. 

Referring back to Eq. (51), let us decompose the 
integral into several parts and express them in different 
forms. Thus if we denote 


Pint f E-H*de=3 f E,-H4*do (52) 
71,02 71,92 
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and 
Pmt f E-Htdo=3 f Ey: Hg*deo, (53) 
73 o3 
then 


P=P,+P,. (54) 
For P;, one has 


l 
P\= sf Ev Hstdo=) { 2rr sin6,E,H 4*dr 
1,02 -I 


l 
= f L,*Edr, (55) 
-—l 


. where J, is the total current flowing along the lateral 


surface of the cone. For P;, one has 
1 a1 
Pi=-f E.H,*do= f 2rl* sind EH 4*dé 
2S «3 0 
6 
-1f Ip*Epd@ (56) 
0 


where J¢ is the total radial current flowing on the top 
surface of the cone. It will be recalled that for small- 
angle biconical antennas, J, can be separated into four 
parts given by Eqs. (36), (37), (41) and (42), while J» 
is given by Eq. (34). If we think of each of these cur- 
rents as a source of radiation, then there is an electro- 
magnetic field associated with each of these currents. 
As a result, we may write 


5 
E,=>, Ens (57) 
n=l 
and 
5 
Es=>, Ens, (58) 
n=1 


where, for example, Ei, is the r-component of the E- 
field, expressed in a spherical coordinate system, pro- 
duced by J;; and E;,, E59 are the components due to 
Is, which will henceforth be denoted by Js. Substi- 
tuting these formal expressions of J ,, Js, E, and Eg into 
Eqs. (55) and (56), we obtain 


tlepya4 5 
Pn > I,* y Ear) (59) 
2 ill k=1 


n=1 


1 
P=! f (1 > Eas (60) 
0 n=1 


Of all the cross-products involved in Eqs. (59) and 
(60), because of the different orders of magnitude of 
the various currents, it is obvious that the most sig- 
nificant part is due to the integration of /,*E,,. The 
remaining terms, though contributing only negligibly 
to the value of the total complex conjugate power 
radiated from the system, nevertheless are important 


and 
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in the sense that by their inclusion the boundary condi- 
tion is automatically satisfied. 
The evaluation of the integral 


1 I 
Pope= = f I,*E,dr (61) 
=! 


was briefly described by Schelkunoff both in his paper* 
and in his book,’ but since the computation definitely 
involves some approximations, it is perhaps worthy 
of a more detailed derivation here. Moreover, it will 
be seen that the steps adopted here do not precisely 
follow those outlined by Schelkunoff. It is well known*® 
that the field produced by a sinusoidal distribution of 
current 














is given by 
jl veo — gp iBR1 e iBR2 2 
“= _ +— cos$le- on| (63) 
4 Ri R, Ro 
ploeofz—l 4 2z cospl 
1 p= ——| ——e7 #14. —¢~ ib Ra cen | (64) 
4p R, R2 Ro 
fof ) : 
Hy= le ipRi- ¢~ i6R2— 7 conse (65) 
4p 
where 
Ri=((I—2)*+ 6") 
R2= ((l+-2)?+p")! (66) 
Ro=(3-+6")! 


and p is the distance measured from the point of ob- 
servation to the z-axis. According to Schelkunoff,’ one 
can use £,,+£,, sin; for Ei, in Eq. (61) to perform the 
integration. The resulting integral, however, is some- 
what complicated. A simpler method is to perform the 
integration along a cylindrical surface as shown in 
Fig. 2(B), and then to make the necessary correction 
for a conical surface by means of Poynting’s theorem. 
Denoting the lateral surface of the cone by o;, the sur- 
face of the open cylinder by oo, and the volume in 
between these two surfaces by 7, we find? that 


1 1 w 
-{ E-H*de=— { E-H*de+ j- { (uit eF*)dr (67) 
2 oi 2 70 2 T 

or 


1 7’ 1 ¢' 
-{ Is*Exdr=— f I,*E\,dz 
2/_; 2/_; 


@ a l ; 
+j- f 2ep(uH?—eE2)dpdz. (68) 
2 (a/l)z*% —1 


7S. A. Schelkunoff, Electromagnetic Waves (D. Van Nostrand 
Company, Inc., New York, 1943), p. 453. 

5 See reference 7, p. 371. 

*R. W. P. King, Electromagnetic Engineering (McGraw-Hill 
Book-Company, Inc., New York, 1945), Vol. I, p. 221. 
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The two integrals at the right-hand side of Eq. (68) can 
be evaluated as shown in the appendix, giving 


roe | To|2(R+7X) 69 
| Pechaeatita ial (69) 


where 
R= 20/4n[_2(C+ log28l—Ci261)+ (Si48]—2Si2Bl)sin2pl 
+ (C+ log6l— 2Ci2B14+-Ci4fl)cos261] (70) 


X = 20/4e[2Si261—Si4Bl cos26l+ (—C—loggl 
—2 log26l+2 log8a+Ci4gl)sin2sl]. (71) 


For small-angle cones, as treated here, the character- 
istic impedance of the cones, is given by 


K = (20/1)log(2/01) = (zo/m)log(26l/Ba). (72) 
Thus, X can be written as 


X=(—K sin26l/2)+ 20/4x[2Si261— Si4Bl cos2pl 
+(—C—log8l+-Ci46l)sin2g1]. (73) 


Now we have two approximate expressions of the power 
fed into a thin biconical antenna; one is given formally 
by Eq. (50), while the other has been derived from 
Poynting’s theorem and is given by Eqs. (69), (70) and 
(73). If they are equivalent we should obtain 


K°G,= 20/4x[2(C+log26l—Ci28l) 
+ ($i481— 2Si26l)sin2sl 
+ (C+logBl—2Ci2BI+-Ci4Bl)cos261] (74) 


and 


K*B,= 20/4n[2Si261— Si4gl cos26l 
— (C+ log6l—Ci4sl)sin2sg1}. (75) 


As shown by Rice,’® the expression of K°?Y, thus ob- 
tained is identically equal to the real and imaginary 
part of Z, given by Eq. (25), which was obtained more 
rigorously without resorting to the various approxima- 
tions encountered in the e.m.f. method. It is to be no- 
ticed that although we use the approximate expressions 
of the complex conjugate power to determine the value 
of K*Y,, the actual power fed to the system is given by 
Eq. (47) instead of by Eq. (50) or Eq. (69). On the 
other hand, if we identify Eq. (47) with Eq. (69) the 
value of KB, thus determined would not be the same 
as before and it would not be the correct one, even 
though the value of K°G; is unaffected. This is an in- 
dication of the critical nature of the e.m.f. method: 
Unless it is handled with extreme care an error may 
easily result. 


APPENDIX 


1 l 
Evaluation of the Integral > f I,* E,,dr 
—2 


0S. O. Rice, “Sum of Series of the Form 20”A »J n4a(2)Jn+a(2),” 
Phil. Mag. XXXV, 686 (1944). 
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The first integral at the right-hand side of Eq. (68) is known*** 
to have the following value, for a@&P. 


sf it Bve= BL Lol(R+IX0) (A-1) 
where 
R= rm [2(C+log2¢l— Ci2Bl) + (Si4gl—2Si2Bl)sin26l 
+(C+logsl—2Ci2g1+Ci4pl)cos2pl] (A-2) 
X= q_L2Si2pl+ (2Si261— Si4pl)cos2,l 
+(C—loggl+2 log8a—2Ci2g1+Ci4pl)sin26l]. (A-3) 


To evaluate the second integral of Eq. (68), we use the expressions 
of E and H given by Eqs. (63)—(65). It can be shown that 
pl? [ 1 ena Re? 


—P 
pH? —«E? ~ Baty? ee eos(R:—R) 


R?—a-l 
_ 21 _ if )eossl cos6(Ri— Ro) 


0 





Re+2l 
-A1- 2 cos coss(R:— Ro) |. (A-4) 


When pz, (/—z), (/+2), the following approximations are 


*** S. A. Schelkunoff, Electromagnetic Waves (D. Van Nostrand 
Company, Fnc., New York, 1943), pp. 372-373. 





permissible : 
R&P . 2 
i--R. R.— R,=22 
R?—al 
1— =2, Ri—Ro=l—2z 
RoR: 1 0 
Rg +2l 
———=0, R.—Ro=l. 
RoR2 2 0 
Then Eq. (A-4) reduces to 
pH? — cE? = (—p| Io|?/42p?)cos2B(I—2z). (A-5) 


Substituting Eq. (A-5) into the second integral of Eq. (68) and 
performing the integration, one obtains 


@ fa l . 
Fe yad 2x0 ul — eE*)dpdz =} |Io|?AX (A-6) 


where 


AX= — 4,2 [2i2el cos26l+-2(C— Ci2gl+-log2sl)sin2sl]. (A-7) 
Thus, we have 
i * Eudr=}|Io|2(R+jX) (A-8) 
aJ-1 1 £1 0 j 
where R is given by Eq. (A-2) and 
X=X.+AX= jlas 12p81— SiApl cos2pl 
T 


+ (—C—logsl—2 log26l+-2 log8a+Ci4gl)sin2p1]. (A-9 





BOOK 


Atomic Theory for Students of Metallurgy, 
Institute of Metals Monograph No. 3 


By Wri1t1AM Hume-Rortuery. Pp. 282, Figs. 124, 7s.6d. The 
Institute of Metals, London, 1946. 


This monograph presents a summary of the mathematical 
theories of the electronic structure of atoms, molecules, crystals 
and alloy systems in language understandable to the advanced, 
non-mathematically minded student. The accomplishments of 
wave mechanics as applied to the fundamental questions that ex- 
plain the structure of solids, or that show promise of ultimately 
doing so are presented with remarkable clarity. The author has an 
unusual gift for this kind of writing, and has a breadth of first 
hand information on the nature of alloys that comes from years of 
research directed to understanding the fundamental principles 
governing the structure of alloys. Although not a mathematical 
physicist himself, he has studied this field extensively and has been 
guided past many pitfalls in drawing conclusions from the pub- 
lished theories by discussions with Professor Mott and Dr. Fro- 
lich, as he says in a footnote. 

The book outlines many of the metallurgical problems, the 


VOLUME 20, JULY, 1949 


REVIEW 


successes and some of the inadequacies of theory in a way most 
helpful to the metallurgist and the average graduate student in 
the sciences; it does not attempt to derive any of the principles, 
but merely to explain, interpret, apply, and test them. More 
attention is given to the successes than to the failures of theories, 
as is understandable. 

There is no extensive overlapping of the contents with the con- 
tents of Hume-Rothery’s “Structure of Metals and Alloys,” an- 
other monograph of the same series (revised 1947). The subjects 
treated include the Bohr theory, the wave equations and their 
interpretation, potential boundaries in wave mechanics, the elec- 
tronic structure of the atom and assemblies of atoms, the inter- 
actions between atoms in molecules and crystals, Fermi-Dirac 
statistics and the electron gas, the Bruillion zone theory, and the 
application of the theories to the metals, particularly copper, 
silver, gold, the hexagonal metals, and the elements of the first 
long period. The fields of elasticity, anelasticity, plastic deforma- 
tion, and preferred orientations are not covered. 

C. S. BARRETT 
Institute for the Study of Metals, 
The University of Chicago 
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Letters to the Editor 








Ionospheric Absorption 
R. N. BRACEWELL AND K. WEEKES 
Cavendish Laboratory, Cambridge, England 
April 11, 1949 


OUBT has been cast, by your correspondent L. G. 
McCracken,' on the formula given by Best and Ratcliffe? 
for the absorption of a radio wave which penetrates an ionized 
region of the Chapman type. He states that an approximation has 
been made, as a result of which the formula is correct only when 
the sun is near the horizon. He states that when the sun is over- 
head, the formula is wrong by a factor of approximately one 
quarter. The formula is, however, well-established and is in agree- 
ment with the results obtained in a different way by Appleton.’ 
The formula given by Best and Ratcliffe, under the assumption 
that the wave frequency is much greater than the collision fre- 
quency, comes out in terms of the error function as follows: 


H  4névo (2) L ’ 
2mc (p+ | pr\) C secx 


x [) exp(— u?) — f exp(—w')du] (1) 


§C secx Le 


—} loge= 


McCracken obtains the following formula under the same as- 


sumption: 
H 4rre?v9 (7) 
ne C secx)! 
2mc (p+| pr!) G, 4 ~~ 


where (x, y)! is the incomplete factorial function. Thus where the 
first formula has 


—} loge= 


(4C secx)* 
f exp(— u*)du—(4C secx)! exp(— $C secx) (2) 


the second has 
(4, 4C secx)! 


Now it is a property of the incomplete factorial function that 
(4, a) !=4(—4, a)!—ate. 
Further, since 


(n'—1, x") != nfo exp(—!")dt, 
it follows that 


i(—4, 0) != f™ exp(—)adt. 


Hence: 
(4, a) =f exp(—?)dt—ate. 


The expression (3) is therefore equal to (2) if we identify a with 
4C secx, and McCraken’s formula is identical with that of Best 
and Ratcliffe. 

The numerical discrepancy mentioned above arose because 
McCracken erroneously took (4, $C secx)! equal to (4, 4)! when 
x=0, thus implying that C=1. Now (4, 4)! is smaller than the 
correct value by a factor 0.199 which gives rise to the discrepancy 
stated to be approximately one quarter. 

Best and Ratcliffe do not, after setting up their integral, 
“immediately go to an approximation,” but correctly perform the 
integration, obtaining an expression in the form (1), and this form 
is simpler than that of McCracken. To evaluate their expression 
they discard the second term in (2) and take the upper limit of the 
error integral to be infinite. The change introduced is small— 
about ¢~78-#, 


These adjustments were numerical, but have the following 
physical significance. The variable u in (1) comes from the sub- 
stitution «= ($C secx)! exp(— 4$4/H), where h is the height above 
the ground, H is the scale height of the atmosphere, x is the sun’s 
zenith distance, and C secx=exp(/m/H), where hp is the height 
of maximum ionization density. The limits of integration, 
(4C secx)? and 0, in (1), signify that the integral is taken from the 
ground level to infinite height. If the integration is taken from 
negatively infinite height upwards to infinity, as throughout 
Appleton’s treatment, the absorption varies as (cosx)!. But in the 
general case of integration upwards from some finite level, the 
rising and falling of the region, as the sun’s zenith distance changes, 
leads to a departure from the (cosx)! law. In the practical case of 
the ionospheric layers, the height of maximum ionization for 
region E is about twelve scale heights, and C secx thus has the 
large value e”. The result is almost identical with that obtained by 
assuming that the integral extends over all values of height, and 
this is what Best and Ratcliffe get in replacing (2) by 42%. 

Appleton’s treatment makes the same assumption, in that he 
takes the integral over all values of height. He also evaluates the 
integral numerically for any relative value of wave frequency and 
collision frequency. 

It may be noticed that the case of C=1, to which McCracken’s 
conclusion applies, corresponds to a layer whose maximum ioniza- 
tion density is at ground level when the sun is overhead, and it is 
to be expected that such a case would lead to a severe departure 
from the (cosx)! relation. 

IL, &. McCracken, J. App. Phys. 20, 229 (1949). 


2j. Best and J. A. Ratcliffe, Proc. Phys. Soc. London 50, 233 (1938). 
3E, V. Appleton, Proc. Roy. Soc. A162, 451 (1937). 





Errata: On Theoretical Signal-to-Noise Ratios in 
F-M Receivers: A Comparison with 
Amplitude Modulation 
[J. App. Phys. 20, 334 (1949)] 

Davip MIDDLETON 
Cruft Laboratory, Harvard University, Cambridge, Massachusetts 


Equation (4.2): 6 should read @. 

Equation (5.20): the left number is (s/m)?77;-1in. Omit the 
word “‘where” immediately preceding Eq. (II.2). 

Equation (II.7) : 2 should read x*. In the sentence immediately 
preceding Eq. (I1.8), @ and y should be replaced by x and », 
respectively. 

Equation (I1.10): replace ¥ by v. Write (B; (R)6) in the sen- 
tence following Eq. (II.10). 

Equation (II.13): Replace ¢ by x, and in the equation imme- 
diately following write » for ¥; sin(é—y) should read sin(x—»). 
In the first sentence following Eq. (II.31) we should have 2¢? vs. 

p?/ac, instead of 2G? vs. p?/ac'. 





Erratum: The Thermal Behavior of Evaporated 
Layers in Vacuum Devices 
{J. App. Phys. 20, 8 (1949)] 


H. A. STAHL 
Laboratory for Material Physics, Telefunken Tube Factory, 
Berlin, Germany 


HE author has pointed to an uncompleted reference. On 
page 9, the following should be added to footnote 2: L 
Mathes, Zeits. f. tech. Physik 22, 232 (1941). 
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